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I.  INTRODUCTION 


The  contract  N00019-78-C-0064  entitled  "An  Investigation  on  Characterizing 
Mutual  Coupling  Between  Two  Antenna  Slots  on  a Cone"  was  awarded  to  the 
University  of  Illinois  by  Naval  Air  Systems  Command  for  the  period  of 
16  November  1977  to  15  November  1978  and  with  funding  of  $53,000.00.  The 
contract  was  later  extended  to  15  January  1979  with  no  additional  cost. 

The  contract  monitor  is  Mr.  J.  Willis  of  AIR-310B. 

This  is  the  final  report  of  the  contract,  covering  personnel  (Section  II), 
technical  results  (Section  III  and  attachments)  and  publications  (Section  IV). 
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R.  Hittra,  Professor  of  Electrical  Engineering,  Principal  investigator 
E.  Yung,  Research  Associate 
L.  Grun,  Research  Assistant 
C.  L.  Law,  Research  Assistant 


III.  TECHNICAL  RESULTS 


Our  study  of  the  mutual  coupling  effect  in  a conformal  array  has  been 
successfully  concluded.  In  the  present  contract,  the  following  two  tasks 
have  been  carried  out: 

(i)  GTD  solution  of  self-admittance  of  slot  on  a cone  or  cylinder. 

In  applying  GTD  formulas  to  calculate  the  self-admittance  of  a slot  on  a 
cone  (cylinder),  there  is  a difficulty  which  was  not  previously  presented 
in  the  calculation  of  mutual  admittance,  namely,  the  GTD  Greens'  function  G 
for  the  surface  field  has  a r ^-singularity  at  the  source  point  r * 0.  IVe 
removed  this  difficulty  by  subtracting  from  G,  where  G^  is  the  Green's 
function  for  a planar  conducting  surface.  The  contribution  of  G^  to  the 
self-admittance  can  be  calculated  by  a Fourier  transform  method.  The 
remaining  function  (G  - G ) is  of  order  r 1";)  as  r - 0,  and  is  therefore 
integrable.  When  the  slot  is  on  a cylinder,  our  GTD  results  of  the  slot 
self-admittance  are  in  excellent  agreement  with  those  calculated  from  the 
known  exact  solution.  Details  are  given  in  Attachment  A. 

(ii)  Justif icat ion  of  the  transverse  curvature  term  in  the  GTD  solution. 
In  our  GTD  solution  for  the  magnetic  field  on  a convex  conducting 

surface,  there  exists  a rather  peculiar  term.  Contrary  to  all  previous 
GTD  theories,  this  term  depends  on  the  surface  curvature  in  the  transverse 
direction  of  the  ray.  For  many  practical  situations,  the  inclusion  of  this 
term  is  of  critical  importance  in  getting  accurate  numerical  solutions. 

This  term  was  first  introduced  by  us  in  November  1976  as  a conjecture.  Now, 
we  have  shown  through  a rigorous  asvmptotic  expansion  of  an  exact  solution 
that  our  conjecture  is  indeed  correct.  Details  are  given  in  Attachment  B. 
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Our  work  on  the  GTD  calculation  of  mutual  coupLing  is  summarized  in 
a review  article  (Attachment  C),  which  will  be  included  in  the  forthcoming 
book  entitled  Principles  and  Applications  of  Antenna  Design  sponsored  by 
IEE  (London). 
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S.  W.  Lee,  E.  Yung,  and  R.  Mittra 
Department  of  Electrical  Engineering 
University  of  Illinois  at  Urbana-Champaign,  1978 


ABSTRACT 

The  input  admittance  of  an  elemental  radiator  on  a curved  surface, 
e.g.,  a slot  on  a cone,  plays  an  important  role  in  the  design  of  conformal 
arrays.  A search  through  the  literature  reveals  that  at  present,  there  is 
no  reliable  theoretical  method  available  for  computing  this  admittance. 

The  objective  of  this  paper  is  to  provide  a solution  to  this  problem  using 
a surface  ray  approach  — within  the  framework  of  GTD.  The  solution  is 
verified  for  the  limiting  case  where  the  cone  degenerates  into  a cylinder 
and  it  is  shown  that  the  GTD  results  compare  extremely  well  with  the  exact 
modal  solution  to  the  cylinder  problem.  Extensive  numerical  results  are 
presented  in  the  paper  for  the  input  admittance  of  a cone  as  a function  of 
varioin? design  parameters. 


*This  work  was  supported  by  Naval  Air  Systems  Command  under 
Contract  N00019-78-C-0064 . 
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1 . INTRODUCTION 

Because  of  its  simplicity  of  geometry  and  ease  of  flush-mounting, 
the  slot  radiator  is  one  of  the  most  frequently  used  elemental  radiators 
in  the  design  of  conformal  arrays.  A crucial  design  parameter  of  a slot 
is  its  input  admittance  Y,  whose  value  depends  on  the  slot  dimensions  and 
the  geometrical  property  of  the  conducting  surface  on  which  it  is  mounted. 

In  the  simplest  case,  the  mounting  surface  is  an  infinite  ground  plane,  in 

which  Y of  a thin  slot  can  be  related  through  the  duality  relation  to  the 

input  impedance  of  a thin-wire  antenna  in  free  space.  The  latter  quantity 

was  first  calculated  by  P.  S.  Carter  in  1932  [1]  - [3].  Direct 

calculations  of  Y of  a slot  on  a plane  were  reported  in  [4].  When  the 

mounting  surface  is  a cylinder,  the  solution  of  Y has  been  expressed 

exactly  in  terms  of  cylindrical  modes,  namely,  an  infinite  series  in 

the  azimuthal  direction  and  a spectral  integral  in  the  axial  direction  [51,  [6]. 

This  modal  solution  is  suitable  when  the  radius  R of  the  cylinder  is 

small  in  terms  of  wavelength  (kR  < 10);  otherwise, its  numerical  evaluation 

is  extremely  laborious. 

In  the  present  paper,  we  consider  the  calculation  of  Y of  a slot 
on  a cone  (or  cylinder)  using  surface  rays  in  GTD.  The  general  concept 
of  surface  rays  was  introduced  by  J.  B.  Keller  in  1956  [7].  The 
explicit  formulas  for  surface  rays  adopted  here  are  those  reported 
recently  in  [8].  Our  solution  of  Y is  an  asymptotic  solution  which  is 
valid  when  the  radii  of  curvature  at  all  points  on  the  cone  are  large 
in  terms  of  wavelength  (or  high  frequency  solution) . Its  calculation 
is  relatively  simple, and  its  accuracy  is  surprisingly  good  as  verified 
by  the  comparison  with  the  exact  modal  solution  for  the  case  of  3 slot 
on  a cylinder. 
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FORMULATION  OF  INPUT  .ADMITTANCE 


Referring  to  Figure  Id,  let  us  consider  a slot  on  the  surface  of 
an  infinitely  large  conducting  cone  with  a half-cone-angle  0^.  The  slot 
is  relatively  small  when  compared  with  its  surrounding  cone  surface, 
and  the  shape  of  the  slot  is  assumed  to  be  rectangular  on  the  deve loped 
cone  (Figure  lb).  Note  that,  depending  on  the  exact  manner  in  which  the 
feeding  waveguide  is  fitted  into  the  cone  surface,  the  shape  of  the  slot 
mapped  onto  a developed  cone  can  be  quite  irregular.  The  assumption  of 
rectangular  shapes  represents  a good  approximat ion  for  practical  cases; 
at  the  same  time,  it  simplifies  the  subsequent  calculations.  The 
location  and  dimensions  of  the  slot  are  described  bv 
(c,  m) , and  (a  x b)  , 

where  c is  the  radial  distance  of  the  center  of  the  slot  from  the  cone 
tip,  and  uj  is  the  angular  deviation  of  the  slot  axis  from  the  cone 
generator  (co  - tt/2  for  a circumferential  slot  and  u)  * 0 for  a radial 
slot).  Throughout  this  work,  we  assume  that 

(i)  kc  >>  1 (slot  awav  from  cone  tip)  , 

(ii)  ka  - tt  ar.d  kb  <<  1 (resonant  thin  slot) 

Assumption  (i)  is  necessary  because  near  the  tip  the  cone  is  highly 
curvec  and  GTD  is  not  valid  there.  As  a consequence  of  assumption  (ii), 
the  aperture  field  of  the  slot  can  be  adequately  approximated  bv  a 
simple  cosine  distribution,  i.e.,  the  so-called  "one-mode  approximation": 

(2-la) 


E » z V / cos  (— *-)  , 

/ ab  a 


, TV. 


H * — v I / -t-  cos  (— — ) 
► ab  a 


( 2-lb) 


where  V 4ml  I are,  respectively,  the  modal  voltage  and  current  of  the 
slot.  The  Input  admittance  of  the  slot  is  defined  by 


T 


Alternatlvely.lt  can  be  calculated  bv  the  expression 


Y 


K ds 


(2.3) 


Here  A Is  the  aperture  la  x b)  of  the  slot.  K is  the  equivalent  surface 
magnetic  current  density  due  to  the  application  of  a field  E described 
In  v2.1a),  and  is  given  by 


K-Exn*Exx.  (2.  a) 

The  magnetic  field  H in  (2.3)  is  the  field  produced  by  K when  the  slot 

— ► 

Is  short-circuited.  Under  the  one-mode  approximat ion , this  H is  identified 
with  the  radiation  of  K located  on  a completely  filled  cone  (.without  slot). 
We  emphasize  that  (2.3)  is  an  approximat ion  of  i2.2)  because  of  the  manner 
in  which  H is  calculated,  and  is  valid  only  under  the  one-mode  approximat ion 
in  (2.1). 
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3.  CTD  SOLUTION  OF  INPUT  ADMITTANCE 

To  calculate  H in  (2.3),  we  first  determine  the  Green's  function 

h(I,2)  which  represents  the  magnetic  field  at  point  2 due  to  a magnetic 

dipole  at  point  1.  Both  points  1 and  2 are  on  the  surface  of  the 

completely  filled  cone  (after  the  slot  is  removed),  as  shown  in  Figure  2. 

According  to  GTD,  h ( 1 , 2 ) has  the  following  two  dominant  contributions 

-*d 

at  high-frequency:  h due  to  the  direct  ray  12  going  from  the  source 

-*t 

to  the  observation  point,  and  h due  to  the  tip-diffracted  ray  102. 

Thus , 


-*  ‘■d  "*t 

h(l ,2)  ~ h + h ; k - * . 

Accordingly,  input  admittance  Y also  has  two  parts 

Y - Yd  + YC  ; k - *> 


(3.1) 


(3.2) 


The  calculation  of  hd  is  detailed  in  [8],  and  that  of  hC  in  [5],  [8]. 
Applying  the  above  results  in  (2.3),  Yd  and  YC  on  a cone  are  found  as  follows 
(for  expjuit  time  convention): 

Direct  contribution: 

, , f a/2  fb/2  fa/2  fb/2 

Y “ ‘ ^Tj  dyl  J d2l  J dy2  dz2 

-a/2  -b/2  -a/2  -b/2 


(3.3) 


ra/2 

rb/2 

W 2 

J 

dyi 

J dzL 

-a/2 

-b/2 

-a/2 

a 

cos 

*y<> 

(~T)  8(yr  2 

where 


g(y1,  z^;  y^,  z,)  - sin"  9 + Hc  cos'  9 

% " G(s)  ( 1 • i Tv(°  • (iij:  r3u(° 


+ j(r'2kRc)  [tv’(0  + (R  /I^)  T 


tV(dA  , 


(3. A) 
(3.5) 


Ht  - G(s) 


(i)  Tv(5)  + 1 " ^ T3u<£)  + j (<^2kR~)_“/3  t3u’(0 

L 


,2  - j ks 

k e 
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S * /(y2  " yl)2  + (z2  " Zl)2 


, (3.6) 


(3.7) 


(3.8) 


0 ■ tan 


■1  Z 2 ~ 21 
y->  - yi 


(3.9) 


> ■) 

ft  » /c“  + s“  - 2cs  COS  (u)  - !jl)  ) , 


n n 


(3.10) 


* i>  sin  9 * sin  \ ^ sin  (u  - u)  ) 

n n 0 \ SJ  n 


(3.11) 


3 • f\  + Z 

n ' n n 


(3.12) 


<a  a tan  (z  /v  ) , 

n n n 


(3.13) 


sin  |$'  - 4>*  I 
T - — ; vr  " ;<t  1 

1 ^2  “ ^l' 


(3.14) 


kr^  sin  ft^  1/3 
2 tan~  0„ 


i*;  - 


(3.15) 


- ^ >rlr2  tan  e0 
t ” sin  ft^  sin  ft0 


(3. in) 


'rlr2  tan  9Q 

cos  ft^  cos  ft? 


(3.17) 


Fock's  functions  (u,v)  and  their  derivatives  are  described  in  14] . 


Tip  contribution: 


..t  _ .2  (1  + Dab  f Can  ^0 

Y >a0sin  ."T  [ tt 

60*t  c sin  0g  v 


sin  (kb/ 2)  ~ -j 3kc  . 

kb/2  • (JllS) 
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Here  is  the  ieroth-order  tip-dif fraction  coefficient  and  is  a 
function  of  the  half-cone-angle  It  is  approximately  given  by  [8] 


J-  » Ae 
0 


JB 


(3.19) 


where 


A - 1.3O570*1 
0 


1.753  2 . 77200  - 1.4590“ 


B « 2.7195  + 1.46O80Q  - 1.12956“  + 0.656600 


in  which  both  00  and  B are  in  radians. 


(3.20) 

(3.21) 


If  the  slot  is  on  a cylinder  and  is  oriented  along  the  circumferential 
direction  (rig-ore  3),  the  same  expression  in  (3.3)  can  be  used  to 
calculate  the  input  admittance  Y (there  is  no  tip  contribution  for  the 
cylinder  case).  The  GTD  solution  of  the  Green’s  function,  which  is 
determined  in  [4],  is  given  by 


sty,.  *1:  y , . 2,)  - G(s)  (v(S)[sin“  4 ^ 201 


uu)  [cos'  J ( 1 - f^-)  + A sin“  0] 
ks  ks  ks 


+ r^-  s [ v ’ ( £)  sin“  0 + u’(^)  cos'  0 (tan"*  0 + r^-)  1 } , 

ks  ks 


(3.22) 


,0  - 1 1/3 

kx  cos  r I s 


where 


R * radius  of  the  cylinder. 

The  functions  G(s),  s,  0 are  defined  earlier  in  (3.7)  through  (3.9). 

In  summary,  the  GTD  solution  of  the  input  admittance  of  a slot  on 
a cone  is  given  in  (3.2),  (3.3),  and  (3.18);  while  that  of  a cylinder 
is  given  in  (3.3)  and  (3.22).  These  solutions  are  approximately  valid 


( *.23) 
( 3.24) 
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4.  FINITE  PART  OF  DIVERGENT  INTEGRAL 
The  integral  for  calculating  Yd  in  (3.3)  is  in  fact  a divergent 
integral.  This  is  due  to  the  fact  that,  as  point  1 approaches  point  2, 
(Figure  2)  the  Green's  function  in  (3.4)  becomes  infinite  as 


g(yL.  zL;  y-,.  z2)  ~ cs 


s 0 , 


(4.1a) 


where  s,  defined  in  (3.8),  is  the  distance  between  the  two  points,  and 


the  parameter  C is 


(2-3  sin3  8) 


(4.1b) 
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It  is  well-known  that  the  singularity  of  cubic  power  is  non-integrable 
with  respect  to  a surface  integral.  This  difficulty  can  be  traced  back 
to  the  derivation  of  the  Green's  function  g.  Strictly  speaking,  g is 
a distribution  and  can  be  written  as 


g = Dg  , (4.2) 

where  D is  a second-order  differential  operator  with  respect  to 
coordinates  of  point  2,  and  g is  the  Green's  function  of  a vector 
potential  component.  A "legitimate"  expression  corresponding  to  (3.3) 


should  read 


Yd  = - -k  If  dy2  dz2 cos  (~r)  {D[  ||  dyi  dzi cos  (~r)g]}  * (4-3) 

which  is  convergent,  and  Yd  has  a well-defined  finite  value.  However, 
in  writing  (3.3),  we  have  interchanged  the  differential  operator  D and  the 
second  surface  integration  operator  in  (4.3).  This  interchange  is  not 
permissible  and,  therefore,  leads  to  the  divergent  integral  in  (3.3). 

Since  (4.3)  contains  a differential  operator  and  is  net  suitable 
for  numerical  evaluation,  we  prefer  to  work  with  (3.3),  provided  of  course 


< 


AT 


we  can  extract  the  correct  finite  part  from  the  divergent  integral.  To 
this  end,  we  rewrite  the  Green's  function  in  (3.4)  as 


g0  + gl 


(4.4) 


The  first  term  g^  in  (4.4)  is  the  Green's  function  of  an  infinite  ground 
plane,  and  is  given  by  the  well-known  expression 


gQ  = G(s)  [sin2  6 + ^.(1  - •j^-)(2  - 3 sin2  0)  ] 


(4.5) 


Note  that,  as  s -*•  0,  g^  has  exactly  the  same  singular  behavior  in  (4.1) 
as  g.  This  is  expected,  because  in  the  sufficiently  small  neighborhood 
of  a point  source,  the  cone  can  be  approximated  by  its  tangent  plane. 
The  second  term  g^  (g^  = g - g^)  in  (4.4)  is  the  difference  between  the 
Green's  function  of  a cone  and  that  of  a plane.  Near  the  source,  it 


can  be  shown  from  (3.4)  and  (b.5)  that 
■3/2 


gl  ~ C1S 


s -*■  0 


where 


— -1  -1/2  ? 

Cj_  = (1920Rc)  k tt  (1  - j ) (2  - 3 cos'  0)  . 

„d 


(4 . 6a) 


(4.6b) 


When  (4.4)  is  substituted  into  (3.3),  the  admittance  Y on  a cone  is 
decomposed  into  two  components,  namely, 

yd  = Y0  + Yi  • (i<7) 

The  singularity  of  g^  at  the  source  point  specified  in  (4.6)  is 
integrable.  Thus,  there  is  no  difficulty  in  evaluating  Y^  numerically. 

The  first  term  Y^  is  the  admittance  of  a slot  on  a plane.  It  is 
defined  by  (3.3)  after  replacing  g by  g^  in  (4.5).  From  (3.7)  and 
(4.5),  we  recognize  the  following  identity: 


(1  + -r  ^7>  G(s) 
k"  3y“ 


(4.8) 
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Note  Chat  (A. 3)  Is  in  the  form  of  (4.2).  Substituting;  (4.3)  into  (3.3) 
and  interchanging  integration  and  differentiation  operators,  we  obtain 

ff  itv  , . . 2 r r fly, 

Y1!  - - -=-  ! ' dv,  dz,  cos  ( — -)  {(1  +~— )[  ! dv  dz  cos  (— — ) G(s)]} 
0 ab  Jj  1 2 a k-  ^ JJ  I 1 3 


The  Integral  in  (4.^)  is  now  convergent,  and  (4.9)  can  be  considered 


(4.9) 


as  the  "finite  part"  of  the  divergent  integral  in  (3.3).  For  numerical 
evaluation,  (4.9)  is  converted  to  that  in  the  Fourier  transform  domain. 
Following  Rhodes  [10],  it  is  simplified  to  become  (Appendix  A) 


15fl\ 


r ~a  ) 

ImY  . * t~ 1 S 

0 , , s.  1 

15fl  k 


r3b 

C(ct)3  ( J()(t)dt  - Jx(Bb)  1 > 


(4.10a) 


da  C(a) 


+ — da  C(a)y 

TT 


">  2 1/2 

where  3 » (k-  - a“)  ~ 


■“  i-pb 

3 j YQ ( t ) dt  - Y1(Sb) 

- o 


K0(t)dt  + Kt(Yb)  - ^ 


(4.10b) 


3 0 1 'I 

Y =■  (a"  - k-)t,_ 


, and 


cos~  (aa/2) 
l - (aa/n)' 


(4.11) 


In  summary,  the  direct  contribution  Y on  a cone  as  given  in  (3.)) 
is  divergent,  due  to  an  "illegal"  Interchange  of  integration  and 
differentiation  operators  in  the  derivation  process.  The  (correct) 
finite  part  of  the  divergent  integral  Is  given  in  (.-*.7),  where  Y^  is 
given  in  (4.10),  and  Y ^ in  (3.3)  after  replacing  g by  g^.  The  same 
difficultv  arises  in  the  case  of  a cylinder,  and  it  is  treated  bv  the 
same  procedure  as  in  the  case  of  a cone. 


5.  NUMERICAL  RESULTS 

We  have  derived  the  input  admittance  Y of  a .slot  on  three  types 
of  surfaces:  (i)  For  an  infinite  plane,  the  final  solution  Y » Y ^ is 
given  In  (4.10).  (ii)  For  an  infinite  cylinder,  Y - Y^  Is  given  In  (4.7), 
where  Y^  Is  calculated  numerically  from  (3.3)  after  replacing  g by  g^. 

The  explicit  form  of  can  be  gathered  from  (4.4),  (3.32),  and  (4.5). 

(ill)  For  an  infinite  cone,  Y has  two  contributions  as  described  in  (3. 2'. 

. i l 

Y is  given  in  (3. IS)  and  Y in  (4.7).  To  calculate  Yj , we  use  (3.3) 
after  replacing  g by  g^,  where  g^  can  be  gathered  from  (4.s),  (3.4), 
and  (a. 5).  Numerical  results  of  Y on  the  above  three  surfaces  are 
presented  below. 

Slot  on  a plane.  As  a function  of  slot  length  a , we  plot  (a/2b)Y 
in  Figure  4 for  three  different  values  of  slot  width  b.  Those  curves 
are  practically  linear,  and  can  be  described  to  a good  accuracy  by, 
for  0.4  _<  (a/\)  O.b  , 

Y 7:  --{[  1.029  + JO. 59b]  + (3.75  + ]B)(-r  - 0.5)]}  millimho  (3.1) 

& ' \ 

where  l)  * 44,  33,  and  21  for  b • 0.000 1\,  0.001 \,  and  0.01X,  respectively. 

Relation  to  dipole  impedance.  As  discussed  in  [11),  there  is  an 
alternative  definition  for  the  (input  or  mutual)  admittance  of  a slot. 
Instead  of  (2.1),  a modal  voltage  V can  be  defined  by 

-*  ' • l 7 

E - z V ~ cos  l-y)  , (5.2) 


or  equivalently, 

_ ;b/2  , v 

V - ! (E  . *)v.Q  dz 

-b/2 


(5.3) 
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Then  a different  input  impedance  Y is  defined  by  (2.3)  after  replacing 
V by  V.  It  is  easily  shown  that 

Y - (a/2b)Y  . (3.4) 

Conventionally,  Y is  used  if  the  slot  is  fed  by  a waveguide,  while  Y 
is  used  if  the  slot  is  centrally  fed  by  a pair  of  transmission  lines. 

From  the  duality  relation  in  Maxwell's  equations,  it  can  be  shown  that 
(p.  519  of  [2]) 

2 - -~( 120n ) “ (2Y)  , (5.5) 

4 

where  2 is  the  input  impedance  of  a centrally  fed  dipole  radiating  in 
the  free  space  (not  in  a half-space  as  in  the  case  of  a waveguide-fed 
slot).  From  (.5.1),  (5.4)  and  (5.5),  we  find  chat  for  a half-wave 
length  dipole, 

2 « 73.12  + j 42.3b  ohm  , (5.n) 

which  agrees  with  the  results  in  [2], [10]. 

Slot  on  a cylinder.  Consider  a circumferential  slot  of  dimensions 
0.4"  x 0.4"  on  an  infinitely  long  cylinder  whose  radius  is  3.S".  Figure  5 
shows  Y calculated  by  the  presenc  OTD  solution  and  that  bv  the  modal 
series  solution  in  [5],  These  two  solutions  are  in  agreement  within 
0.3“  in  magnitude  and  one  degree  in  phase.  Note  that,  under  the  "one- 
mode approximation,”  the  modal  series  solution  [5]  is  exact . It  is 
amazing  that  the  present  CTO  solution  gives  such  an  accurate  result 
for  kS  - IS. 

Slot  on  a cone:  variation  with  radial  distance.  In  all  the  following 
cone  calculat ions,  the  slot  has  the  dimensions  of  0.5\  * 0.051,  except  where 
stated  otherwise.  In  Figures  b and  7,  the  slot  is  circumferentially 


lo 


| Y | in  millimho 


Phase  in  degree 


m rmlltmho 


jtar*  5.  Irpuc  aJaiiceance  or  a sloe  on  a cone  as  a function 


radial  distance  c. 


Phase  in  degree 


) 


oriented  on  a cone  with  0^  • 30°,  and  the  variation  of  Y with  the 
radial  distance  c is  presented.  We  observe  two  effects:  (i)  As  c 
is  increased,  the  radius  R = c sin  0^  of  the  "equivalent"  cylinder 
becomes  larger  and  larger.  The  magnitude  of  Y decreases  and  approaches 
the  asymptotic  value  of  the  slot  on  a plane.  (ii)  At  c = 2X,  the 
tip  contribution  |YC|  is  less  than  IZ  of  the  |y|,  and  this  contribution 
diminishes  as  c increases. 

Slot  on  a cone:  variations  with  cone  angle.  As  0q  is  increased, 
the  cone  surface  becomes  flatter.  Therefore,  Y in  Figure  8 approaches 
its  value  on  a plane. 

Slot  on  a cone:  variations  with  slot  length.  It  is  interesting 
to  observe  from  Figure  9 that  the  minimum  values  of  | Y ( for  both  cone 
and  plane  cases  occur  roughly  at  a - 0.45X,  not  at  the  resonant  length 
a * 0.5X. 

Slot  on  a cone:  variation  with  slot  orientation  angle.  Figure  10 
shows  that  there  is  about  a 10*  increase  in  |y|  as  uj  varies  from  0 
(radial  slot)  to  tt/2  (circumferential  slot'). 
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APPENDIX  A S IMPLICATION  OF  YZ 


The  input  admittance  Y of  a slot  on  an  infinite  plane  is  given 
in  (4.9).  In  this  appendix,  it  is  transformed  into  a form  more  suitable  for 
numerical  evaluation.  Making  use  of  the  identity 


ou  r . 

-k  ff  ..  ..  exp  H[yy2-yi>  + k2(vzi)]} 

30tr3  [I  y 2 (k2  - k2  - k2)1/2 

_oo  ' V Z 


(A-l) 


Equation  (4.9)  may  be  rewritten  as 


-a-k-  ff  a: 


60tt  k 


5,  ii  “y 


dk  dk  C(k  ) S(k  ) 


2 ? 

k - k 

v 

2 2 1/2 

(k“  - k“  - k )l/ 
v z 


(A-2) 


where 


S(kz) 


sin  (kb/2) 

(k  b/2) 2 
z 


(A-3) 


cos'  (k  a/2) 

C(ky)  = * 

1 - (k^a/TT) 


(A-4) 


This  result  is  identical  to  that  obtained  by  Borgiotte  [4],  which  is 
derived  from  a different  method.  Following  Rhodes  [10]  who  has  studied 
a similar  integral,  we  separate  (A-2)  into  real  and  imaginary  parts. 


Consider  first  the  real  part 
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2 2 1/2 

By  a change  of  variable  k = (k  - k ) cos  n.  the  inner  integral  of  the 

z y 

preceding  equation  is  transformed  into 


2 2 1/2 

r(k  - k 

I V7 


S(V 


v 2 ’l/1  : / / i 

(k'  - k‘  - kV7'  b'(k'  - k ) 1 


9 99  1/9 
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cos  n 
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(A-6) 


Differentiating  the  integral  on  the  right  twice  with  respect  to  the 


2 2 1/2 

parameter  b(k  - k ) , we  obtain 


i r 


T/2 


2 J 


9 ’ 1 / ■’ 
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, 2 2 l'\ 

which  is  equal  to  t/4  J,  (b(k  - k~)  Then,  Integrating  twice  and 

0 y 0 

applying  the  recurrence  formulas  of  Bessel's  functions,  the  integral 
in  (A-6)  becomes 


f fb(k--k-) 
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Substituting  this  result  into  Equation  A-5,  the  real  part  of  is 
reduced  to 
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Next,  we  study  the  imaginary  part  of  V*?,  which  is 
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Bv  another  change  of  variables,  * (k*  - “ cosh  n in  the  first 

. 2 2 1 / '1 

inner  Integral  of  the  above  equation  and  k “ fk  - k ) sinh  n in  the 

2 V 

second  Integral,  we  have 
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It  can  be  shown  that  the  integrals  *.  n the  right-hand  sides  of  C A— 11) 

and  (A-12)  are,  respectively,  related  to  Neumann's  functions 

(Y-,  Y, ) and  the  modified  Bessel's  functions  of  the  second  kind 
0 1 


(Kq,  K^)  by 

9 221/2 

sin"  (cosh  n(k'-k“)  ” b/2)  ? n l /? 

dn  7 - - 4 b(k--k-)  ; j 

0 cosh"  n l^o 


•>  ^ 1/2 

) 

v' 


/b(k“-k') 


dt  YQ(t) 


■■>  9 in 

- Y. (b(k  -k")  )-- 

^ nb(k~-k-)  | 

y J 


2 , 27172/  (A"13) 


and 


sin'  (sinh  n(k'-k")  “ b/2)  0 

dn  b(k"-k~)  ' ' 

sinh“  n 


u 


2 9 1/9 

cb(k  -k“) 1 

y 


dt  k0(t) 


9 ) 1/9 

+ K (b(k'-k')  ' ■)  - 
l v 


1 


1 9 1/9 

b(k'-k')  “ 

V 


(A- 14) 


. d . 

Substituting  these  results  into  (A- 10),  the  imaginary  part  ot  Y^  is 


reduced  to 


l Yg 

m 0 


a f fk 


2 ,2. 1/2  T fb(k"~kv 1 


1/2 


. dk  C(k  ) (k'-k' ) , 

i^li  y y v i 


dt  Y()(t) 


9 9 1/9 

- Y (b(k'-k') A/  ) - 


9 9 1/9 

Tb(k'-k')  ■ ~ 


I 

k 


9 2 1/2 

dk  C(k  ) (k'-k  ) ' 
n i v y v 


2 2 1/9 

/btk  -k  ) i ^ / 9 

dt  Kn(t)  + K (b  (k"-k  ) ' ) 
j 0 1 v 

0 


9 2 1/2 

b(k~-k  ) 


( A- 15) 


The  final  results  in  (A-9)  and  (A- 15)  are  duplicated  in  *.4.10)  after 


an  obvious  change  of  notations. 


E&WZZZ”**** 


APPENDIX  31 

COMPUTES  PROGRAM  OP  Y FOR  CYLINDER 


PROGRAM  MAIN  ( INPUT  > OUTPUT  » T APE3“QUTF'UT 
COMPLEX  Y 1 1 » YPL . YDIF • PY 1 1 , PYPL . PYD IP •€ J * CONA 

COMMON  A.B.R.AK.AKA.  AKS  .AKR.3B.SC.3E.THIRD.pl 'CONST r CJ ' CON A 
CJ=<0. >1.) 

PI=3. 1415926536 
THIRD=1 ,/3. 

READ*  .AIN.BIN.RIN.  PREC1 
YR-20.*Pt*FRECl/3. 

A=A IN*0 . 0254 

B=SIN*0.0234 

R = RIN*0  * 0254 

AKA“AK*A 

AKB=AN*B 

AKR=AK*R 

CALL  PLANE  C YPL . AN. A, B J 
CALL  DIFF  < YD  IF . 5 ) 

Y1 1 = YPL  ryD IF 
PYPL-YPL 
PYDIF=YD IF 
PY1 1*Y1 1 

CALL  XTOP  ( PYPL  ) 

CALL  XTOP  ( PYlt  ) 

CALL  XTOP  ( PYDIF  ) 

'JRITEI3. 1 J FREQ. A. AIN.B.BIN.R.RIN. yPL.PYPL. YDIF.pYDIF > Yll .PY1 1 
F0RMATC/T10. ‘INPUT  ! * • /T10. "FREQ  =‘.F3.3.‘  JHZ*. 

*/T10.‘A  * ‘ . F3 . 3 • " METER  = ‘ • F8 . 3 . ‘ INCH'. 

t/TlO.’B  * * . F3 . 3 » ‘ METER  =‘.F3.3.‘  INCH'. 

»/T10.*R  = * • F3 . 3 . " METER  *’.r3.3.‘  INCH1. 

T///T10. ‘OUTPUT  ! ‘ ./T10* ‘ Yll ‘ . 

ST23.  " REAL  ‘ .T33,  • IMAGINARY’  . T53.  ‘MAGNITUDE’  • To3.  ‘="MASE‘  . 

t/TlOf ‘PLANE* . T20.3E13. 7.E12.4. 

t/T 10* ‘DIFFERENCE’ . TOO  • 3E 1 3 . 7 . E 1 2 . 4 . 

S/T10. ‘CYLINDER* . TOO . 3E13 . 7 • £12 . 4 ) 

STOP 

END 


SUBROUTINE  SELF  ( RESULT  ' 

COMPLEX  FTH.FV.RESA.RE3D.RES1 .RE32.RE33 »RESY. RESULT. CJ. COMA 
EXTERNAL  FTH.FA.FY 

COMMON  A . B • R * AK  > AKA  * AKB  • AKR  * SB  . SC  . 3E  * TH  l RD  • F‘  I * CONS  T » C J * CON A 
COMMON  /DA T A2/  y l . 21 , YO . CL . CU 
Y1  -0  • 

21*0. 

38=0. ?*SC 
SD*0 . ' *SC 
7ML=P I/O. 

CALL  RGO  f FA  * 0 . * THL  . -3  * RES  ) 

RE3A*C0NA*SGRT<2.*A  ■ «R£3 


CALL  CGQ 

( ETH 

.0.  p 

7HL.-3*RE3D 

CALl.  cgc. 

< 1-  Y 

.0.  . 

Of 32  * k£3 1 

L A L «.  £30 

p y 

. *ru 

» .52  ♦ *‘E32 

CALu  CoQ 

•;  *=■  r 

* 3C » 

3 £ » !>  2 • a E 3 3 

<>;  **:;•£  *j  1 

•*S32 

♦K  £3 

3 

- i l L’L  7 = *• . 

« - *£ 

E 3 0 ♦ R E 3 Y 

' ^ 

PRECEDING  PeGB  hi-AilC 
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SUBROUTINE  DIFF  i COIF. NY  ) 

COMPLEX  Fy; , c J . CONA , SPLINE » YD  IF 
COMPLEX  TOTAL/RESULT. RES1 > RES2 • RES3 t RES4 
REAL  SPtlC) 

EXTERNAL  FYC 

COMMON  A.&.R.AK.  ANA  r ANB  .AKR»SB*SC»3E»  TH  I RD  * P I r CONST  » CJ  » CON A 
COMMON  /DATA2/  Y 1 / Z 1 • YC  > CL  > CU 

DATA  5P/7.*lA.fl4.»lo.>14.*13.fi4.»l3..14.ilA.*l4.»16./ 
CONA«(i.~CJ>/<  1920.#R*PI*SaRT( AK*PI  > ) 

AN4>4*NY 
N2»C*NY- l 
DY=AYAN4 
0Z»B/  10. 

3C=AMIN1 (DY/3. > DZ  < .IE-1) 

3E=3.«SC 

CALL  SELF  RESULT  ' 
f DIF  = < 0 . > 0. ) 

DO  COO  IZ“l • 3 > 2 
IFIIZ.EQ.l)  EZ«t. 

IFCir.NE.l)  £2=0. 

Z1=0Z*< IC-1 > 

SF'LINE  = ( 0 . -0.  ) 

DO  100  IY=l. NC,0 
I3P’(IYK)/2 
Yl=DY*( I Y - 1 > 

RES1*RES2=RE33=RE34»(0. >0. ) 

ZL--B/C . 

ZU=  B/2 . 

YL=-A/2. 

YU=Yl-SE 

tF(YU.GT.YL)  CALL  CGQ  ( FYC . YL . rIJ , <3 , RES  1 ' 

YL=Y1+SE 

YU=A/C. 

IFIYU.GT.TL)  CALL  CGQ  < FY2fTL*YU><?»R'£S2  > 

YL*Y1-S£ 

YU=  Y 1 
ZL»-8/2. 

ZU*21-SC 

IFCU.GT.ZL)  CALL  CGO  « FYC  * TL  . YU  • *3 » RES3  > 

ZL*Z1PSC 

ZU-8/2. 

IFCZU.GT. ZL)  CALL  CGQ  < FY2 > YL • 'U . Sf RES4  ) 

TQTAL=RESUL  T TRESl+RESCfRESSTRESA  *-R£33  TR  £54 

total»total«cos<p:*yi,  a; 

SPLINE -SPLINE  + SP  < ISP '« TOTAL 
CONTINUE 

spl:ne=2.*by*spl:ne.'13. 

YDIF»YDIF+EZ«SPLINE 

CONTINUE 

YDIF=DZ  HYDIF 

YDIP*-9. RYDIF/A/B 

RETURN 

END 


• TiJ . S » RES  1 


t— -4_J 


IHIS  PA®  IS  BEST  QUALITY  FRACTICABI4 
fftOJI  OOPY  7UFBISHHD  10  W>Q  -o — ' 


F JNCTICN  -A<TH>  * 

COMPLEX  CJ.CONA 

COMMON  A * E<  * R * A N .AKA.AK8.ANR.S8.  3 C » S E . ~H  I R 0 • F*  I .CONST  .CJ.CONA 
COSX-COS<  TH> 

COSX  *A8S  V COSX ■ 

:■  -PKCGSXFSP.- A 

CALL  FRESNEL  < C > 3 * X ) 

F A»C*COSX«i E. -3. *COSX  FCOSX ) *SQRT . COSX ) 

a E TORN 

END 


COMPLEX  FUNCTION  FTH  ( TH  > 

COMPLEX  P'S. Cj.CONA. CONY 
EXTERNAL  F3 

COMMON  .DATA/  S INC > COSO . S IN4 , C0S4 . COSOX . XCCN . YCON . CONY 
COMMON  A.D.R.AK .ANA.AK8.AKR.SP.SC.se. TH IRO . P I . CONST . C J • CONA 

COSX  *COS  < TH  ) 
caso*cosx*cosx 
3 INC  * l . -COSO 
3 tH4»S IN2«$ tNO 
;2S4=00S2*C0S2 
:052X«  l . - J . «S  tNO 
.1  C)N»AKR*C0S4  o. 

CON*XCON**TH IRD.  R 
t CUN*F  1 «C  JSx.  A 

cONsr.oos-i.  . o.*akr«akr’ 

: CNY*CJ*CGNST««THIRD 

Call  '000  ( : 3 . 0 . . 58 . 30 . FTh  ' 

RET  UK  N 
END 


complex  function  fs<s> 

COMPLEX  AX.LX.YPX.OPX.OS 

'"OMPLEx  CCYL  .GPL.  QD  IF*  APP.CJ.  JOKS.  CONA.  CONY 
COMMON  'j  A'A.  SING » COSO  .SIN4.C0S4.C030X.XC0N.  'CON.  CONY 
COMMON  1 • y • K . AK .aka. AK8. AKR » SB . SC . SE*  THIRD • P I • CONST . CJ.CONA 
'3  = 0.  -0.  ’ 

lt.  ;.e-io>  return 

AK j - AK  *S 
JOKS  *C  J >K3 
<*XCON«S 

CALL  cOCK  x. jx.ux.ypx.upx  > 

: ;»cmplx  cos i aks ' • -s  r,x(  aks ) ) aks 
CCrL*'. i»(SINOtJOKS»COSOX) 

> -OX «jCK3*< COSOt ( 1.-2. *J0KSH--0KS«SIN2> 

. »YPX*CCNY*Si'NO 

t »gpx  «CONY*wOKS»COSO 

;f  ; jso  . ce  . t . .* -so  ocyl«scvl*.jpx*conyrs:n4  ;csg 

VCSO.L*.  t .£-50)  JCXL  *OCYL  *0 . 3~3*  < 1 . -CJ  ' «SOR  T -I«Ak3'  ANR 

:ipl  -s:n:p_uks*<o. -3.  «sino' tt  i .- joks' 

3D  IF  *3* ( OCYL-OPL ) 

IF  3. or.  1 .£-!)>  30TC  10 
-.E*REAL.GDIF) 

'•DIF  »C  MF  . RE  > - RS  > 

;;i:e-"is*oi'lF 

\P,:  --  '0na«C03C«  . 0 . -3  . *COSO  ' SJF  ’ S' 

-•i*>3o;p“Apc 
F3«‘  5 * C OS  > v CON  ) 

■ E'ajRN 
1(1 


1 


,vri 


,cU^ 


couple < ;|  notion  ryiyj1 

COMPLEX  -C.CJ.CONA 

external  fz 

COMMON  A . B . R > AN  . ANA.  ANB  . ANR.5B.SC.3E.THIRD.PI  . CONST  . Cj -CONA 
iummcn  data’.-  Tt.2i.rr.2L.JU 

TYs/C 

IF; SB .GE. Y2>  ZLsSC1RT ( 3B*SB-Y2*Y2 ' 

IF  < 3B.LT . Y2  ' ZL  = 0. 

CALL  C0Q2  < FZ.JL-SC.S.FY  ) 

FYaFY«COS(PI«Y2/A 

RETURN 

END 


COMPLEX  FUNCTION  F2<22> 

COMPLEX  UX.UX.UFX.UPX.LiS 

COMPLEX  i3CYL.UPL.CJ.  jQN5.C0NA.C0NU 

COMMON  A . B • ft . AN . ANA . ANB . ANR > SB • SC  > SE . THIRD .P I « CONS  T . C J . CONA 
COMMON  DATA  J/  T l . 21  • TO  . JL  . JU 
DY»Yt-Y 1 

d:*:j-zi 

5 = SllR T DY«DY-0J*D2» 

:f  s.lt. t .s-io'  return 

C03X»DY  3 
COSJ»COSX*COSX 
U I NC - 1 . -COSC 
RINATS INCxS INC 
CUE  A'-COSQxCOSC 
i '00  OX  * l . -C.  «S  INC 

»(\  j - An  y‘3 

jON3*CJ  anS 
, anr«C0SA/2 . 

=><*<THIRD*S/R 

ALL  FOCN  ■ X.UX.1JX.UPX.UPX  > 
iS-*CMPLX<  COS  CANS)  .-3IN(ANS>  > ANi 
,:ONsr*cos.i/  cc.  «anr«anr> 

,'OMU -Cj*CC)NST  *x Th  IRD 
3CYL*UX*< SINC+ JONSpCOSJX > 

» .UX*jON3*l  COSCtL 1 . -C . «JCNS 1 **CNS  «SI  -C  ' 

s v.'px*ccnu«s:nc 

i .UPXXCONU* jQn3*C0S2 

IF  COSO.C-E.  l.E-IO'  0CYL*'3CYL*‘UPXYC0NVXSIN4  COSC 
I F ( COSC . C T . ! .z-ZO  ' GC YL *GCYL FO . 375*  i : . -Cj  ' «Si)RT>  F I HANS'  AnR 
3PL *3IN2+ JQN3#< C . -3 . «S INC > X ( t . -JONS ' 

"0  *03  * jCYI.-CSPL  ' 

RETURN 

•:nd 


COMPLEX  FUNCTION  FY2l  VC ) 

COMPLEX  FJ.CJ.CONA  . 

EXTERNAL  - C 

COMMON  A. B.R .AN .ANA. AND .AnR . SB "SC • SE • THIRD  <P I ■ CONST . Cw .CONA 

common  data:.  rt.ct. ty.cl.cu 

TY  * Y 0 

DY  * ABS  ( Y2-  T I » 

Nt*A 

:F,D'..Lr.  . IE  - 1 NI*la 
IF  DY.lT. . IE-1 ' ni*j: 

CAL*  COQJ  ffJ.CL. JU.NI  .F>2 
- ■ C = FYC*COS . I «YC  A i 
\ E ** 1 R N 
END 


3 


i 

ji 


' 


APPENDIX  3i 

COMPUTER  PROGRAM  OF  Y FOR  CONE 


program  main  , input. output. tape3»output  > 

COMPLEX  (PL  • YDIF • YTIP . YD IR . Y V l . PYPL  > PY  D IF • PY  T IP » Py9 IR • PY  1 1 

COMMON  P I . AN  > A.  B • BE » S3  > Ba  > C < JMEO  • 3 INO  > COSO  > 'AnO  • DEL 

PI  »3. 1 4l5«0iS3a 

REAP  *.pREQ.A.B.C.OME» THETA 

an*CO.  *P I HFRECl/  3. 

TH£TAO»  THET  A*P  t 130. 

OMEO*OM€*PI  130. 

SINO-SINlOMEO' 

CQSO»C03tQME0> 

T ANO" ( AN ( '4ETA0' 
fC«S/C . 

S3  ’9  5. 

>. 

CALL  PLANE  , YPL.AN.A.B  ' 

CALL  DIFFER  v rPIF • 3 ' 

r ANO"  THET  AO 

CALL  TIP  ( (TIP  ' 

AIN-A/0.0034 
S CN"6- 0.0E34 
CtN»C  O.0CI4 
Y D IR*YPL  ► YDIF 
•ii*yoirpytip 
PyPL  *YPL 
r (OIF *YD [F 
PYDIR-YOIR 
,-YTIP.YTIP 
F'H*YU 

CAL_  XTOF  v PYPL  ' 

CALL  XTOF  ( PYPIF  ' 

CALL  XTOF  v PYPIR  ' 

CALL  XTOF  v PY  T IF  > 

CALL  XTOP  . PyU  > 

'JR  I TE%  3«  I "REG  • THETA.OME.A.  AIN  • B . 9 IN . C .0  IN . 

«YPL .PYPL»  YHIF. PYPIF. YPIR.PYOIR.YTIP.PYTIP. 'll .PY’. 1 
S'GP 


FORMAT  4 TlO.‘ INPUT  1 

•.  TIO.’FRElJ 

• •.!=■  u • jh: • . 

’”.0. 

* THE  T A »• .F • E. ' 

. • DEGREE  •• 

r io. 

•OMEGA  * ' • F 1 0 . ’ 

.•  DEGREE  ’ • 

TV  J. 

*A  •••PIC.’ 

.•  mE'ER  *• 

inch*  • 

'10. 

• 9 •••PIC.’ 

.•  mE'ER  «* 

• r 12.  "•  * INCH*  * 

' 10. 

•C  •••PIC.' 

.*  METER 

• F12  • ■**  • l^CH*  • 

' t 

■ >. -output  r 

10. * Y 1 1 * . 

i 

* 3 * * 

PEAL1 • '38* • IMAGINARY* . 13  • • 

‘lAlSNlTUDtE*  • ,C,HASE 

10. 

•PLANE* • TC0.3EI 

'10. 

•difference* • TC 

0 . 3E II . '.tic 

. 4 • 

'to. 

•I  IRECT  • . TC0.3E  13.  '.£1C.  •*. 

no. 

•TIP* . 'CO. 3EiI. 

'•eic. a. 

no. 

• :gne' . too. 3eis 

.'■•e:c.-»- 

END 
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_ IS  best  9UU.I"  re*cnctB« 
^Fo^^^ru“C  ' 


SUBROUTINE  DIFFER  < YDIF.NY  > 

COMPLEX  YDIF.C J.  51, ’MY  • RESUL  T 

COMMON  PI  > AN  » A • B • 80  > B3 . BA  « C > GMEO • 3 INO • COSO  » T ANO  » DEL 

da  r a cj/ ( o . . i . > . 

NO»C*NY- l 

dy  ^a/ float  . o*ny  > 

Y 0 1 F « ( 0 . -0. ' 

DO  100  IZ«  1.3 
Z2»83*i IZ-O' 

SUMY»*0. .0. > 

DO  OOO  IY«1 .NO 
Y1«0Y*< IY-NY) 

CALL  HYO  ( RESULT. Yl.Zl  ) 

RESUL T»RESULT«CCS*P  C * Y 1 / A ) 

rs»<  i iy+i  > o)-'  rv.  o> 

SUMY»3UMY+RE5ULT*<  14+15*0) 

OOO  CONTINUE 

YDIF»rDIF+SUMY*OY  15. 

I 00  CONTINUE 

YDIF=»CJ*AN*AN*B3«YDIF  f 1 00 . *P [ *P IXAxS ) 

RETURN 

END 


i 


SUBROUTINE  TIP  1 YTIP  > 
complex  cj • yttp. r .sigma 

COMMON  PI.AK,A.B.B0.B3.3A.C.0ME0.SIN0.C0S0. theta. [ El 
CJ-<0. . 1 . ) 

3A«1 .3037/THETA-l.  .■’33  + 0.  '7"0*THE'A- 1 . 4J9* ThET A* T JETA 
3B*0  . ’193  + 1 . 4a08*  THETA-  l . 1 0«3*  THE  T A*  T -l£”  A +0  . ,5*o*  rMETA»*3 
SIGMA=SA*<  CCS< SB > +CJ*SIN( SB > 1 

ANBCaANXB/ 0 . 

5INX*SIN(  AKBO/ANBO 
AKCO*C. *AN*C 
r »COS ( ANCO ' -CJ*SIN( AKCO) 

T»3IGMA*A*8*SINX*SINX*  1 l . +C„ ' *T  „ 30  . «F  t «X4*C*C  ' 

T + T /SQR T < S IN  < tm£T  A ' *COS  ( THET  A t4.  «P  r > 

YTIP«T*SIN0*SIN0 

RETURN 

END 


IOMPLEX  FUNCTION  “A.TH1 
COMPLEX  Cj 

COMMON  P I . AN , A . B. 90 . 63 . 9A  »C  . OMEO . 3 I NO • COSO . * ANO . CEL 
COMMON  /DA TA/  Y1 .01 ,Rl .pNl.rO* TTH 
DATA  CJ/10. . 1 . )/ 

C0ST«AB3'C0S<  TH' ) 

ARG-PI«DEL«COST  A 

CALL  FRESNEL  ( CF.SF.AfiG  . 

CON-SORT.  0.  *A,  COST  ' 

-A«S0RT<  DEL ) «SIN( AfiG  > -CON  0. «SF 
F4«-CJ*AN *DEL/ARG«FA+C0n«CF 
5INA»SIN(0ME0  + c'Nl-TN) 
r»*.4li0.-3.  «CCST«CCS'r'  «S  INA*S  INA 
RETURN 
END 
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SUBROUTINE  HY2  i RESULT,  n.zi  I 
COMPLEX  CJ»Fv2.FTH.FA.RESl .RES2.RESULT 
EXTERNAL  FY2.FTH.KA 

COMMON  PI .AR.A.B.B2.B3.BO.C.0ME0.3IN0.C0S0. 'ANO.OEL 
COMMON  /DATA/  TY t . TZ I . R1 . PHI . Y2 • TH 
DATA  ZERO.CJ/ . IE-50. (0. • 1 . ) / 

TY 1 * Y t 
TZ1*C1 

DEu'O. 95«AMIN1 (( 92-Z1 ) . < Z1 TB2 ) > 

SI SQ*Z 1 tZl +Y 1 *Y l 

IF ( S l SCI . L T , ZERO  > SOTO  1 00 

3 1 *3QR  T i 3 * SO  ' 

3IN1-Z1/S1 

C0S1*Y1/S1 

3IN0«3IN0*C0S1-C0S0*SIN1 
cosd*coso«cos:*si.vo«s:ni 
RlSQ»C*C+3l 3Q-2  . *0*31  *COSD 
R1»SQRT(R13Q) 

PH  1 JAS I N ( SI HSIND/Rl > 

SOTO  200 
ICO  R1»C 

DH1*0. 

200  CONTINUE 

THL-OMEOtPHt 
THU*  THLtP  t,'2. 

CALL  CGQ  F A . THL • THU . 3 • RES  1 > 

THL*  THU 
THU*  THL  YP I / 2 . 

CALL  CGCI  i FA. ’HL.  THU. S.RES2  > 

RESULT* l t . tC  J > *SGRT  <■  P I AR ’ ' ( A . HAN dANtRl «TANO > « i RE31 *RES2 ' 

THL»OMEOtPHl 

THU*THL+PI 

CAL.  COO  ( FTH. THL.THU.3.RE31  ' 

THL*  THU 
THU* THLtP I 

CAL.  CGCI  ( FTH.  THL.  THU.S.RES2  > 

RESULT  *RESUL ” *RE3 1 tRES2 
RES  I * . 0 . -0. ' 

YL*-A/  2. 

YU»Y1-DEL 
Nl*Io. *< YU-YL  > / A 

tF< YU.OT.YL)  CALL  CGQ  . FY2 > YL  > VU • N I > RES l ' 

RE3UL’’*RE3ULTtR£Sl 

NI*12.«DEL  B2 

■'L-.YU 
VIJ*T  1 

CALL  CGQ  ( FY2.YL.YIJ.NI.RESl  > 

YL*YU 
YU»Y 1 +DEL 

CAL.  CGQ  C FY2.YL.YU.NI.RE32  ’ 

RE3UL  T*RESULT+RE  3 1 tR£52 
RE52  * l 0 • .0. ) 

VL»TU 
YU-A.  2. 

Nl*la.  t(y'J-'t)  / A 

IF  Y'J.OT.TL'  CALL  CGQ  FY2 . Y|. . YU  • N I . RE32  ' 

RESULT "RESULT *RE32 

RETURN 

END 
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complex  function  fth.tn) 
complex  -5 

EXTERNAL  ps 


COMMON  PI  . if. . A.  B.  BO.  S3. SA.C.GME 0 • 3 INC.  COSO  • r ANO  • DEL 

COMMON  DATA/  Y 1 . C*.  . R l . t HI . 1 0 . TTH 

TTM*TH 

CALL  CGGO  ( F3 . 0 . » DEL  * l A » FTH  > 

RETURN 

END 


COMPLEX  FUNCTION  -SiS> 

COMPLEX  F 20 

COMMON  PI . AK .A. R.  SO. S3. Ba . C . OMEO . a INO . COSO • TANO-DEL 
COMMON  'DATA/  Yt.2I.RI.PHl.YO.TH 
YO»S*COS<  TH) 

00»S*SIN( TH> 

PS-*F00vC0) 

FS»3*FS«C0S(PI*Y0.  A) 

RETURN 

END 


COMPLEX  FUNCTION  FYO<YO' 

COMPLEX  FCo.AESl •RESC.RE33 
EXTERNAL  FOO 

COMMON  eI . AK.  A.  8.  BO.  63.  BA.C.  OMEO  • i I NO.  COSO.  TANO  > DEL 

COMMON  . DATA/  Y1 . 21 . ft 1 .PHI . 1*0 , rH 

TYO-YO 

RES1»RES0»RES3«'.0.  . 0.  > 

DY»ABSCY0-Y1 ) 

N 1 *N0*  3 

IFIZl.OT.  DEL)  N0*la 
IF( 21 .LT. -DEL  1 Nl»la 
IF' DY.oT.DEL’  00 TO  1000 
20  «SOR  T ( DEL*DEL-DY*DY ) 

OU-  BO 
OL’OIPOC 

IF . OU. OT . 2L ' CALL  CGQO  < FOO . OL • ZU « N1 . RES1  ’ 

0U-01-2C 

0L»-SO 

IF ' OU . OT . 2L ’ CALL  0300  i FOO * OL » OL » NO » PESO  ' 
FY0-RES1 PRESS 
GOTO  5000 

1000  IF' Dr.GE.S3'  GOTO  OOOO 
OU«RO 

0L«2I+C£L  0. 

CALL  COUO  . F00-2L.2U.N1.RES1  ' 

2U*2L 

2L*21-DEL/0. 

CALL  CGOO  ; FOO • OL . OU • 1 a . RE50  ' 

0U*2L 

2L*-B2 

CALl  CGQO  . FZO.OL.OU.NO.RE53 
FYO*REo 1 yRSSOtRE  33 
GOTO  3000 
2000  2U-S2 

OL*-BO 

CALL  CG HO  ( FOO.OL.O’J.  J-Fy;  . 

3000  Fy;sPv2hcoS<P[«yO  A' 

RETURN 

END 
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COMPLEX  FUNCTION  FZ2(Z2) 

COMPLEX  U.U.UP.UP.GS 

COMPLEX  C J . J0K3  * HB  • HT  » GC  * GP  * GA 

COMMON  P I . AK » A . 8 . 82 . B3 . 86 . C » OMEO . 3 t NO  . COSO  . T ANO  • DEL 
COMMON  /DATA/  Y 1 # Z 1 < R 1 » PH 1 . Y2 • TH 

DATA  ZERO. THIRD .CJ/. IE-30. . 3333333333333.(0.  . 1.  > / 

32SQ-Z2FZ2+Y2FY2 

IF (S3SQ.LT. ZERO)  GOTO  1000 

S2»SQRT<S2S0) 

StN2*Z2/32 

C0S2»Y2/S2 

3INO»SINO*C3S2-COSO*SIN2 
COSD=*COSO*COS2YSINO«SIN2 
R2SQ*C*C+S2SQ-2 . *C*S2*C0SD 
R2»SQRT(R2SQ) 

PH2»ASIN(S2*SIND/R2) 

GOTO  2000 
1000  R2=C 

PH2=0. 


2000  ANG*ABSv PH2-PH1 ) 

DZ=Z2-Z1 

0Y»Y2-Y1 

SSQ*DZ*OZ+DY  *0Y 

IF i S3Q . L T . ZERO ) GOTO  3000 

S«S0RT(S3Q) 

IF < ANG . LT . ZERO ) GOTO  LOO 
SINA-SINC ANG) 

SIN01 =R2*S INA/3 
3IN02*RHSINA/S 

RT»SQRT(R1*R2)*TAN0/SIN0:/SIN02 

RTQRB=SQRT( C 1 . -SINOKSINOl ) *( 1 . -SIN02*SIN02 ) ) 'SIN01/SINC2 
CONP=  < 2 . *AK*AK*RT*RT > ** ( -THIRD ) 

tausq»sina/ang 

TaU-SQRTCTAUSO) 

TAU3»TAU*TAUS0 

ZETA3=AK*R1 *SING1*ANG**3/ < 2 . *TAN0#TAN0 > 

ZETA^ZET  A3** THIRD 
GOTO  200 

100  3IN01=SIN02=Z£TA»CONP=0. 

TAU=TAU3=1 . 

R8-TAN0*SQRT<Rl*R2> 

200  CONTINUE 

CALL  FCCK  < ZETA.O.U.OP.UP  ) 

AKS»AK*3 

JOKS=CJ/AKS 

H8=< 1 . -jOKSXTAUtO+jONS^JOKSFTAUGFU+CJ^CONpurAUFOP 

IF < ANG . GE . ZERO ) HB=HB*C-*CONP*RTORB* TAU3 PUP 

IF< ANG. LT. ZERO)  HB=HB+ . 375* ( 1 . -CJ ) *SQRT ■; P I *3 > R8 

HT=TAU*Y+ ( 1 . -2.*J0KS)«TAU3*U+CJ#CGNF*TAU3*UP 

HT=JCNS*HT 

3INTSQ-DZ*DZ/SSQ 

C0ST5Q=1 . -SINT3Q 

GC=HB*S INTSQ+H  T PCOSTSQ 

0PaSINT3G+ JONS* ( 1 . -JONS ) * ( 2 . -3 . FSINTSQ  > 

GS=COS ( AK3 ) -CJ*SIN  < AN3 ) 

G3=GS/AKS 
FZ2=GS*  < GC-GP ) 

IF( 3 . GT . DEL  ) RETURN 
3A*  < 0 . «0. > 

IF' ANG. GE. ZERO) 

t GA= ( 1 . -CJ*AK3 >*< I .+CJ)*<3.*SIVTSQ-1 . > «SQRT < P I / AK3 ) 

* *5IN01 H3IN01/ (3.PANS*AK*ftl kT  ANO i 

FZ2=FZ2-GA 
RETURN 

3000  FZ2*GA*<0. .0. ) 

RETURN 

END 
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APPEND  I X 33 


SUBROUTINES  OF  Y 


SUBROUTINE  PLANE  ( YPL  > I AK  > TA  » TB  ) 
COMPLEX  YPL.CFt .RES1 .RE30 
EXTERNAL  CF1»RFI 

COMMON  PLANAR.  ='t.AK.A.B,AKA.AKB.CQNC 

OAT  A PI. '3.  lAiJ9Ca53«S/ 

ak*tak 

a»:a 

AK  A*  AN  Xt  A 
»K&*AK  KB 

conc*r  : *x3*b*o  . s,  a**4 

=ERIOB»C. *P  C/A 

R£3 1 = XL  =0 . 

xU-PERIOO/3. 

if  xu.ge.ak>  goto  coo 

call  osar  < cfi  » xl  » xu • s • resi  ) 

CIO  <L-XU 

<U*XL  -PERIOD 

IF.XU.GE.AK  GOTO  COO 

call  CG03  ; CFI  • XL.  XU.IJ.RESO  > 

R£Sl »RE3i RRESC 
G0T0  CIO 

COO  CALL  CGQC  ( CFI i XL t AK * 32 > RESC 
R£3i *RES1 +RES2 

call  sic:  c si>c:.aka  > 

R£33*C0S l AkA/C • ) 

RS33-C . «RES3*RE53-AKA*SIN< AKA  -AKA  XAKAPC 
RE33*RE33*C0NC/  < •> . *Ak*AK  ) 

NU*C . XAK. PERIOD 
x'J*;  NU+C>  FPERIOD/  2 . 

CF . XU.LS.AK>  XU=XU+PSRIOD/C. 

CALL  RGQC  ( RFI  - AK  , XU  > e) « RES4 
R£33*RES3+RESA 
IOO  XL*  XU 

<U*XL+RER:on 

CAL-  RGC1C  RFI  • XL  > XU • 3 f RES  J ' 
:;E33*RE33+RE3A 
VATIO*AeiS<R£SA/ RES3  > 

:f<rat:o.3t. i.e-5)  goto  no 
tpl*res:+cmplx(o. <-c . »res3/p:  ' 

'PL  * TPL  «A/  l IS.  HARBOR'  I PX4  / 

RE’URN 


AJ^r 


THIS  PAGE  IS  BBS!  QUALITY  PRACTICABLE 
PBOJI  OCkPY  PUKAISrlfcU  TO  LLC 


COMPLEX  FUNCTION  CFIARY) 

COMPLEX  CFO 

COMMON  /PLANAR/  PI .AK/A»S. ARA. AND.C0N2 
X*B*SORT( Ar «An-ARY*AKY ' 

CF  l =CY  AK  Y > *CF2  ( X > 

PE TURN 
END 


FUNCTION  RFUAIxY) 

COMMON  /PLANAR/  PI / AK , A > 9 > AKA » AKD / C0N2 
X=e*SQRT ( AKY*AKY-AK*AK ) 

P1»CY(ANY) «RF2(X> 

P2=C0S< AKY*A/2. > 

P2»C0NZ*P2*P2/ AKY* *3 

rfi*p:-p2 

RETURN 

END 


FUNCTION  CY'.AKY) 

COMMON  /PLANAR/  PI t Ai\»A»&r AK A , AK D » CON2 

CY=C03< ARYPA/2. ) / ( L . -<  ARY* A/ P I > «*2 ) 

CY=CY*CY 

RETURN 

END 


COMPLEX  FUNCTION  CF2(X> 

CALL  CHOI  ( ZJOI  • Z'fOl  • X > 

REAL*  <* ( ZjOI-CJI(X>  ) 

AIMAG=<*(  ZYOI-ZYl(X)  > -2 . / 3 . 1 4 l3i>2o536 
CF2*CMPLX<  REAL > -A  I MAG  > 

RETURN 

END 


FUNCTION  RF2.X) 

RF2*X*(  ZNOKX)+ZKl(X)  )-l. 

RETURN 

END 
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APPENDIX  C 


SUBROUTINES  FOR  SPECIAL  FUNCTIONS 


SPECIAL  FUNCTION 


SUBPROGRAM 


Bessel's  Functions  Jg(x),  J^(x) 

ZJO(X),  ZJ1(X) 

Neumann's  Functions  Y^(x)  , Y^(x) 

ZYO(X) , ZY1(X) 

Modified  Bessel's  Functions  Kq(x),  K^(x) 

ZKO(X),  ZK1(X) 

[X  rX 

j JQ(t)dt,  YQ(t)dt 

CALL  ZHOI  (ZJOI,  ZYOI , X) 

KQ ( t ) d t 


ZKOI(X) 


Fresnel  Sine  and  Cosine  Integrals 


CALL  FRESNEL  (C,  S,  X) 


1 fX  sin  t „ 1 fX  cos  t 

S (x)  = dt  , C(x)  = dt 


/2tt  £ /t 


/2?  £ /t 


Sine  and  Cosine  Integrals 


CALL  SICI  (SI,  Cl,  X) 


Si(x)  = dt  , Ci(x) 


Fock's  Functions  V(x),  U(x) , V' (x) , U' (x) 


CALL  FOCK  (X,  V,  U,  VP,  UP) 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
FR)UM  OQFY  JTUWUSHJSU  TO  COO  — - 


"UNCTION  3-iO.X) 

if. x. 37.3. > goto  1000 
«.3*X. 3 . 

;_o» ; . -3 . 349699 -*x3**3 

1 -l.3ala£08*x3**4 

3 -0.31a3Sao*X3**a 

' +0.0444479*X3**3 

4 -0.0039444*X3**10 

5 rO • 3003100<X3<<13 
RETURN 

1000  X 3 * 3 . / X 

F0»  0.79~3343a 
1-0.0000007"*X3 
:'0.00533740*X3*«3 
3-0.00009513*X3**3 
4+O.0013?337«X3**4 
■i-0 . JOO  "3'303-«X3**3 
3+0.0001 44Ta*X3**o 
+HET4-X-0. "353981a 
I -0 . 04 1.ja39"*x3 

3- o.oooo3954*x3«*3 
3+0.  003a3573*X3**3 

4- 0. 000541 35*X3**4 
3-0. 0003 5 333 *X 34*3 
■> +0.0001 3553*X3**a 

I *C3S  . TNE"A ) 3GR7.X> 

5ET URN 

END 


FUNCTION  3N0iX> 

IF ( X . lT . 1 . £-30 ) GOTO  3000 
IF  1 X . GT . 3 . ) GOTO  1000 
T»X,'3. 

:no»-4logo  .3*x>  «zxov  :o  -0.37-31:00 
140.433"3430*r«*3 
3+0.'330o9"5a*T«*4 
3+0.03488S90*T«*a 
4+0 .003a3a99*T**8 
3+0 . 000 10750*T**10 
6+0.00000740*7**13 
RETURN 
1000  T-O./X 

OK 0*1. 35331414-0 . 0 7933333 XT 
t+0.031395a3*T«*3 

3- 0.0l0a344a*T«*3 
3+0.0033"373*"**4 

4- 0.00351540*7  «*3 
3+0. 00033308  *T*.xo 

ZK0»3K0«EXPl-Xl  -’SGR  T < X > 

RETURN 

3000  3N0*  1 . £30 

RETURN 
END 


function  :uiix> 

IF. <.GT»3.)  GOTO  lOOO 
<"*  x . 3 . 

3 . 1 *0.3-0 . 3o349983«X3«*3 
•0 .31093573**3**4 

3 - 0 . 0393 4339* x 3**o 
5+0 . 004433 1 ’ «X3**8 

4- 0 . 00031 7ai «x3«* 10 
5+0 . 00001 109*x3«* 13 

3..  1»3-1*X 
RETURN 
50  X3»3 . ■ X 

- 1 *0 . " 9"3343a+0 . OOOOO 1 3o*x3 

: -O,oia39ao7*x3**3 
3 + 0.  >001  "105«X3**3 
3-0 . 003495 1 1 «x3«*4 
4+0 . >011 "o5J*x3*  «3 
3-0 . 0OC30O33*X3«*6 
THET4»X-3.33o19449 
1+0. 13449a;3iix3 
3+0 .00003850  *X3  **3 
3-0 . 0063-37 9«X3*  *3 
4+0 . 00074348«X3**4 
I+O.  >00"<’934«-.(3*«; 

3-0. >0039 1 aa*X 3**a 
3-1*+'  1 «CCS.  THET4)  SORT i <> 
\ETURN 
END 


FijNCT  ION  3N . . < > 

IF.X.LT. 1 .£-50  GOTO  3000 
IF  1 X . GT  . 3 . > GOTO  1000 
T*X/3. 

3N1»X*8UGG v O.S*X> *31 1 . <>  *1  . 
1+0. 154431 44*T**3 
3-0.a7373579-*T x « 4 

3- 0 . 1 3 15689" *7 x*o 

4- 0.01919403*7**8 
3-0.001 10404*7**10 
3-0 . 0000 4o 3a  *’'**13 

3M»3N1  X 
RETURN 
l 000  7*3.  \ 

ON  1 * 1 . 3533 1 4 i 4 +■> . 33  4'’Sa  1 9 « + 
1 -0 . 03835o30*  ' « *3 
3+0.015043a3*T**3 
3-0.007303S3*T««4 
4+0.00335al4«T«*5 

5- 0.0 00a3343*T**a 

ON  1 » 3N 1 *EXP i -X  • SORT . X ) 
RETURN 

3000  OM  * 1 . £30 
RETURN 
END 
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function  ;yo<x> 

IF< X.lT. I .£-30)  joto  2000 
tFiX.OT.3.)  30T0  1000 
x3-X.  3. 

;yo»2..  3. i4i3^:a34»^LJOi0.3*x)*:-o.x) 

1 TO . 3*T4**4  •. 

2T0.a0S343ao*x3**2 
3-0. T4330384*x3**4 
4 TO . £13001 1 ’«x3**a 
3-0. J42al2l4*X3«*8 
a tO. 004274 ia«x 3**10 
’-0.0002484a*x3«*l2 
RETURN 

1000  X3-3..X 

FO-  0.  "»~9843a 
t-0.0000007~»X3 

2- 0 .003327 40 *X 3**2 
J-0. 000043 12*X3**3 
4+0 .00 1 372374x3**4 

3- O.OOOT£303«X3«*3 
*T0.000144?a«X3**a 

THETA«X-0.'9S349U 

1- 0.04lao34’XX3 

2- 0 . 000034344X3**2 
3»0.002a23-3*X3**3 

4- 0.000541 2S*X3**4 

3- 0. 00024333 «X 3* *5 
af0.00013333*X3**o 

2Y0«F0*3INI THETA). 3QRT.X) 

RETURN 

2000  2Y0» - 1 . £30 

RETURN 
END 


FUNCTION  2Y l i X ' 

IF.X.uT.l.E-SO)  30T0  2000 
IF  . X . 3T  . 3 . ) 30 TO  1000 
\3«x  3. 

2T 1-2.  3. I 4i342s34*X*4L0Gv0 .5<X ) «2 j 1 i X ' -0 . a3aa 1 40 
l tO . 22 1 204 i »\3**2 
2*2 . I a82'"09*X3**4 
j-i  .oi.j-ao'xxo*** 

4 +0 . 3 1 2 345 1 * x3**8 
3-0.04004'a*X3**lO 
at0.0027873»X3**l2 
2'r  1 »2V  l ■ X 
RETURN 

1)00  X3»3 . X 

Fl«0. "4 '9843*  TO . 000001 1*  *X3 
1 f0.01*34aaT*x3**2 
2t0.0001'103*X3**3 
1-0. 0024451 > fx3**4 
4T0.001 1 3a53*X3*  «5 
3-0 . 00020033  *X3  »*a 
’HET4-X-2 .35*1  ’444 
If). 12444o;;»X3 
2 TO . 00003*30  *X 3* *2 
■-0.00*3’9'4»x3»*3 
4 to. 000 ’4343* X3« *4 
ITO.OOO'4<324*x3«*5 
*-0 . 00024 1 *o«x3*t* 

2X 1 «F 1 «SIN» TSE’A'  3GRT . \ N 
RETURN 

2000  2V 1 ■ - 1 . *30 
RETURN 
END 


THIS  PAGE  IS  BEST  QUALITY  FRACTICABld| 

fro*  oopy  njKaisHjfiD  to  ooq 

GUBK'QU  TINE  CHOI  v 0*01 -OYOC • X > 

INTEGRATION  OF  ,0  AND  *0  FROM  0 TO  X 
EDWARD  R.  fUNG  AND  DONG  NGUYEN 
AUGUST  19”9 
COMPLEX  Cj.COMP 


REAL 

DATA 

Av7).B(7).Ci3>-0l3> 
A - B - c . 0/ 

» FN ( 1 1 ) • GN ( l<) 

. 10000000000E+01 . 

- . 0049999-0006+01 • 

. lCa3oC080006+01  - 

- • 3 1 o33ao0000E+00 • 

. 444479000006-0 l > 

-.  i94440oooooe-oo- 

. Cl OOOOOOOOOE-03  > 

,3o74aa91000EY00» 

. *0539 iaoOOOE +00  - 

- . 7 43S0384000E+00 • 

. 03300  U'OOOErOO- 

-.40610 1 40000E-0 1 . 

. 40791 aOOOOOE-OO  > 

-.048460000006-03- 

.60334730400E-01 • 

-.404033390006-00- 

. 100898700006-00- 

-.336616900006-03- 

■ 39923030000E-03 » 

- . 073S0370000E -03  • 

. 107003900006-03- 

-.060480000006-04. 

, '9  ■’984560006  + 00  • 

-.  1 036 4 0 403006-0 1 . 

. 1 '9709440006-00- 

-.67401 4 90000E-03- 

. 41006'60000E-03- 

DATA 

-.034393300006-03- 
FN - GN/ 

. 1 lO709«OOOOE-O3- 

- . ccsc  jaooooot-04, 
-.o 79 134031 09E+OO • 

. O4331304a37E+Ol - 

-.603813309016+01  - 

. 10991343974E+0C - 

14633o393G0E+00. 

. 140'9864404E-00. 

-. .0040113339E+00. 

. 496171 15oGaE+01 - 

-.  Ia33993 06396+0 1 * 

.302045753136+00- 

-.C973733914'E-01 . 

.044133863796+00- 

,30634O9”0l”E+O0. 

-.G4-a0G10a3»EY0l. 

. 37156168U9E+01  - 

- .Gal  431 ’60806+01  • 

.903098073a'E+0! . 

-.  666660600” 1 E +0 1 • 

.34O">460314o£  + O1  • 

DATA 

1 13031050966+01 . 
PI.CJ/3. 141390a33a. 

.0311C4aC40”E+00- 
.0. - 1 . '/ 

-. 00936 130099E-01 

3joi*zyoi=o. 

:f<x.le.o.)  return 

XT*DEN»l . 

IF( X . GT . 3 . > GOTO  1 000 
con»;., pi 
XL0G»AL0G(X/0.  ) 

X C * X 4 X / 9 . 

DO  100  N*t.’ 

TEMP-XT.  DEN 
Z JO  I =Z  jO  I +A  ( N ' «T£MP 

CYOI  *ZYO  I +CON*A  i N ) «TEMP«<  XLGG-1  . . DEN ' +B l N ' * TEMP1 

xt*xt«xo 

DEN=DEN+0 . 

100  CONTINUE 

c„oi=x*zjoi 

:yo:»x*:yo[ 

RETURN 

1000  IFiX.GT.3.'  GOTO  3000 
x3»X  3.3 
F = o*0  » 

DO  COO  N*l-U 
r*F+FNlN l «X* 

G*G+GN\ N l « < T 
COO  XT  *XT«X5 

CjOI=1 . +F«COS IX  > +G*StN<X> 

CYOI*  F«SIN<X> -G«CGS  < X ) 

RETURN 

3000  \3*8.  x 

X0*X8*X3 
DO  300  N*  1 .a 
cjoi«:jc i pc  > n i «xt 

ZYOt*CYOC+D<N! «XT 
XT  *XT*XC 
>00  CONTINUE 

COMP*CmFL-<  . <3«Z-,0I . CY 0 1 ' 

COMP*  l . -CCMPxCZ  <F  Cw  X i <-P  I J.  3JR7i<' 

C JO  I “REAL  COMP  ' 

CYO  I *4  I MAG  > 0 7MF  ' 

RETURN 

End 
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FUNCTION  ZNOI(X) 

C INTEGRATION  OF  NO  FROM  0 TO  X 

C 

C COWARD  K.  TUNG  AND  DONG  NGUYEN 


AUGUST  t9-’3 

REAL  A<7> .8(7) »C(7) »D<  11) 
DATA  A.B.C/ 

. looooooooooe+oi » 

. 1Z067492000E+01  • 

. 43813000000E-00.  - 

.Z30o973a000E+00. 

. 10730000000E-03. 
-.111 900890 OOEtOO* 

. 41743400000E-02.  - 

DATA  0/ 

. 93449344234E  + 01 . - 

-.3860a846431EYOC> 
.70401397713E+O1,  - 

ZKO I -0 . 

IF<X._c.O.>  RETURN 
0EN»AX»BX»CX»XT»1 . 

IF  < X . GT . 2 » ) GOTO  1000 
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ABSTRACT 


The  paper  studies  the  asvmptot ic  solution  of  the  surface  magnetic 

field  due  to  a magnetic  dipole  on  an  infinitely  long  cylinder  whose 

radius  is  large  in  terms  of  wavelength.  Starting  from  the  exact  modal 

series  solution,  we  extract  a dyadic  Green's  function  for  the  magnetic 

field  which  is  valid  for  all  points  on  the  cylinder.  In  particular, 

-1/2 

our  solution  justifies  for  the  first  time  the  (ks)  behavior  of  the 
field  propagating  along  the  generator  of  the  cylinder,  where  s is  the 
distance  between  the  dipole  and  the  observation  point. 
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1.  introduction 


Let  (r,  z)  be  cylindrical  coordinates.  An  infinitely  long, 
perfectly  conducting  circular  cylinder  is  described  by  r = a.  A tangential 
magnetic  dipole  with  magnetic  point  source  current  density  M is  located 
at  ?(r  = a,  > = 0,  z = 0)  on  the  surface  of  the  cylinder.  A time 
dependence  exp(juJt)  is  assumed  throughout  and  suppressed.  The  problem  is 
to  determine  the  surface  magnetic  field  H at  an  observation  point  Q = (a,p,z) 
on  the  surface  of  the  cylinder  under  the  assumption  that  ka  is  large 
(k  = 2tr/\  is  the  wave  number).  P and  Q are  connected  through  a surface  ray 
(geodesic)  which  makes  an  angle  9 with  the  f-direction;  the  distance  from 
P to  Q is  denoted  by  s;  see  Fig.  1.  The  present  high-frequency  diffraction 
problem  was  studied  by  Chang,  Felsen  and  Hessel  [1]  and  by  Lee  and  Safavi- 
Naini  [2],  however,  their  results  differ  in  various  ways.  In  particular, 

Lee  et  al.  [2,  Eq . (2.13)]  predict  a rather  peculiar  term  in  the  approximat ion 
of  the  component  H (Q)  when  9 = tt / 2 , i.e.,  when  P and  Q lie  on  the  same 
generator  of  the  cylinder.  This  term  behaves  Like 

— — exp(-jks) 
kav  ks 

for  large  ks  and  is  introduced  in  a rather  arbitrary  manner.  It  is  the  aim 
of  this  report  to  clarify  the  appearance  of  this  peculiar  term  and  other 
points  in  which  the  asymptotic  solutions  [1]  and  [2]  differ. 

We  shall  start  from  the  exact  modal  solution  for  the  surface  magnetic 
field  as  presented  in  [1].  Then  the  quotients  of  the  Hankel  function  and 
its  derivative  are  replaced  by  a Debve-tvpe  approximation.  As  a result, 
we  find  a two-term  approximation  for  the  surface  magnetic  field.  The 
leading  term  is  equal  to  the  so-called  planar  solution,  that  is,  the 

1 


Figure  1.  A surface  ray  (geodesic)  from  source  point  P to  observation 
point  i)  on  an  infinitely  long  circular  cvlinder. 


solution  for  the  surface  field  due  to  a tangential  magnetic  dipole  on  a 
flat  ground  plane.  The  present  approach  differs  from  that  of  [1]  where 
the  Hankel  functions  are  replaced  by  their  uniform  asymptotic  expansions 
in  terms  of  Airy  functions.  Then  the  resulting  approximation  for  the  surface 
magnetic  field  is  expressed  in  terms  of  Fock  functions.  The  Debve-tvpe 
approximation  of  the  Hankel  function  is  discussed  in  the  next  section. 
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2.  DEBYE- TYPE  ASYMPTOTIC  EXPANSION  OF  THE  HANKEL  FUNCTION 

(2) 

Consider  the  Hankel  function  H ( (z)  with  both  v and  z being  large 

and  positive.  From  Watson  [3,  Secs.  8.4  and  8.41]  we  quote  the  Debve-tvpe 
asymptotic  expansions,  taking  into  account  two  terms: 

H(2)(v  sech  ci)  - - j Y (V  sech  a) 


= -3 


, • e 


v (a-tanha) 


[ i - 3 coth  ot-5  coth  ci  ^ ^ ^ 1^  j 


and 


(2) 

'(v  sec  3)  = 

0 


»'2TVtanha 
a > 0 


v( tan  3 - 8)+  Jt/4 


24v 


(2.1) 


[i  + j + o(4), 

-v  v“ 


v 2ttv  tan  3 

8 > 0 (2.2) 
the  former  approximation  applies  when  the  argument  is  less  than  the  order, 
while  the  second  approximation  applies  when  the  argument  is  greater  than 
the  order.  Replacing  v sech  a and  v sec  3 by  z,  we  find 


Hl3)(z)  > 

V v TT 


1 /V, 


n 


exp  [v  cosh  (— ) - v\>  - z 

~2  21/4 

(v  - z ) 


[i  + 


2v“  + 3z‘ 


0 T 1 / O 

24 (v  - z~) 


Z < V 


and 


(O)  /z 

t (2)  - fz 


exp  [-j/z  - V + jv  cos  (— ) + “] 


✓ TT 


, ,2,1/4 

( z — V ) 


[1  + j 


2v“  + 32 


+ o(-7 
(2.3) 


TT y + O(-^r)  ] 


j / _ 


24 (z  - v-)- 

z > v . (2.4) 

(2)  1 

It  can  be  shown  that  the  asymptotic  expansions  of  (z)  are  obtainable 

from  (2.3)  and  (2.4)  by  a term-bv-term  differentiation;  see  Abramowi tz-S tegun 


XA 


[4,  pp.  366-367].  Through  division  we  obtain  the  asymptotic  expansions  of 
the  quotient  H^”\z)/H^  (z),  viz., 


/ 2 - T 

/v  z 


! (v 


9 

z“) 


2 + 0(~} 

V 


and 


h<2)U) 


JZ 


Z < V 


(2.5) 


H. 


(2)' 


(z)  /T 


f Y-T  + 0(4)  , z > V 

2(z  - v-)“  v“ 


(2.6) 


These  asymptotic  approximations  are  essentially  the  same.  Notice  that  for 

(2)  (2)' 

fixed  z > 0,  the  quotient  (z)/H^  (z)  is  an  analytic  function  of  V in 

(2)  ' 

the  whole  complex  v-plane  except  for  poles  at  the  zeros  of  (z) . These 

zeros  lie  in  the  second  and  fourth  quadrants  of  the  v-plane  and  are 
approximately  given  by 


v = ± [ z + e 2^/3  (z/2)^^  + 0(z  '*‘4  ] 

P P 


p = 1,2,3,...  (2.7) 


where  -q  are  the  negative  zeros  of  the  Airy  function  Ai(x),  i.e.,  Ai(-q  ) = 0; 
P P 

see  Keller,  Rubinow  and  Goldstein  [5].  In  order  to  uniquely  define  the 

fl  T 

square  root  /z  - v , we  introduce  branch  cuts  from  v = z downwards  and 
from  v * -z  upwards  in  the  complex  v-plane.  Then  it  is  easily  seen  that 


n T /2  T 

/z  - v in  (2.6)  passes  into  -j/v  - z“  in  (2.5),  when  v passes  above 
the  branch  point  z or  below  the  branch  point  -z.  This  shows  the  equivalence 
of  (2.5)  and  (2.6).  In  the  sequel  we  shall  use  the  approximation  (2.6), 

2~  /9  2 i i 

whereby  it  is  understood  that  /z  - v -*■  -j/v  - z when  ]v|  > z. 

The  approximations  (2.5)  and  (2.6)  apply  when  either  [ V ' , or  z,  or 

ii  -2-2 

both  vj  and  z are  large,  the  error  being  of  order  v " or  z , whichever 
is  smallest.  The  approximations  are  not  valid  in  the  transition  region 


r 


|v  i z|  < 0(z2/3)  and  z large.  Furcher  difficulties  appear  when  both  v and 

(2) 

z are  small.  For  example,  from  the  power-series  expansion  of  (z). 


one  may  find 

<2)(z) 


»r<*> 


log 


4 2)  + (Y  + 4?> 2 + 0(zj  log2  z)  , z - 0 (2.8) 


where  y denotes  Euler's  constant.  The  latter  result  does  not  agree  with 

(2.6)  when  V = 0.  On  the  other  hand,  when  v = 0 and  z is  large,  we  find 

(2) 

from  the  asymptotic  expansion  of  (z): 


-JTy = J ~ TT  + 0(— ) , z ~ » 

Hq  ' (z)  "Z  z‘ 


(2.9) 


which  is  in  perfect  agreement  with  (2.6)  when  v = 0.  Similarly,  when  z 

is  small  and  |vj  is  large,  we  may  find  from  the  power-series  expansion  of 
P) 

H 4 (z): 


H^2)(z) 


J - + 0 (~7  ) , 

V 2v3  V4 


(2.10) 


in  agreement  with  (2.5). 

In  subsequent  sections  the  approximation  (2.6)  is  going  to  be  used 
in  integrals  of  the  form 
.(2), 


r®  , . , * H. , (k  a) 

L »„  <V> 


(—•ti) 


where  J)  > 0,  k = /k  - k^,  k^  is  real  when  k~  < k“,  and  kt  negative 
2 2 

imaginary  when  k"  > k~.  In  order  to  justify  the  replacement  of 
(?)  (?)< 

H * (k  a)/ H (k  a)  by  the  approximation  in  (2.6),  we  distinguish  three 
V C V t 

cases:  (i)  If  k_a  is  real  and  large,  the  integration  contour  in  (2.11)  is 

6 


W 


deformed  by  introducing  semicircular  indentations  above  k^a  and  below 

1/3 

both  of  radius  0((k^a)  ).  Along  the  deformed  contour  the  approximation 

(2.6)  is  certainly  valid.  After  making  the  replacement  of  the  integrand 
by  (2.6),  the  contour  is  deformed  back  to  the  real  axis.  Notice  that  in 
both  deformations  no  poles  or  branch  points  of  the  integrand  are  crossed. 

(ii)  If  k^a  is  negative  imaginary  and  large,  the  approximation  (2.6)  is 
valid  along  the  real  v-axis  and  no  deformation  of  contours  is  needed. 

(iii)  If  kca  is  small  due  to  k_„  = k,  then  the  approximation  (2.6)  is  not 

valid  when  |v|  is  also  small;  it  is  valid  though  when  jv|  is  large.  It 
is  not  clear  what  effect  this  will  have  on  the  error  in  the  resulting 
approximation  to  the  integral  (2.11).  This  case  certainly  needs  further 
consideration.  For  later  use  we  also  establish  an  approximation  for 
the  quotient  (z)/H^  ;(z): 


a. 


«?><*> 


H>  T 
J v _ 


2(z"  - v") 


(2.12) 


rr 


rr 


where  it  is  understood  that  v'z  - V-  ■+■  -j/v  - z~  when  |v|  > z. 


3.  SURFACE  MAGNETIC  FIELD  DUE  TO  A CIRCUMFERENTIAL  MAGNETIC  DIPOLE 


(2 

i.I  Magnetic  field  component  H . 
2_ 

In  the  case  of  a circumferential  dipole 


M 


(3.1) 


the  resulting  surface  magnetic  field  components  are  denoted  by  H (Q) 

9 

and  H3(Q).  Neglecting  the  contribution  of  creeping  waves  which  have 
travelled  around  the  cylinder,  it  is  found  in  [1,  Eqs.  (18)  and  (19)]  that 


-jk  j k'  f* 

3 ! dk  e 2 

,2  3 l 3 z ,3 

4t r cju^a  ; k 

0 — » t -°° 


, -is*  2 

dv  e v 


IT  ' (k  a) 
v t 


ik" 


-jk  z , f°o 
! z 1 

— “ j dk2e  r- 

“V  loo  t 


dv  e 


. „H(2),(k  a) 
-jv<J  v t 

Hf(ka) 

V L 


where  k 


(3.2) 


■■>  9 9 

k”,  k is  real  when  k“  < k“,  and  k negative  imaginarv  when 
z t z t 


2 7 

k~  > k”.  If  necessarv,  one  mav  think  of  k having  a small  negative  imaginarv 
z 


~ 

part.  Then  k = /k  - k~  has  branch  cuts  from  k = k downwards  and  from 
t z z 

k * -k  upwards  in  the  complex  k -plane, 
z z 

In  (3.2)  we  replace  the  quotient  of  the  Hankel  function  and  its 
derivative  by  the  approximations  in  (2.6)  and  (2.12).  Furthermore,  we 


set  v = k^a  in  the  inner  integrals  in  (3.2).  Then  we  are  led  to  the 
following  approximation  for  h'T(Q): 


8 


rMa^M’^rifiiinrt'i’ni  iii^Ti  ' ’-f/rmi  - *■—- - ^-  ■ 


fc-r 


H‘(Q)  * -5J—  ffexp  Hk  z - Jk  Ml  { ,jKt 

4tt  alUrt  ' ‘ |<  n T 

0 -«-■»  t /k  - k 

^ tv 


/ o o “* 

,2  j /k“  - k“  , k 
k ■ t v 1 t 

{ ^ — - + - — r^— — > dk  dk 

t *t  J 2(k;  - k~)  ' z 

t v i 


Y f 


I -A ; 

3 2(k:-k->“ 

t V 


exp  [ — j k z-jk^aA] 


2 ? 4 o o o o oo 

k - k"  . k - k“k“  - k'k~  - k‘k~ 

V j v z v z 

{2 ~2 ~2  (k  - k"  - k ) 

► k - k - k v 2 

v 2 


dk  dk 

V 2 


(3.3) 


0 0 o 


where  Y * Jz  I\x  ; the  square  root  vk  - k“  - kw  is  positive  when  k“  + k“  < k“ 

oo  v z v z 

2 2 2 

and  negative  imaginary  when  k“  + k~  > k~.  Notice  that  (aj>,z)  are  just  the 

y z 

rectangular  coordinates  of  Q on  the  developed  cylinder. 

The  result  in  (3.3)  can  be  expressed  in  terms  of  the  derivatives  of 
the  following  two  key  integrals 


r=°p°  dk  dk 

I,  (a<J>,z)  = i exp  [ — j k_  z - jk  a<j)]  - - 

1 J 2 y p 2 

-CO-CD  /k“  _ k”  - k“ 

v z 


(3.4) 


rr 

I-(aJ>,z)  = | ; e> 


f^f00  dk  dk 

i j exP  Hk,2  - Jk  ■ a*]  — , -V  T — , 

LL  " y (k-  - k;  - k-) 


(3.5) 


both  of  which  can  be  evaluated  in  closed  form.  To  that  purpose  we 
introduce  polar-coordinate  variables  k^  = t cos  oc,  k_  * t sin  a, 
ai  = s cos  d,  z = s sin  0.  Then  by  use  of  Watson  [3,  Eqs.  2.3(1)  and 
13.47(4)],  1^  reduces  to 

I,  (a4) , z)  = f t dt  f2'  d,  ^^i.-jgt  cos  jot-e ) ] _ r j /gt) — 5 jt 

i J J /-i r J 0 — - 

0 0 *k“  - t 0 ► t“-k“ 


9 


Mo t ice  that  we  made  the  replacement  /k  - t -*•  -j/t  - k , and  that  in  the 
latter  integral  the  path  of  integration  passes  above  t = k. 

In  a similar  manner  I9  reduces  to 


I^CaO , z)  = J t dt  j da 


f"  cos  (,-e>i  , 2,  |“  j (st, 

i (k  ' tV  J,  (c"-k2)- 


(3.7) 


According  to  [3,  Eq.  13.6(2)]  and  [4,  Eq.  11.4.44],  we  have 


t J (at) 

i — 7- 7-7  dt  = ~ K , (az)  = j-  K (az) 

l (t2  + Z2)2  2z  -1  22  1 


(3.8) 


valid  for  a > 0,  Re  z > 0.  Remember  that  k has  a small  negative  imaginary 
part,  if  needed,  hence  Re(jk)  > 0.  Setting  a = s,  z = jk  in  (3.3),  we 
obtain  for  I„: 


I2(a<J>,z)  = KL(jks)  = - H<2)(ks) 

by  use  of  [4,  Eq.  9.6.4]. 

Returning  to  the  result  (3.3)  for  H*" (Q) , we  have 


(3.9) 


»:<Q>  * - -V  ih2  * -S]  V-*.*  - j-A±‘  + k2  ~ 2 + ^ 


4tt  k ! v.  3 (ad) 


B(a<p)' 


Bz2d(a4>)y  2 


l2(a*'z)  = 2^jk]f2  + cos2  8 


32  . .21  A.e‘ji 

— — 2 -r  sin 


^ + ^(5  + * &)  • 


2 2 3 2 2 1 ■ 

cos  9 sin“  9 — 7 — (1  - 6cos  9sin~9)—  - 

3s4  3 ; 


+ (2  - 15  cos29  sin2  0)4r  - (2  - 15c 


2 2 1 3 ( (?) 

os  3 sin  9)—  tt— ) (sH' (ks 


f 


The  derivatives  of  the  Hankel  function  H^-/(ks)  can  be  evaluated  and 
simplified  by  means  of  the  well-known  recurrence  relations  for  Bessel 
functions*;  see  e.g.,  [3,  Sec.  3.2]  and  [4,  Eqs.  9-1.27  and  9.1.30]. 

c 

Thus  we  obtain  as  our  final  result  for  H.(Q): 

0 

->  -jks 

«*«»  i*1"' 6 + i<2  - 3sl":  *> + ~h  (2 ' 3311,2  9)1 

k s 
o 

k~Y  (2)  1 (2)  (2)  2 2 

- [3Hq  (ks)  - ~ H^ (ks)  + ksH^  '(ks)  cos  9 sin“  9]. 

(3.11) 

Notice  that  except  for  the  Debye-type  approximation  to  the  Hankel  function 

quotients,  no  further  approximations  were  involved  in  the  derivation 

of  (3.11).  The  Hankel  function  H^~y(ks)  in  (3.11)  can  be  expressed  in 

p)  (2) 

terms  of  Hq“  and  H^  through 


(2)  (l)  (!)  3(1) 

ksH,  (ks)  = -4H^  '(ks)  - ksH,  ' (ks)  + p-H,  '(ks) 
3 U 1 ks  1 


(3.12) 


The  first  term  in  (3.11)  is  exactly  equal  to  the  planar  solution,  that  is, 
the  solution  for  H due  to  a magnetic  dipole  M = on  a flat  ground  plane; 
see  [1,  Appendix  D]  and  [2,  Eq.  (2.18a)].  The  second  term  in  (3.11) 
represents  the  effect  of  the  finite,  but  large,  radius  of  curvature  of  the 
cylinder.  For  large  ks  the  Hankel  functions  in  (3.11)  can  be  replaced  by 
their  large-argument  asymptotic  expansions,  thus  leading  to 


( 2 ) 

* In  (3.10)  the  differential  operator  in  front  of  sH^  (ks)  can  be  rewritten  as 

4 222  2 23  2 44  2 2 

k +ksD  + 2k  D - s D -D  -sD  cos  A sin  9 

. n 1 3 

where  D = — — . 

s 33 


11 


1 


h*(q>  » r e 


.K  ,-;ki  r 9 


2mj  ks 


sin-  0 + ^ (2  - 3sin-  0)  + — 


- (2  - 3s in-  3) 


k-s 


1 t 1/2  -Jit/4,.  1/2  r-  2 _ 2 ....  35,  . ..  1 

+ Tj^(— ) e (ks)  'i3  - cos  3 sin  0(jks+— ) + 0(— ) ■ . 


ks 


(3.13) 


Lee  the  term  of  order  7 — in  (3.13)  be  denoted  bv  W,  then  for  8 = t/2 

ka 

one  has 


..  3 or.  1/2  -Jir/4 .1/2 

W “ 4kl  e (k3) 


(3.14) 


which  is  in  exact  agreement  with  the  peculiar  term  in  Lee  et  al.{2,  Eq.  (2.18c)] 
We  shall  now  compare  our  approximate  results  in  (3.11)  and  (3.13)  to 
the  solutions  presented  in  [1]  and  f2].  Chang  et  al.  [1]  have  two  different 
formulas  for  H^,  namely,  the  asvmptotic  formula  [1,  Eq . (124)] 


, 2..  -jxs 
„c,„,  k \ e 

VQ)  “ 


r . - ,,  , . ; , 2 , 61  . 2 1 sin-6v_  ,.-N 

[sin  0 vrt(0  + £*r(cos  9 _ W4Sln  ~ g ~yQU) 


0 


¥:  VC)] 

cos  0 


(3.15) 


not  valid  when  r .gets  close  to  t/2,  and  the  "full  formula"  [1,  Eq . (130)] 


<(h»  = r e 


? -j  ks 


, [sin-  8 vn(0  + T2r(2  cos-  0 - sin-  8 - l)v  (5) 


. sin-  0 v, (l)  - u^(i;) 
ks 


(3.16) 


COS~  0 

which  remains  finite  on  9 = t/2.  Lee  et  al.[2]  present  the  formula 
[2.  Eq.  (2.16b)] 


+ On  extending  the  result  in  (3.13),  it  is  found  chat  the  term  of  order  0(~) 
is  equal  to  KS 


i 5 945  22  1 

ks'-  8 + TVJ005  J sin  3)  + 0i  — — -) . 

k-s- 


(3.13a) 
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f 


. 2 -jks 

H^(Q)  * ~ ~— 

$ 2irj  ks 


sin  0 v(0  + ^-(cos‘  0 - sin-  0)v(£)  + -L  cos-  0 u(£) 
ks  ks 


-1/3 

1 2 2 2 1 4/3  a 

+ ~ 9(2  cos  6-  sin  9)u(£)  + -^-pr  cos  9{sin-  9 v ' (£) 

k-s-  '■ 


(ka) 


+ &-iny  u'(£)  + r^~  cos-  9 u'(£)} 
cos-  8 


(3.17) 


In  (3.15)  - (3.17), 


.-1/3  ks  4/3  „ 

= 7 cos  e 


(k  a) 


2/3 


and  v = v,  11  = u stand  for  certain  Fock  functions  as  defined  in  [1] 

0 0 


and  [ 2]  . 


We  shall  now  re-expand  the  solutions  (3.15)  - (3.17)  in  the  case  of 


large  ka,  up  to  and  including  order  — . For  large  ka,  £ is  small  and 

Ktl 


we  replace  the  Fock  functions  by  the  leading  terms  of  their  small-argument 
expansions  quoted  from  [1]  and  [2],  viz., 


r 


<r,  _ ,fv  „ , /TT  Jir/4r3/2 
v0(O  - v(>)  * 1 - — e £ 


, Sir  jTr/4_3/2 
u (£)  = u(£)  - 1 - — eJ  t. 


L /n  j7T/4.3/2 

1(i.)  - 1 + — eJ  £ 

3 /if  j it/4.1/2 


^ v'(4)  = - g 


4 


3.18) 


Then  the  results  (3.15)  and  (3.16)  of  Chang  et  al.  [1]  become 


2 -jks  r 

f.C(Q)  = S. 1 

? 2-rj  ks 


sm 


2 ^ . i /.  >95  ,2,8  sin-  3, 

8 + ks(2  ‘ 24  Sin  8 + o — 


COS 


L -JTT/4  v 1/2  r ~ 10  . 2 a 2 ' 2 w..  . 35. 

+ — — ( r)  e (ks)  ( 3 - •—  sin  0 - cos  * sin  (jks  + ~) 

-*Ka  2 9 2-f 


(3.19) 
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i 


and 


„c  ,,  A «-|ka 
!,0l!>  ‘ ST  kT 


jin-  9 + r^-(2  - 3sin~  0) 

K.S 


1 tt  1 / ° - \ 7 / 4 l/"  1 > > 

+ -rf-(-r)  -e  J tks)  3 + sin-  0 - cos"  6 sin"  J ( Jks  + 
-*  ka  J 


(jks  + 3)} 


(3.20) 


k-Y  e 


-jks 


1 


H- (Q)  a —•  ~ — r — ,s in-  9 + r*-(2  - 3sin2  0)  + — ~(2  - 3sin-  0) 

>J>  .’Tj  ks  ^ <s  k--'- 

+ T7-(^)1/:e*J,T/4(ks)l/2{3  - cos2  0 sin2  3 (jks  + ) - rLcos2.M  1- 3sin20) } 

■tka  2 2 ks 

(3.21) 

The  asymptotic  formula  (3.19)  due  to  Chang  et  a 1 . [ 1 ] agrees  with  (3.13) 

only  at  0 = 0.  Even  the  leading  term  in  (3.19)  agrees  with  the  planar 

solution  only  at  0 = 0.  The  "full  formula"  (3.20)  of  Chang  et  al . and  the 

solution  (3.21)  of  Lee  et  al.  do  have  the  planar  solution  as  their  leading 

term.  The  nest  term  of  order  . — in  (3.20)  is  in  agreement  with  (3.13) 

ka 

only  on  0=0  and  certainly  not  on  6 = t/2.  Ihe  corresponding  term  in 

(3.21)  agrees  with  (3.  13)  both  at  9 = 0 and  at  = i/2,  however,  in 
between,  the  agreement  is  only  partial. 

3.2  Magnetic  field  component  h\ 

In  this  section  we  consider  the  z-component  H of  the  surface  magnetic 
field  due  to  a circumf erent ial  magnetic  dipole  '1  * *>.  The  contribution 
of  creeping  waves  that  have  travelled  around  the  cvlinder  is  again 
neglected.  Then  according  to  [1,  Eq . (20)),  is  given  by 

- j k 


H (Q)  = sgn( f)sgn(z) 


sTT  <ou  a 


d k e 


14 


-Hi!  (0-\v> 

Hj-'V 


(1.22) 


r:  A* 


where  k^  = /k“  - k~  as  before.  In  (3.22),  Che  quotient  of  the  Hankel 

function  and  its  derivative  is  replaced  by  the  approximation  (2.6), 

and  we  set  V = k a.  Thus  we  obtain 
v 


r® 

H"(Q)  - — — I j exp[-jk  z - jk  ao] 

4TT~k  J j 2 y 


k k . k k (k"-k~)  I 

y 2 ,i  v z z \ ,,  ,, 

+ TT  i f dk  dk„  • 


/ 2 2 2 
/k  -k  -k 
v z 


(k~-k“-k~) 
v z 


(3.23) 


The  latter  integral  can  again  be  expressed  in  terms  of  the  derivatives  of 
I^(avp,z)  and  I?(a4>,z),  as  defined  by  (3.4)  and  (3.5),  viz., 

r 1 2 9 \ 

•>  * - 717  slfuio  V**-2*  k'  +fi  ¥**•«>  ■ <»•**> 

-*TT  k 3z  I w 

On  substituting  the  explicit  values  of  and  I0  from  (3.6)  and  (3.9),  we 


,.c  . . V (T  3“  1 A fe"1ks 

Hz(Q)  - 27jk  cos  e sin  ehT2  - 7 TsJ  |~1- 


it  r ’22 

+ 8k7  fk  S D + sin 


' 3sV  + 3s2D3}(shJ:) (ks)) 


(3.25) 


where  D = — t— . Bv  use  of  the  well-known  recurrence  relations  for  Bessel 
s )S 

functions,  we  ultimately  obtain 

..  k“Y  e ^ ks  , . , r,  3 ] 3 , 

Hz(Q)  - - 27j  cos  a 3in  e [1  - kl  - r^1 


cos  3 sin  3 [H.,”\ks)  - ksH^"\ks)  cos"  3’ 


(.3 .26^ 


In  deriving  (3.26)  only  the  Debye-type  approximation  to  the  Hankel  function 
quotient  is  involved.  The  first  term  in  (3.26)  is  again  the  planar  solution 
for  H^Q)  as  given  in  [1,  Appendix  D].  For  large  ks  the  Hankel  functions 
in  (3.26)  can  be  replaced  by  their  large-argument  asymptotic  expansions. 


vielding 
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. 2„  -jks  r ,,  a 

„c , ,,  k \ e % , , i , 3 1 3 

Hz(Q)  = - 2ij  ksT  °0S  6 3in  0 I1  " 'k sTTz 

J U k S 

+ 7T-(?)I/2e"J,7/4(ks)1''",{l  - cos2  0(jks  + ~)  + 0(~)} 


(3.27) 


On  extending  the  result  in  (3.27),  it  is  found  that  the  term  of  order 


0(— ) is  equal  to 

1 . 15  945  2 a.  , 1 . 

i(-  T + 123  003  9)  + °(^2) 


(3.28) 


£ 

Notice  that  the  approximate  result  (3.26)  for  H^(Q)  vanishes  when  6=0 

or  6 * tt/2.  The  same  holds  true  for  the  exact  value  of  H*"(Q)  in  (3.22). 

z 

The  present  approximate  results  in  (3.26)  and  (3.27)  are  now  compared 
to  the  solutions  derived  in  [1)  and  [2],  Chang  et  al.[i]  present  the 


asymptotic  formula  [1,  Eq.  (125)] 


,2  -jks 


cos  0 sin  0 (v  ( 


..  , 1 , 23  , 5 sin-  0.  ,rx , 

O + ±(-  T + s r-)vn(5)l 


0W/  ks  8 9 


COS  0 


(3.29) 


and  the  "full  formula"  [1,  Eq.  (105)] 


, , -JP.S  , 

He(Q)  = - cos  0 sin  0 [v  (O  “ 7^  v ( E]>  ] 

z 2ttj  ks  0 ks  0 

From  Lee  et  al.[2,  Eq.  (2.6)]  we  quote  the  solution 


(3.30) 


' (Q)  = - cos  0 sin  0 [H  (Q)  - Hc (Q) ] 


. 2 -jks 
k Y e 

2^j"  ks 


,-1/3, 


— cos  0 sin  0 (1  - fl)v(Z)  ~ + -y  — )u(C) 

L k s 

2 “N 

s4/J  9{v’(£)  + -3'i-n-— 0-  u'(S)  - - u’(-)  • . (3 

COS~  0 ‘S  _ 


(3.31) 


The  “olutions  (3.2Q>  - (3.31)  are  re-expanded  for  Large  ka  by  means  of  (3.18). 


Then  the  results  (3.2q)  and  (3.30)  of  Chang  et  al.(i]  become 
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7 - i ks  r 9 

“>>  * - sf  ns-  -■°8  s si" 9 1 + ii(-  f + 1 si!T7) 

w COS  0 

1 , IT. 1/2  -Jit/4,.  . 1/2,5  2 ....  247,? 

+ 4ka(I)  e (ks)  " cos  e^Jks  + “72)j 


and 


-jks 


f 


(3.32) 


2ttJ  ks 
{-  cos”  6 (jks  + 3)} 


ks  4ka  2 


(ks) 


1 /•> 


(3.33) 


whereas  the  result  (3.31)  of  Lee  et  aL[2]  becomes 


c k2Y  P~jks 

H ,(Q>  ” ~ Irr  cos  e sin  e 

z 2irj  ks 


1 - - 3 
ks 


•">  o 

k s 


+ — -(-^)1/”e_j'T  '4(ks)1/2{3  - cos2  0 ( j ks  + ~)  + ^ cos2  0} 


(3.34) 


The  asymptotic  formula  (3.32)  agrees  with  (3.27)  only  at  0 = 0.  Both  the 
"full  formula"  (3.33)  of  Chang  et  al.  and  the  solution  (3.34)  of  Lee  et  al 
have  the  planar  solution  as  their  leading  term.  As  for  the  next  terms  of 
order  in  0-33)  aru*  3.34),  there  is  only  partial  agreement  with  (3.27). 
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4.  SURFACE  MAGNETIC  FIELD  DUE  TO  AN  AXIAL  MAGNETIC  DIPOLE 


In  this  section  we  consider  the  case  of  an  axial  dipole 

M = z . (4.1) 

The  resulting  surface  magnetic  field  components  are  denoted  by  H. (Q) 

Q 

and  Ha(0) . Exact  results  for  Ha(Q)  and  Ha(Q)  are  presented  in  [1,  Eqs.  (26) 

Z y Z 

and  (27)],  subject  to  a neglection  of  the  contribution  due  to  creeping 
waves  which  have  travelled  around  the  cylinder. 

4.1.  Magnetic  field  component  H . 

. £ 

The  result  for  H in  [1,  Eq.  (27)]  is  identical  to  that  for  the 

v 

component  H^,  as  given  in  (3.22).  Thus  when  applying  the  Debve-tvpe 

approximation  to  the  Hankel  function  quotient,  we  are  led  to  the  results 

(3.26)  and  (3.27)  for  Ha.  Also  the  asymptotic  formula  [1,  Eq.  (127)]  and 

* a 

the  "full  formula  [1,  Eq.  (112)]  of  Chang  et  al.  for  H are  the  same  as 

$ 

those  for  H^.  Hence  Eqs.  (3.29),  (3.30),  (3.32)  and  (3.33)  also  hold  for 
Ha(Q).  Also  the  solutions  (3.31)  and  (3.34)  of  Lee  et  al.  for  HC(Q)  hold 
true  for  fT(Q)  as  well. 

4.2.  Magnetic  field  component  H'\ 

a * ?» 

According  to  [1,  Eq.  (26)],  is  given  by 


H2(Q) 


4*  uju0a  _K 


r „„  '3v„ 

dk  e k 

J 2 e 


r 

j 

—CO 


dve 


(2) 

• j.  H ^ (k  a) 

-jvd>  v v t 


(2)  ’ 

' (k  a) 
v t 


(4.2) 


* ?3  ^3 

Minor  printing  errors:  in  Eq.  (127),  replace  -7-  by  in  Eq.  (112), 

■4  ' O 

replace  Ha  bv  Ha. 

z >> 

Notice  that  there  is  a misprinc  in  [1,  Eq.  (26)]:  -j  in  front  of  the 
integral  should  be  +j . 
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In  (4.2),  che  quotient  of  the  Hankel  function  and  its  derivative  is 
replaced  by  the  approximat ion  (2.6)  and  we  set  v = k^a.  Thus  we  obtain 


Y I"”!"”  ) k2  " k2 

H1(Q)  — ( exp  [ -j  k_z  - jk.a$]  < 

2 4TT~k  J J z y /~T~rj 

_«_oo  / k -k  -k 

V.  y z 


0 0 0 
. (k“  - kV 

J 2 

2a  2 2 2.  2 


T-J-J  -a  (kz-k"-k-) 


dk  dk  . 


(4.3) 


The  latter  integral  can  be  expressed  in  terras  of  the  derivatives  of 
I^(a$,z)  and  I0(a$,z),  as  defined  by  (3.4)  and  (3.5),  viz., 

2 2 2 

H^(Q)  ^ — V [{It2  + (a^.z)  + -Ult2  + I (a4>,z)  ] • (4.4) 

4iT"k  9z“  1 -a  3z  Z 

On  substituting  the  explicit  values  of  1^  and  In  from  (3.6)  and  (3.9), 

we  find 


a y r a a q-  213  -^—jks  r 9 9 2 7 

H (Q)  ' 9-41,  l^k  + sin-  0 — - + cos'  0 — t— }( ) + — { k:  + 2k~  sin-  0s-D- 

z 2>rjk^  ,^g2  s 3s  s aka 


4 , 4 4 


+ 2k-D  + sin  0s  D +6  sin- 0 s-D  + 3D-}(sH^“  (ks))| 


(4.5) 


where  D = — -j^-.  Using  the  well-known  recurrence  relations  for  Bessel 
functions,  we  ultimately  obtain 


„a.  k Y e 

V,!)  * 3^-e 


k“Y  e~^S  7 i i 

-r-7  — r~  [cos-  0 + rM2  - 3 cos-  0) 


+ ~-t(2  - 3 cos-  0)] 


k“Y  7 (2\  TO-)  C7-)  A 

- Y~bka^~2  cos""  ® H9-  (ks)  - :j^7H^-tks)  + ksH^-  (ks)  cos  0]. 

(4 . 6) 

Notice  that  except  for  the  Debye-type  approximation  to  the  Hankel  function 
quotient  in  (4.2),  no  further  approximations  were  involved  in  the  derivation 

of  (4.6).  The  first  term  in  (4.6)  equals  the  planar  solution,  that  is, 

\ 

the  solution  for  due  to  an  axial  magnetic  dipole  M = z on  a flat  ground 
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plane;  see  [1,  Appendix  D)  and  [2,  Eq . (2.17)].  For  large  ks  the  Hankel 
functions  in  (4.6)  can  be  replaced  by  their  large-argument  asymptotic 
expansions,  yielding 


,,  . kfv 

1 ^ 27Tj  ks 


cos'  9 + rf-(2  - 3 cos'  3)  + -~(2  - 3 cos'  0) 

. — *■ 

. K S 


+ ^(f)1/2e"j7T/4(ks)i/2{2  cos2  d - cos4  0(jks  + ~)  + 00— )} 


ks' 


(4.7) 


Let  the  term  of  order  ; — in  (4.7)  be  denoted  bv  W,  then  for  0 = T / 2 one 

ka 

has  W = 0.  Hence,  there  is  no  term  of  the  form  (3.14)  in  the  approximate 
result  for  Hd.  This  agrees  with  Lee  et  al.[2,  Eq.  (2.17)]  where  there  is 
no  such  term  either.  On  extending  the  result  in  (4.7),  it  is  found  that 
the  term  of  order  0(^)  is  equal  to 


^j(-l  - -M-  cos'  0 + cos"4  0)  + O(-mp^) 

k s 


(4.3) 


The  present  approximate  results  in  (4.6)  and  (4.7)  are  new  compared  to 
the  solutions  obtained  by  Chang  et  al . [ 1 ] and  Lee  et  a 1 . [ 2 ] . In  [1]  two 
different  formulas  for  H __  are  presented,  namely,  the  asymptotic  formula 


[1,  Eq.  (126)] 

v.2v  -Jks 


k \ e 


Hz^>  " 2rj 


i ■’0  "* 

[cos'  3 v0(O  + J-(_.  sin- 


cos'  0)v (£)]  . 


(4.9) 


and  the  "full  formula"  [1,  Eq.  (Ill)] 

, U“V  n i 1 2 

Hd(Q)  - -r— - — : [cos'  0 v (5)  + r~(2  sin'  0 - cos'  8)v  (£)]  . 

z 2TT1  K.S  0 ks  0 


2tt  j ks 

From  Lee  et  al.[2,  Eq.  (2.15b")]  we  quote  the  solution 


(4.10) 
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v(?) 


H (Q)  = Trr  — r 
z 2tj  ks 


'COS 


“(cos-  9 - sin-  9)v(?)  + . sin  9 u(^) 

KS  KS 


-1/3, 


+ — 7^1-  sin-  9 - cos-  9)u(£)  + — cos"^39{cos-  9v'(£) 


(ka) ' 


+ sin-  9 u' (^)  + r^-  sin " 9 u' (£) } 


ks 


(4.11) 


The  solutions  (4.9)  - (4.11)  are  re-expanded  for  large  ka,  that  is,  for 
small  by  means  of  (3.13).  Then  the  solutions  (4.9)  and  (4.10)  of 
Chang  et  al.[l]  become 


,2V  -jks 
k Y e 


H (Q)  = ^rf  — - 
c x 2-tj  ks 


2 a ^ i ,20  247  2 

cos  9 + r*-(— rr  cos  9) 


ks  9 72 


1 ,tt  1/2  -jTT/4,,  . 1/2 f 20  2 Q 4 a/J,  ^ 247.J 

+ 7^(~)  e J (ks)  {—  cos  9 - cos  9(jks  + -~) } j , 


J 


(4.12) 


and 


Hc(Q)  = 2^  TT~  [C0S“  6 + -(2  ~ 


KS 

i !7 


3 cos-  9) 


1 ,TT  1/2  -1TT/4  I / i „ 1 a 

+ ^j^(~)  “e  (ks)-'  -{2  cos-  9 - cos  9(jks  + 3)} 


(4.13) 


whereas  the  solution  (4.11)  due  to  Lee  et  al.[2]  becomes 

3 cos-  9)  + .,(2  - 3 cos' 


Hz(Q)  a M jL4r  !cos2  6 + ii(2  • 3 cos:  9)  + " 3 cos“  ?) 

~ J L k s" 


i /? 

. -w  *- r c. 


n. 


(ks)  16  cos- J - cos  9 (jks+— ) + ^jcos-9(2-3sin-9) 


(4.14) 


The  asymptotic  formula  (4.12)  due  to  Chang  et  al.[l]  does  not  agree  with 
(4.7),  even  its  leading  term  does  not  agree  with  the  planar  solution. 

The  "full  formula"  (4.13)  of  Chang  et  al.  and  the  solution  (4.14)  of  Lee  et  al 
do  have  the  planar  solution  as  their  leading  term.  The  next  terms  of  order 
—■  in  (4.13)  and  (4.14)  agree  only  partially  with  the  corresponding  term  in  i- 


REPRESENTATION  OF  THE  SURFACE  MAGNETIC  FIELD  IN  DYADIC  FORM 


3. 


According  to  Lee  et  al.[2,  Eq.  (2.6a)],  the  asymptotic  solution  t'or 
the  surface  magnetic  field  at  Q due  to  a magnetic  dipole  M at  P can  be 
expressed  in  the  following  dyadic  form: 

A A A \ 

H(Q)  = M • fb'bH^Q)  + t'tHt(Q)]  . (5.1) 

A A 

Here,  t'  and  b'  are  the  unit  tangent  and  unit  binormal  of  the  surface  ray 
at  the  source  point  ?,  and  similarly,  t and  b are  the  unit  tangent  and  unit 
binormal  of  the  surface  ray  at  the  observation  point  Q;  see  Fig.  2 for  a 
picture  of  these  vectors  on  the  developed  cylinder. 

In  Chang  et  al.[l,  Eq.  (128)],  the  surface  magnetic  field  is  represented 
by  a different  dyadic  form,  namely, 

\ A \ N N A 

H^Q)  - M • Ib'bA  + t'tB  + «j'?C]  (5.2) 

A A 

where  ' and  are  unit  vectors  at  P and  Q in  the  direction  of  increasing  f. 
The  result  in  (5.2)  is  based  on  the  "full  formula"  for  the  magnetic  field 
components.  The  main  difference  between  (5.1''  and  (5.2)  is  that  the  result 
of  Chang  et  a L . [ 1 ] contains  a cross  term  C,  whereas  no  such  term  is 
present  in  the  result  of  Lee  et  al.[2]. 

We  now  examine  the  possibility  of  expressing  our  results  for  the 
surface  magnetic  field  in  a dyadic  form  similar  to  either  (5.1)  or  (5.2'. 

N 

— 

Starting  from  (5.1)  in  the  case  of  a circumferential  dipole  M * f , 
the  surface  magnetic  field  has  components  H^(Q)  and  Hu))  given  by 
„ •*> 

H^iQ)  = H,  sin-  0 + H cos-  0 , (5.3' 

0 b t 

and 

H^(Q)  * - (H^  - H^)  sin  0 cos 


( 5 . ' 


r 


Similarly,  for  an  axial  magnetic  dipole  M = z the  surface  magnetic  field 

S' 


components  and  H^(Q)  become 


H"lQ)  * - iHb  - H£)  sin  ^ cos  6, 


(5.5) 

(5.6) 


,c  „c  „a 


H (Q)  * H.  cos"  9 + H sin-  9 
z ^ b t 

On  substitution  of  the  actual  values  of  H",  H'"  = H“  and  H"1  as  given  bv 

0 z v 2 

(3.11),  (3.26)  and  (4.6),  we  can  consider  Eqs.  (5.3)  - (5.6),  as  a system 
of  four  equations  for  the  two  components  and  Ht.  It  is  easily  seen 


that  this  system  is  incompatible  and  does  not  have  a solution  for  H and 


H^.  Indeed,  when  taking  the  difference  of  (5.6)  and  (5.3).  we  find 


H - H 
b t 


H*(Q)  - H*(Q) 


(5.7) 


whereas  according  to  (5.4),  we  have 


H®(Q) 

H , - H 


b t sin  a cos 

Now  it  can  easilv  be  verified  that 


(5.3) 


H*(Q)  - H*(Q) 


H (Q) 


4 i'  n < 4 


sm  a cos 


(5.9) 


cos-  9 - sin-  9 

In  conclusion,  it  is  not  possible  to  express  our  results  for  the  surface 
magnetic  field  in  the  two-component  dyadic  form  (5.1),  as  found  bv 
Lee  et  al . [ 2] . 

Next,  we  try  to  express  our  results  for  the  surface  magnetic  field 
in  a four-component  dyadic  form,  viz.. 


H(Q)  = M • [b’brt  . + t’tH  _ + b’cH,  _ + t’bH  . ) / 

JO  1 1 DC  tD  / 


/ 


/ 


(5.10) 


which  is  similar  to  (5.2).  Then  the  magnetic  field  components  for  the 
cases  of  a circumferential  or  an  axial  dipole  become 
C c o -> 

H (Q)  = H,  , sin  <3  + H cos  9 + (H  + H ) sin  9 cos  9 , 

< p DD  tt  bt  tb 

HC(Q)  = - (H,  , - H ) sin  0 cos  9 + H,  sin-  0 - H , cos-  9 , 

I z bb  tt  bt  tb 

HU(Q)  = - (H,  - H ) sin  9 cos  9 - H,  cos-  9 + H sin-  9 , 

i>  bb  tt  bt  tb 

a 2 2 

H (Q>  = H,,  cos-  9 4-  H sin-  9 - (H,  + H . ) sin  9 cos  9 . (5.11) 

V*  2 x bb  tt  bt  tb 

Since  H^(Q)  = H°(Q)  as  found  in  Sec.  4.1,  we  have  H,  = H , . Then  the 
p z ’ b t tb 

system  of  equations  (5.11)  can  be  readily  solved,  yielding 

( Hlu(Q)  = H^(Q)  sin2  9 + Ha(Q)  cos2  9 - 2H°(Q)  sin  9 cos  9 , 
bb  p z z 

( H (Q)  = KC(0)  cos-  9 + Ha(Q)  sin-  0 + 2H^(Q)  sin  9 cos  9 , 

1 1 p 2 2 

. H,  (Q)  * H , (Q)  = [H^(Q)  - Ha(Q)]  sin  9 cos  9 - H^(Q)(cos-  9 - sin-  9). 
w bt  tb  P z 2 

(5.12) 


We  now  substitute  the  actual  values  of  H^(Q)  , h'"(Q)  and  H (Q) , taken 

V 2 2 


from  (3.11),  (3.26)  and  (4.6).  Then  we  obtain 

1 , 


. 2 -jks 

H (Q)  =»  (1  - -L 

bb^;  2 Tj  ks  1 ks 


2 2 ‘ 

k-s- 


1.  — / O',  1 1 / 

—Min;-  (ks)  sin-  9 - (ks) 

16ka  0 ks  1 


( ■? \ ? (2) 

2H.,  (ks^  cos  9 + ksH^-  (ks)  cos" 


(5.13' 


and 


HttM) 


k-Y  1 i 2 k2Y  (2}  7 1 (i) 

Itt  + -=rrl  - rCT^ri  <k*>  cos-  e - — h|-  (ks)]. 


2~\  ks  ks  ,22'  16ka  0 

J K S 


(5.1-*) 


and 


H (Q)  « H„b(Q)  - - -j— ~ sin  - c°s  9[3Hq-^  (ks)  + H^-\ks)] 


(5.  IS"' 


f 


Notice  that  the  leading  terras  of  (5.13)  and  (5.14)  agree  with  the  planar 
solution  for  and  H , as  given  in  [1,  Appendix  D]  and  [2,  Eqs.  2.8a,b)] 
For  large  ks  the  Hankel  functions  in  (5.13)  - (5.15)  can  be  replaced  by 
their  large-argument  asymptotic  expansions,  thus  leading  to 


k2Y  e"jks 

Hbb^>  : 


. J_  1,1  .tt.1/2  -j-r/4.,  ,1/2/,  2 43, 

1 ' ks  ' 7T7  + 4ka(7  e (ks)  \3  ' cos  e(jKS  + T} 

k s k. 


+ ks(“  8 + 12§  C°S"  9)  + °N4 

Ik  s i 


(5.16) 


k2V  e~3^s  2i  9 1 tI/2-W/4 

Htt(Q)  = 2?j  ~ ki  + T7  + 4ka(3)  "e  (ks)  ~\3  cos' 
(_  k s V. 


+ k^-1  + S co 


cos2  0)  + 0 -— jjl 


(5.17) 


,,,  ..  k2Y  e_jks  1 ,77,1/2  — i tt/ 4 ,,  ,1/2  , . 

Hbt(Q)  = Htb(Q)  - ITT  “TT  4ka  7}  e (ks)  cos  0 sin 


•{2+-+ote!} 


(5.18) 


Alternatively,  the  surface  magnetic  field  H(Q)  can  also  be  represented 
by  the  form  (5.2).  From  a comparison  of  (5.2)  and  (5.10),  we  find  that 
A,  B,  and  C are  given  by 


A = H - H - -f-  , B = H - H.  ~ , C = ■— "7-  ' . (5.19) 

bb  bt  cos  0 tt  bt  sm  9 sin  9 cos  9 
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MATCHING  TO  AN  EXPANSION  IN  TERMS  OF  FOCK  FUNCTIONS 


6. 

In  Sections  3 and  4 we  derived  approximations  to  the  surface  field 
components  and  H'j1 , due  to  a circumferential  or  an  axial  magnetic 

dipole  on  the  surface  of  a cylinder;  see  (3.11),  (3.13),  (3.26),  (3.27), 
(.4.6)  and  (4.7)  for  the  final  results.  These  approximat ions  were  obtained 
by  starting  from  the  exact  solution  in  [1]  and  replacing  the  quotient  of 
the  Hankel  function  and  its  derivative  bv  the  Debye  approximation  as 
given  in  (2.6)  and  (2.12).  Although  the  use  of  the  Debye  approximation 
has  not  been  fully  justified,  it  is  believed  that  our  approximate  results 
for  the  surface  field  components  are  valid  for  large  ka  and  small  j where 


,-1/3 


ks 


(ka) 


2/3 


cos 


(/3 


For  large  i.e.,  in  the  deep  s’nadowi  the  surface  field  decays  exponentially 
as  a function  of  Z;  see  [11  and  [2].  This  ^-dependence  of  the  surface 
field  is  properly  described  in  terms  of  Fock  functions  u ( ) and  v(£). 
Following  (2,  Appendix],  the  Fock  functions  u(q)  and  v(£)  are  defined  by 


and 


v(0 


no-;) 


ej'T/4  /2  r V, 


(t) 


, r-  ’ j W.(t) 

2v  t 4 2 


e"^C  dt 


. e3jTT/4  -3/2  f -jqt 


Vl 


J 

r 


dt 


(6.1) 

(6.2) 


where  >--(t)  is  an  Airy  function,  viz., 

u,(t)  3 /r[Bi(t)  - j Ai ( t ) ] = 2*"r  e 1 ' ^ Ai(te  ' 3) 
and  the  contour  T is  sketched  in  Fig.  3 (see  Page  23). 


(6.3) 


The  "hard"  Fock  function  v(i)  arises  when  approximating  integrals  with 

(2)  (2)  ' 

an  integrand  containing  (k^aWH  ~ (k  a) ; similarly,  the  "soft"  Fock 


( > ) ' ( 2 ) 

rune  don  u(.%)  arises  when  the  integrand  contains  H " (k^aJ/H  ~ (k^a). 

In  (L,  Appendix  B J the  notations  v^(£)  and  u I"')  are  used  for  v(q)  and 

u(.\),  respectvely.  By  closing  the  contour  T at  infinity  in  (o.l)  and  (6.2), 

we  arrive  at  the  following  residue-series  representations  for  v(.\)  and  u('): 

-ii/i  _ i / ■>  * exp  Mac’  ] 

v(0  » e vTTt,  l p 

n*l  n 


...  , J"/4  ■-  .3/2  r. 

k-j  ■ :e  ri  ; ) e> 


exp  [ - jy  t 


see  [1,  Eqs.  (B20)  and  (B37)|  and  [2,  Eqs.  (A-7)  and  (A-8)].  Here 


( o . -t  ) 


n n 


i ' 1 ' ^ and  t 1 =■  1 t ' ' e ^ \ with  Ail,- 1 1 ) » 0 and  Ai  ' (-  t ' ! ) = 0; 

n n 1 n n 


see  [2,  p.  34]  for  a table  of  the  zeros | t land  ]t'!.  It  is  clear  from 
' n n 

ib.-O,  that  the  Fock  functions  v(£)  and  u(C)  decay  exponentially  as 

t,  - •».  Notice  that  u(Q  decavs  faster  than  v ( O , since  t'1  < t . F 

n n 

small  , the  Fock  functions  can  be  represented  by  the  power-series 
expans  ions 


, ..  , /rr  it/4,3/2  7j  .3  7 ► "r  -jz/4  .9/2 

v(-.)  - 1 - - ♦ *nSm  • * 

u(-)  - 1 - 4 ♦ ■fit’**?  .-t’'-  5,;;  ♦ o(f6) 

_ 1 _ OS 


0(T) 


1,0 . 5) 


quoted  from  [1,  Eqs.  (319’>  and  (.336)  ] , and  [2,  Eqs.  i, A — 1 2 ) and  (A-13)  ] . 

We  shall  now  match  the  previous  approximations  for  the  surface  field 
components  = Ha  and  Ha,  to  a new  set  of  approximations  in  terms  of 

Fock  functions  v(;)  and  u(^).  More  specifically,  we  construct  new 
approximations  involving  Fock  functions  in  such  a manner  that  for  small 
the  new  approximat ions  reduce  to  those  obtained  in  Sections  3 and  wherehv 
the  Fock  functions  are  replaced  bv  their  power-series  expansions  in  i.n.5). 

It  is  hoped  that  these  new  approximations  are  valid  uniformlv  in  £. 
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r>.l  Magnetic  tleld  component  H, 


!'he  exact  solution  for  Ht  is  given  bv  (.3.2)  and  consists  of  two  terms 

Following  [l,  Eos.  03)  and  (19)],  we  denote  these  terms  bv  H 1 1 ' and  !i,c. 

P J) 

both  terms  are  double  integrals;  the  Integrand  of  the  term  H ’ ' contains 

l ' ) l ^ ) * C* 

the  quotient  1P“  (k^a',11  “ (k^a),  whereas  the  terra  H ' has  an  integrand 

which  contains  the  quotient  H*“'(k  a)  11*'“' (k  a).  Thus,  H’1'  should  be 

v t v t P 

matched  to  an  approximat ton  involving  the  hard  Fock  function  v("),  whereas 

il1'  is  matched  to  an  approximat ion  which  involves  the  soft  Fock  function 
*) 

u i r. ) . In  the  exact  solutions  for  H ’ 1 ^ and  H,v  , we  replace  the  Hankel 

P 

function  quotients  bv  the  Debye  approximations  in  (,2.b)  and  (.2.12), 

respectiveiv.  Furthermore  we  set  v - k^a  in  the  inner  integral.  Then  we 

are  led  to  the  following  approximations  for  and  H'  (.Q)  : 

P ' P 


H”  iQ) 


exp  [-jk  2 - jk  a<J>] 
j , 2 V 


exp  [-jk^z  - jk  a$] 


+ J a*)dk  dk 
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» k ~ k - k 


- , - V 


( r> . o 1 


It’  (O'  ' - + 


~ 1 k j l a 


Si  .cu.  ■ ' 
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Vk“  I 


xn  [ - ] k _ z - jk  a$] 


exp  1 - ) k z - jk  a*J>] 


i>kt~kv  1 _J_ 

2 2a  2 2 

k k -k 

t tv 


Jk  dk 

V 2 


* it  -k  -k  . , 

. .J l 

> ^ > . ■* 


i t >dk 


Notice  that  the  sum  of  (6.6)  and  (.0.7)  Is  equal  to  (3.3). 

The  double  Integrals  In  (o.o)  and  (o.7)  can  he  expressed  In  terms  of 
l^laf.c)  ind  l, (a-?,;),  as  defined  bv  (3.-*)  and  (3.5),  and  another  kev 
Integral  I, tad,;)  defined  by 


Ijta^.c) 


j i 


exp  [-Jk#z  - jk  avfl 


k~  - k 


t dk  dk 


lb. 8) 


I'he  latter  integral  can  he  evaluated  in  closed  form  by  use  of  the  relation 


r > rr 

l ; + k 'I,(a?.!'  - , exp  [-.jk  .. 

)z~  1 • 


- jk  a$|»k”  “ k~  - k~  dk  dk 

V V Z V 


" — - 1 + —7  + k"  I (a«?,z)  . (6.1)) 

JU4>)'  )2"  J 

On  substitution  of  the  value  of  I la^.z)  fram  (3.t>),  we  are  led  to  the 
following  differential  equation  for  I . lag>, z) : 


k")  I laJJ 


z) 


k~)2n  1 


-jss  jk 


— 3/2 

2t  k 


3/2  (2) 

H,,;  Iks) 


to. 10) 


/ j 3 

where  s ■ > tar">)  + z is  the  distance  from  the  source  point  P to  the 
observation  point  Q along  the  surface  ray;  see  Fig.  1. 

The  differential  equation  (r>.lf>)  can  be  solved  bv  Fourier  transformation. 
On  defining  the  Fourier  transform  of  bv 


il  j(.O.z)1 


l j ini'.  z)e'  " 1 d; 


aOO 

l^taO.z)  cos  tzt)  dz 
0 ( o . 11 ) 


we  readllv  find  from  (o.ld) 


.'0 


F„ll3(a$,e)} 


T~ — ” • [ (a<f ) " + z"  ) H “ + z ) < 


' j Uv' 


_L  K UC/t"  - k") 


(6.12) 


where  the  Fourier  transtorn  of  the  Hankel  function  was  quoted  from 


o,  Eq.  1.13(42)).  In  to. 12)  it  is  understood  that  »'t“  - k~ 


k“  - t 


when  t < k,  in  accordance  with  k having  a small  negative  imaginary  part. 
Bv  Inverse  Fourier  transformation  we  have  from  to.  12), 

n - 

’rT,  l f”  K.(atp*t“  - k“>  , 

I.(a£,z)  « - -^~(a>?)  2 cos  tzt)  dt  «-2irk  (a£)  ”e 

0 /t”  - k" 


2-t -jks 

k s cos 


to. 13) 


where  the  Fourier  cosine  transform  of  was  obtained  t rom  [o,  Eq . I.13ls4)] 
"he  present  result  for  l (a>7,z)  was  checked  bv  back-substitution  into  tb.10) 
In  fact,  the  general  solution  of  the  differential  equation  (0.10)  is 
given  by 


[ ,(al>, z) 
3 


- 1 - ’ - i ks 

-2Tk  (a2>)  e 1 ■+■  A cos  kc  +■  3 si 


(6.14) 


where  A and  3 are  arbitrary  constants  which  may  depend  on  ao. 

Now  since  [ (a->,z)  is  an  even  function  of  one  has  3 •*  0.  By  a direct 
calculation  of  I (a^'.c  3 0)  it  is  found  that  A 3 0 as  well.  Thus  the  result 
in  (b.13)  is  cor rec  t . 

Ke  now  return  to  the  magnetic  field  components  H''1"  and  H'1,  as 

I 

given  bv  (o.o)  and  to. 7).  As  mentioned  before,  the  double  integrals  car. 


There  are  some  obvious  misprints  in  this  formula:  in  the  first  part  ot  the 

i 2)  i ' i 

transform  result  H should  be  H . : the  range  or  validitv  of  the  sec  >nd 
v v i l 

‘ > 

part  should  be  b < y<  v instead  of  "0  < v^b.  All  results  have  been  checke 
in  an  independent  manner  bv  means  of  [3,  sec.  13. s']. 


;e  expressed  in  terras  of  derivatives  of  the  key  integrals  1^.  I,  and  1^. 


Thus  we  find  for  H''  , 


H ' ' ° (Q)  Vlik"  + I (a.p,z)  - k:i  (ao,2)  +^7 -7 7 I7(a$,z)] 

4"”k  3(a,p)~  ~ ?'(ai)“3z' 


(6.15) 


and  for  H'\ 


H’C(Q)  = - (1,(31),  z)  - -f-'  k ~ + y + I (a*,z)].  (6.16) 

4r“k  ~a  3 (dll)"  3z~ 

'.»'e  now  insert  the  actual  values  of  I ^ , I,  and  I?  as  given  by  (3.6),  (3.9) 

and  (6,13).  The  derivatives  of  the  Hankel  function  H|~\ks)  are  evaluated 

and  simplified  by  means  of  the  well-known  recurrence  relations  for  Bessel 

functions;  see  e.g.,[3.  Sec.  3.2]  and  [4,  Eqs.  9.1.27  and  9.1.30].  Thus 

c c 

we  obtain  as  our  final  result  for  the  constituents  H''  and  H'  : 

} $ 

I • 1 , . 

t.  ”y  ->  J ^ 1 "*  i 1 •■) 

H'^CQ)  = ~ ~ [sin-  d + ri-(2  - 3 sin"  9 - ■— -7  -)  + -7-7(2  - 3 sin"  0)] 

2 - *T J K S n c ~ 4 1.  ~ 0 ~ 

v_  O b K.  ^ 


s)  - — H,  (ks)  + ksH,  7(ks)  cos“  £ sin* 

10K..1  u pCS  1 J 


h;c(Q)  ~ ~ (tfr  — k—'\  - r-f-  [2H^2>(ks)]  . 

yb  -^1  Kb  k*  - x 0 


(6.17) 


(6.  IS) 


As  before,  we  observe  that  except  for  the  Debve-type  approximation  to  the 
Hankel  function  quotients,  no  further  approximat ions  were  involved  in  the 
derivation  of  16.17)  and  (6.18).  The  sum  of  (6.17)  and  (o.iS)  is  equal  to 
the  field  component  H'" < Q)  as  given  by  (3.11). 

f 

For  large  ks  the  Hankel  functions  in  (6.17)  and  (6.18)  can  be  replaced 
bv  their  Large-argument  asvmntotic  expansions,  thus  leading  to 


H'  (Q)  = ~*f  — ; sin-  9 + r^-(2  - 3 sin- 

_ ” 1 rK  S N ^ 
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and 
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(6.20) 

We  now  natch  these  approximations  to  a new  set  of  approximat ions  of 
a form  similar  to  that  of  Lee  et  al.[2]  (see  (3.17)),  involving  Fock 
functions.  The  new  approximation  for  H’’1-  involves  the  hard  Fock  function 
vie)  and  its  derivative  v'(-),  whereas  the  approximation  for  H ' contains 
the  soft  Fock  functions  u(%)  and  u'(%).  Therefore,  guided  by  i3.17),  we  set 


.--v  ,~Jk3  ,-1/3,  /, 

H',CiQ)  = 777  rAv(£)  + -=  w4  COS  7 j • Bv’(*) 

**  - T 1 I..S  , , - , J 
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(6.21) 


and 


,-1/3, 
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’'‘(Q)  ~ 4rt  l^uic)  + — ^/4  cos'4  9 Du’ (5)1  . 

ika)-/3 


io 


2) 


where  the  constants  A,  3,  C and  D are  determined  hv  matching  to  (0.19) 
and  (6.20)  for  small  5* 

For  small  % we  replace  v(5)  and  u(5)  and  their  derivatives  by  the 
approximat ions 


...  , j tt/4  _ 3/ 2 

v(^,)  = l — r 


..  , vT  j T / 4 _ 3 / 2 

n ( ., ) * 1 r * 
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v ' ^ ) 3 e 


3k_t  Jtt/4.1/2  , _ 3/ir  j^/4_l/2 
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(6.23) 


taken  from  (6.5).  Then,  the  results  (6.21)  and  (6.22)  become 
2 - iks 
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, ,c  _ k \ e 
ri.  vQ)  = -7 


1 i 1/9  -iit/i  1/2  ">  3 2 

[A  + 7j— -(— ) ' “e  J (ks)  {- j ksA  cos-  9 + — B cos  9}], 
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and 


? -»  i ICS 

H;C(0)  - tC  + T^-(^)L/2e~r/4(ks)1/2{-2jksC  cosi 2 *  9 + 3D  cos-  9}] 

(6.25) 

5v  identifying  (6.24)  and  (6.25)  with  (6.19)  and  (6.20),  respectively, 
we  can  readily  determine  the  constants  A,  B,  C and  D,  viz. , 


A = sin-  9 , (2  - 3sin‘ 
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1 2 
— — ) + --~(2  - 3 sin  o) 


n _ i 11IL 

3 ' 3 ■ 

cos 
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sm 


cos  9 k s 
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1 .3  7 sin- 


0 + ks(4  " 12 


cos 


+ sin-  9)  , (6.26) 


and 


C = 7 


ks 


cos 


12ks 


(6.27) 


cos 


The  values  of  A,  B,  C and  D,  thus  obtained  are  to  be  inserted  into  (6.21) 
and  (6.22).  Then  by  addition  of  (6.21)  and  (6.22),  we  obtain  the  following 
approximation  in  terms  of  Fock  functions  for  the  surface  field  component 
H*(Q): 


2 — j r 

H>>  = M si*2  9 v(^  + i(2  - 3 Si"2 
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I~)v(D  + ks  — u(^) 

cos  v cos  8 
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(6.2S) 
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"he  present  approximation  is  of  a similar  form  as  that  of  Lee  et  al.[2], 


as  given  In  (3.1'3.  Notice,  however,  that  the  coefficients  of  v(t),  u(.y), 
v ' (~.)  and  a ' (,  \ 3 are  somewhat  different. 
r>.2  Magnetic  field  component  » H . 

As  pointed  out  in  Section  4.1,  the  exact  solution  for  the  field 
component  due  to  a circumferential  dipole  is  identical  to  that  for 
the  field  component  due  to  an  axial  dipole.  Then  the  approximate 
solutions  for  and  are  also  the  same. 

The  exact  solution  for  is  given  in  (3.22),  as  quoted  from  [1]. 

The  integrand  of  the  double  integral  in  (3.223  contains  the  quotient 
(2)  , (2)’ 

H tk^a)  H lk  a) . Therefore,  we  match  our  previous  approximation  to 
a new  approximation  which  involves  the  hard  Fock  function  v(~)  and  its 
derivative  v'<,~3.  Thus  we  start  from  an  approximation  for  H (Q)  of  the 
form  (.6.213  with  constants  A and  3 yet  to  be  determined.  For  small  q , 
when  replacing  v<,^)  and  v ' (' ) bv  v6.233,  this  approximation  passes  into 
the  form  (6.243.  Then  by  identifying  this  form  with  our  previous 
approximation  (3.2~3  and  (3.28),  we  find  for  the  constants  A and  3: 


A ” - cos  1 sin  ~[l  - f 


Jj  3 


d » - cos  r stn 


4i_v:  . i + J-(1E 

- , a ks  Ou 


4'ion  inserting  these  values  into  the  form  (.6.213,  we  obtain  the  following 
approximation  in  terms  of  Fock  functions  for  the  surface  field  components 
H*”  (Q3  - H*(Q>: 
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ri  (Q)  = rl.(Q)  = - 777  — : cos 
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,-1/3, 


-n,-;'  - r)v’(V  + V“)v’(^:- 

cos”  3 143  °"t  i cos”  J 


(6. JO) 

The  present  approximation  should  be  compared  to  the  result  of  Lee  et  ai.[2] 

as  given  in  (.3-31)  . The  main  difference  is  that  the  result  in  (3.31) 

involves  both  the  hard  and  soft  Fock  functions,  whereas  our  approximation 

in  1,6.30)  is  in  terms  of  the  hard  Fock  function  only. 

6.3  Magnetic  field  component  Ha. 

z 

The  exact  solution  for  the  field  component  Ra  due  to  an  axial  dipole 

z 

is  given  bv  (-.,2),  as  quoted  from  flj.  Since  the  integrand  of  the  double 

( ’ ) (■•>>' 

integral  in  (4. 2)  contains  the  quotient  H (,kra)  H, ' (k^a)  , we  construct 

an  aoproximat ion  for  H 3 ('' ) in  terms  of  the  hard  Fock  function  v(^)  onlv. 

We  start  again  from  an  approximation  of  the  form  (6.21)  with  constants  A 
and  3 to  be  determined.  For  small  1 we  replace  v(^)  and  v’(2)  by  (6.23), 
then  th  approximation  for  Ha(Q)  reduces  to  the  form  (6.24).  Then  by 
matching  this  form  to  our  previous  approximation  (4.7)  and  is. 3),  we  car. 
readily  determine  the  constants  A and  5,  viz., 

A = cos”  J + A(2  - 3 cos”  2)  + -7— (2  - 3 cos”  6)  , 

K3  k“s” 


11  2 , . 1 , 7 2 .1ST  2 . , 

= - — cos  3 + r^(-  ^ + — cos  :-). 

3 cos 


(6.3D 


On  inserting  tnese  values  into  the  form  (6.21),  we  obtain  the  following 


approximation  in  terms  of  Fock  functions  for  the  surface  field  component 


Ha(Q): 
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i (Q) 


cos  •?  v(^)  + r*y(2  - 3 cos  9)v(£)  + - , ( 2 - 3 cos  ■?)'» 


- 3 , • II  - - .,  . 1 , 11  2 sin"  9 
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(6.32) 


The  present  ipproximation  should  be  compared  to  the  result  of  Lee  et  al.[2] 
is  given  In  (-*.11).  Again  the  main  difference  is  that  the  result  in  (4.11) 
involves  both  the  hard  and  soft  Fock  functions,  whereas  our  approximation 
in  ^.32'  contains  oniv  che  hard  Fock  function  v(”)  and  its  derivative. 
i . -*  Representat  ion  in  a dvadic  form . 

As  in  Section  3,  we  now  represent  che  surface  magnetic  field  at  Q 
iue  to  a magnetic  dipole  M at  ? in  the  dyadic  form  (5.10).  vie., 


K(Q)  - M • [b'hHbb  + t’tH  + b’cH^  + t’bK  fe] 


(6. 331 


^here  the  unit  vectors  b',  b,  t'  and  t are  shown  in  Fig.  2.  Then  the 

ivadic  components  H,  , , H , H,  and  H , are  related  to  the  surface  field 
r bb  tt  bt  tb 

C C 

components  H * H and  H'  through  (3.121.  In  (3.12),  replace  the  field 

components  by  their  approximations  (6.231.  (o. 301  and  (6.32)  in  terms  of 

Fock  functions.  Then  we  obtain  the  following  approximations  in  terms  of 

Fock  functions  for  the  dvadic  components  H,  . . H md  H,  = H . : 

* bb  tt  bt  tb 
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bo  2j  ks 
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c cos  - k s 
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7.  NUMERICAL  RESULTS 


For  the  numerical  calcu lations , we  concentrate  on  H (the  ^-component 

v 

of  H due  to  a circumferential  dipole)  and  (the  z-component  of  H due  to 

an  axial  dipole).  Three  sets  of  approximate  solutions  are  used  in  the 
calculations,  namely, 

(i)  the  "full  formulas"  of  Chang  et  al.  [1],  which  are  given  in 
(3.16)  and  (4.10)  of  the  present  report; 

(ii)  the  formulas  of  Lee  et  al.  [2],  which  are  given  in  (3.17)  and 
(4.11);  and 

(iii)  the  present  formulas  in  (6.28)  and  (6.32). 

The  radius  of  the  cylinder  is  ka  = 9.5325.  The  ray  directions  are 

0=0°,  45°  and  90°.  As  a function  of  ks , numerical  values  of  HL  are 

•? 

presented  in  Figures  4 to  6,  and  those  of  in  Figures  7 to  9.  In 
these  figures,  we  use  the  following  notations: 


Mag  = absolute  value  of  (H  /k*”)  or  (H  /k*’) 

<P  2 

, . , . , , ,„c  -Hks.  ,,a  + iks. 

Normalized  phase  = phase  of  (H^e  ) or  ) 


The  db  value  of  Mag  is  calculated  by  20  log^IMag  in  ampere-meter)  . 

As  another  accuracy  test,  we  calculate  the  mutual  admittance  Y 0 
between  two  identical  slots  on  a cylinder  by  using  the  three  formulas 
described  in  (i)  to  (iii)  above.  The  geometry  is  sketched  in  Figure  10 
with  the  parameters 

Frequencv  = 9 GHz  , X = 1.3123"  , a = 1.991" 

Slot  dimension  = 0.9"  x 0.4" 

Slot  separation  described  by  ($  ,z^) 
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NORMALIZED  PHASE  (deg) 


Lee [2] 


NORMALIZED  PHASE  (deg) 


-80 


MAG  (dB) 


NORMALIZED  PHASE  (deg) 


-92 


Figure  4a  except  fur  0 = 90°  and  larger  ks. 


NORMALIZED  PHASE  (deg) 


NORMALIZED  PHASE  (deg) 


The  steps  for  calculating  Y , can  be  found  in  [i],  [ d ] and  ['].  The 
results  for  two  i identical)  circumferential  slots  are  given  in  Table  1, 


and  those  for  the  two  axial  slots  in  Table  2.  Values  of  Y,  are  listed 

in  (db  * JO  log  Y,  , , ohase  in  degree'  format.  In  addition  to  three 
iu  i- 

asymptotic  solutions  based  on  *.i'  to  (iii>.  an  exact  solution  of  Y , 
calculated  by  the  modal  series  [7]  and  [S]  is  also  listed  in  the  two 
tables  for  comparison  puroose. 


- -4. 


8.  CONCLUSIONS 


The  surface  magnetic  field  due  to  a magnetic  dipole  on  a cylinder 
can  be  found  exactly  in  terms  of  cylindrical  modal  functions  and  Fourier 
integrals.  This  solution,  however,  is  not  suitable  for  computations  at 
high  frequencies  (ka  >>  1)  because  of  its  slow  convergence  rate.  The 
present  paper  is  devoted  to  extracting  an  asymptotic  solution  (ka  -*•  oo) 
from  the  exact  one.  Explicit  results  have  been  obtained  for  the  following 
cases : 

(i)  In  the  penumbra  region  on  the  cylinder  where  £ is  small, 

the  asymptotic  dyadic  Green's  function  is  given  in  (5.10)  and 
(5.13)  - (5.15). 

(ii)  In  order  to  obtain  a solution  uniformly  valid  for  all  points 
on  the  cylinder  (from  the  penumbra  to  the  deep  shadow),  the 
asymptotic  solution  in  (i)  is  matched  to  the  well-known 
creeping  wave  representation  via  Fock's  functions.  The  final 
dyadic  Green's  function  is  given  in  (6.34)  - (6.36). 

The  present  solution  has  been  compared  with  two  previous  ones: 

Chang  et  al.  [ 1 ] and  Lee  et  al.  [2].  Of  particular  interest  is  that, 

through  rigorous  asymptotic  expansions,  we  have  confirmed  the  peculiar 
-1/2 

(ks)  field  behavior  along  the  generator  of  the  cylinder  (sometimes 

known  as  the  "transverse  curvature  term"  in  the  field  solution).  This 
term  plays  a most  important  role  for  rays  parallel  or  almost  parallel 
to  a generator  of  the  cylinder.  Until  the  present  confirmation,  its 
existence  was  predicted  only  through  speculation. 
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1.  INTRODUCTION 


In  the  design  of  a conformal  array,  the  two  most  important 
electromagnetic  parameters  are  the  mutual  admittance  between  elements 
and  the  active  element  pattern.  These  two  parameters  have  been  calculated 
by  the  following  two  techniques: 

(i)  Modal  analysis  [1]  - [3]  applied  to  problems  with  separable 
geometry.  The  solution  is  usually  rigorous,  and  is  in  the  form  of 
infinite  series/ integrals . Because  of  the  convergence  rate,  it  is  suitable 
for  numerical  calculations  only  when  the  radii  of  curvature  of  the  array 
surface  are  small  in  terms  of  wavelength.  In  other  words,  modal  analysis 
is  a low-frequency  technique. 

(ii)  Ray  technique  [4]  - [10]  is  based  on  surface  rays,  first  introduced 
by  Keller  in  his  Geometrical  Theory  of  Diffraction  (GTD).  It  normally 
yields  an  asymptotic  solution  valid  for  high  frequencies.  Because  of  the 
wide  range  of  its  applicability  and  simplicity  of  its  final  solution,  the 
ray  technique  is  a most  attractive  tool  in  solving  conformal  array  problems. 

This  part  of  the  book  will  describe  the  ray  technique  for  calculating  the 
mutual  admittance  and  active  element  pattern.  We  will  concentrate  on 
conformal  arrays  which  have  rectangular  slots  as  their  elements. 


CIRCUIT  DESCRIPTION  OF  SLOT  ARRAY 


">cr^ 


Consider  an  array  of  N slots  over  a curved  conducting  surface  I 
(Figure  I).  Each  slot  is  fed  by  a rectangular  waveguide  (Figure  2), 
where  only  the  dominant  TE^q  mode  propagates  and  all  other  modes  attenuate. 
The  electromagnetic  properties  of  the  array  can  be  conveniently  described 
by  circuit  parameters  detailed  below. 

Let  us  concentrate  on  a typical  element  n in  the  array.  At  a suffi- 
ciently large  distance  i from  the  aperture,  only  the  dominant  TE^  mode 

t h 

is  present.  Then  the  transverse  field  vectors  in  the  n guide  'an  be 
represented  by 

E(x,v,z=-£)  * V e (x,  y)  (2.1a) 

n 

i 

H(x,y ,z»-i)  =»  I (z  x e)  (2.1b) 

n 

where 


^7  1/7 

e(x,y)  = v(— ) “ cos 

■ ab 


V = modal  voltage  in  n 
n 


I = modal  current  in  n 
n 


■ IT 

(- 


x) 
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th 

th 


element 

element 


12. 2a) 
(2.2b) 
(2.2c) 


Note  that  the  field  in  (2.1)  is  the  total  field  consisting  of  waves 
travelling  in  both  +z  and  -z  directions.  Because  of  the  linearity  of 
the  Maxwell's  equations,  the  current  in  the  mt!l  element  is  linearly 
proportional  to  the  voltages  in  all  elements  in  the  array,  i.e.. 


[ 

m 


N 

V Y V 

mn  n 


m 


1,2, 


,N 


(2.3) 


In  matrix  notation,  (2.  3)  may  be  rewritten  as 

I - YV 


(2, 


where  t and  V are  column  matrices  with  elements  tl  } and  (V  },  and  Y is 

n n 

a square  matrix  with  elements  Y ). 

mn 

The  proportional  constant  Y.,  in  (2.3),  tor  example,  is  called  the 
mutual  admittance  between  slots  1 and  2.  By  reciprocity,  Y 0 ■ Y0  . 

'v’e  may  calculate  (measure)  Y , from  the  following  setup  (Figure  3): 

(1)  Element  l is  excited  so  that  the  (total)  voltage  at  the 
reference  plane  a » -i)  is  V . 

(ill  Conducting  planes  are  placed  at  the  reference  planes  of  all 

other  elements  so  that  V'  ■ 0 for  n * L . 

n 


Then  it  follows  from  (2.3>  that 


12  V. 

1 I o 


(2. 


short  all  except  l 


which  may  be  considered  as  the  definition  of  Y . 


th 


As  a transmitting  antenna,  the  n element  in  the  slot  array  in 


Figure  l is  excited  bv  an  incident  TE.  . mode  with  voltage  V , where  the 

10  n 


superscript  signifies  that  the  wave  propagates  toward  the  aperture 
in  the  direction.  The  discontinuity  at  the  aperture  causes  a reflec 


TE,  , mode  with  voltage  V .which  travels  in  the  -z  direction.  Then  the 
10  n 


(totals  voltage  at  che  reference  plane  i.c  = -(.)  is 


V - Y+  + Y~  (2. 

n n n 


while  its  correspond lng  current  is 


where  V is  the  characteristic  admittance  of  TE  mode 

0 IU 

Y • — - — [ l - 

0 1201T  11  vka 


(2.7) 


■f  , 


For  a tiven  set  of  incident  voltages  \V  , we  can  determine  the  reflected 


voltages  (V  ' and  the  (.total)  voltages  (V  ' from  (J.o)  and  (l.~>).  The 
n n 


results  are 


» *—  ] a — ] * ■—  1 m — -f. 

V - (1  + Y Y)  U - Y Y ) V 


U’.  8) 


— a a — 1 * — 1 — + 

V - :(1  + Y Y)  V 


i.:.h) 


where  1 Is  an  Identity  matrix,  and  Y * Y 1, 

c c 

In  addition  to  admittance  tY  },  another  set  of  important  parameters 

mn 

is  the  short-circuited  active  element  patterns  IP  0,1')}.  For  example, 

n 

P 10,>>)  is  defined  as  the  radiation  pattern  (o  t E.  or  E component"'  when 

I O 0 

V « 1,  and  V * 0 if  n f l (.Figure  3).  The  term  "active"  means  that  the 
l n 

radiating  element  is  situated  in  an  array  environment  (.not  in  an  isolated 
environment).  When  an  array  is  excited  by  an  incident  voltage  vector  V , 
the  radiation  pattern  of  the  whole  array  is  then  given  by 


P (f,t>)»  ' VPi\4>)"VP 

arrav  n n 

n- 1 


'.10) 


where  \V  • can  he  calculated  from  ld.l>),  and  T is  a transnose  operator. 

ti 
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5.  RCA  ITER  l Nil  DESCRIPTION  OF  SLOT  ARRAY 


For  the  same  slot  arrav  in  Figure  1,  a different  and  equivalent 
description  may  be  given  in  terms  of  scattering  parameters,  instead  of 
circuit  parameters. 

Parallel  to  (2.3),  the  basic  relation  in  the  second  description  is 


N 

V S V 
“ , " mn  n 


m ■ 1,2, 


(3. la) 


>r  in  matrix  notation. 


V - rv 


(3. lb) 


Here  R * [S  ] Is  a scattering  matrix.  R,  , - R,,,  for  example,  is  the 
mn  l - - 1 

induced  voltage  at  element  2 when 

tl)  element  1 is  excited  with  - l (not  V « 1),  and 
(il)  all  other  elements  are  terminated  with  a matched  load, 
in  the  manner  sketched  in  Figure  . Rometimes,  R is  also  known  as  the 
coupling  coefficient  between  elements  1 anil  2.  I'he  comparison  of  (3.1b) 
with  (2.R)  leads  Immediately  to 


* m 1*  —1  m ■ — 

R - (1  + Y Y)  (I  - Y Y) 


( 3.2) 


which  relates  R to  the  admittance  matrix  Y.  For  the  special  case  N » 2, 


we  have 


"Vu 

* hi’'  - vi; 


(3.3) 


For  a given  Incident  voltage  vector  V , the  (voltage)  reflection 


coefficient  in  element  m is  defined  bv 


_rn 

+ 


m - 1.2, ...,M 


(3.4a) 


and  is  found  from  (,3.1a)  to  be 


v S (V  /V  ) 
, mn  n m 
n«l 


( 3.  -»b) 


The  Input  admittance  of  the  TE^  mode  in  element  m is  given  by 

v V l + R 

,,(in)  _m  m „ m_ 

m I " c 1 - R 

m m 


(3.5) 


L’nlike  Y . we  note  that  R and  are  functions  of  the  arrav 

mn  m tn 

excitations. 

Under  the  condition  sketched  in  Figure  4,  the  radiation  pattern  is 
called  the  match- loaded  active  element  pattern  For  a given  incident 

voltage  vector  V , the  pattern  of  the  whole  array  is  given  by 


P (0,4>)  - y V+Q  (e,4t)  - <V+)TQ 

irrau  u n 


arrav 


n-1 


n n 


( 3 . t> ) 


With  the  help  of  (2.4),  the  comparison  of  (3.b)  and  (2.10)  leads  to 

— * * — 1 « —1  T— 

Q - 2[ (1  + Y Y)  ] P 


1 3.  :> 


which  relates  two  types  of  active  element  patterns.  Note  that  Q^(0,'i')> 
for  example,  depends  on  P^(.0,$),  P,(0*$) , . • • , P^(0,  t)  through  the  matrix 


relation  in  (3.7). 


ONE-MODE  APPROXIMATION 


In  the  discussion  of  Sections  2 and  3,  the  reference  plane  for  the 
voltage  and  current  is  taken  to  be  a distance  l from  the  aperture 
(Figure  1).  Specifically,  l should  be  chosen  sufficiently  large  so  that 
all  reflected  modes  other  than  TE^  attentuate  to  negligible  values  within  i. 
As  an  example,  with  parameters  (Figure  2) 

a = 0.9",  b = 0.4",  f = 9 GHz 

Z should  be  at  least  0.45"  in  order  that  the  next  higher-order  mode  TE,^ 

attenuate  to  one-tenth  of  its  magnitude  within 

For  a finite  Z,  the  calculation  of  {Y  } and  other  scattering 

mn 

parameters  is  quite  difficult.  Hence,  in  practice,  we  often  set 

Z = 0 . (4.1) 


When  (4.1)  is  used,  all  of  the  analysis  in  Sections  2 and  3 becomes 
approximate.  This  approximation  is  valid  if,  despite  the  discontinuity 
of  the  guide  and  the  coupling  in  the  array,  the  aperture  field  of  the  slot 
essentially  contains  ruj  other  modes  than  TE^.  For  this  reason,  the 
approximation  in  (4.1)  is  known  as  the  "one-mode  approximation."  It 
has  been  verified  experimentally  and  theoretically  that  one-mode  approximation 
is  a good  one  if  (i)  the  slots  are  thin,  and  (ii)  their  length  is  roughly 
a half-wavelength. 

Under  the  one-mode  approximation,  the  expression  of  the  mutual 
admittance  in  (2.5)  can  be  replaced  by 


_L_  ( 

V1V2  J 


(4.2) 


where 


A,,  = aperture  of  slot  2 

= magnetic  field  when  slot  1 is  excited  with  voltage  , and 
all  other  slots  are  covered  by  perfect  conductors  at  their 
openings  (i  = 0 in  Fig.  2) 

E,  = electric  field  when  slot  2 is  excited  with  voltage  V , and 
all  other  slots  are  covered  by  perfect  conductors  at  their 
openings . 


Because  = I^h,,  and  E0  = V,e,,  it  is  a simple  matter  to  verify  that 
(•i.2)  and  (2.5)  are  equivalent  under  the  one-mode  approximation. 

There  is  an  alternative  definition  of  mutual  admittance.  Instead 
of  (2.1),  a modal  voltage  (with  a bar)  may  be  defined  through  the 
expression  for  the  aperture  field  of  slot  1 as  follows: 


E = y z V cos  (—  x) 

b 1 3 


(4 . 3a) 


or  equivalently 


V = | (v  • E)  . dv 


(4. 3b) 


Then  a different  mutual  admittance  Y^-,  is  defined  by  (4.2)  after  replacing 


by  (V^,V,).  It  can  be  shown  that 

Y = — y 
12  2b  12 


(4.4) 


Two  remarks  are  in  order:  (i)  In  the  limiting  case  that  b — 0,  Y goes 
to  zero  in  proportion  to  b,  whereas  Y approaches  a constant  independent 
of  b.  (ii)  For  the  special  case  when  a = \/2  and  the  slots  are  arranged 
on  a plane  (planar  array),  it  is  Y^,  not  Y^,,  that  is  related,  by  the 


3 


3abinet  principle,  to  the  mutual  impedance  Z ,,  between  two  corresponding 
dipoles  calculated  by  the  classical  Carter's  formula  [11],  [12].  (iii)  When 
the  slots  are  excited  by  waveguides  (.transmission  lines),  one  often  uses 
Y 'i , (A' 'i , ) . From  here  on,  we  will  concentrate  on  instead  of  Y^,. 

Under  the  one-mode  approximat ion , the  short-circuited  active  pattern 
P^(0,J)  of  slot  l (Figure  3)  becomes  the  pattern  of  a single  slot  when 
all  other  slots  are  covered  by  conductors  at  their  openings  (i  * 0 in 
Figure  3).  Its  calculation  is  thus  greatly  simplified. 
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5.  GTD  GREEN’S  FUNCTION  FOR  SURFACE  FIELD  ON  A CYLINDER 


Under  the  one-mode  approximation,  the  mutual  admittance  between  two 
slots  in  an  array  can  be  calculated  from  (.4.2).  L'e  will  now  applv  it  to 
a slot  array  on  an  infinitely  long  conducting  cylinder.  The  key  step  lies 
in  the  calculation  of  H^,  the  magnetic  field  at  the  aperture  of  slot  1 
due  to  a voltage  excitation  in  slot  1.  To  this  end,  we  consider  the 
following  Green's  function  problem. 

At  point  Q’  on  the  surface  of  the  cylinder  of  radius  R (Fig.  Sal, 
there  is  a tangential  magnetic  dipole  source  described  bv  a magnetic 
current  densitv  (for  exp  +juit  time  convention) 

Kir)  * !l  ; Mr  - R)5i4>)6i.c)  (5.1) 

where  M is  the  magnetic  dipole  moment,  and  ir  = R,  1=0,  z = 0)  are 
the  cylindrical  coordinates  of  Q’ . The  problem  is  to  determine  H at 
another  point  Q = (R,vJ,z)  on  the  same  surface.  The  rav  technique  described 
below  applies  when  kR  is  large  (sav  10  or  more). 

According  to  GTD  113],  [la],  the  dominant  contribution  of  H at  0 is 
the  field  on  the  surface  ray  from  Q’  to  Q.  The  surface  rav  is  a geodesic 
on  the  conducting  surface,  and  in  the  present  case  is  a helical  path 
(Figure  5).  The  arc  length  of  the  surface  rav  is 

•>  > 

s = * \ Rb ) ~ + z ~ . (5.2) 

The  tangent,  normal,  and  binormal  of  the  surface  rav  are  it’,  -n',-b’) 
at  O',  and  (t , -n,  -b)  at  Q.  Thus,  it.  n,  b)  form  a moving  trihedron  along 
a surface  rav,  pointing  toward  tiie  longitudinal  and  two  transverse 
directions.  At  anv  point  on  the  surface  ray,  the  curvature  of  the 

conducting  surface  is  described  bv  two  parameters: 

10 


R = the  radius  of  curvature  in  the  direction  of  t (or  that  in  the 
longitudinal  direction  of  the  survace  ray) , and 
= the  radius  of  curvature  in  the  direction  of  b (or  that  in  the 
transverse  direction  of  the  surface  ray) . 

On  a convex  surface,  both  R and  are  nonnegative.  For  the  present 
case  of  a conducting  cylinder,  one  has 

R „ _ R 


R = 
t 


2 „ 
cos  0 


• 2 a 
sin  9 


(5.?' 


where  8 is  measured  from  the  R<J)-axis  in  Fig.  5b,  and  takes  a value  between 
0 and  2tt . The  large  parameter  for  our  asymptotic  expansion  is 


, 1 , . ,1/3 
m = (—  kR  ) 


(5.4) 


Thus,  the  solution  to  be  presented  is  an  approximate 

asymptotic  solution  valid  for  m -*■  ».  Furthermore,  let  us  introduce  a 

distance  parameter 


? -f  - <k/2R-y/:i  . - -!% 

t 2m 


(5.5) 


which  is  the  arclength  normalized  by  k and  R . Mote  that  C = 0 defines  the 
lit  region  (G  = tt / 2 ) , E,  ~ 1 defines  the  penumbra  region,  and  E,  » 1 defines 
the  deep  shadow  region.  Our  solution  is  uniformly  valid  for  all  5 _>  0. 

Due  to  the  point  source  in  (5.1),  the  final  asymptotic  solution  for 
the  magnetic  field  on  the  surface  derived  in  [9]  is  given  by 

A A y\  /v 

H(Q)  * M • (b’bHh  + t’tH  ) (5.6a) 

where  the  transverse  component  is 

^(Q)  ^ [(1  - j^-)v(C)  - (j^)2  u(C)  + j(/2  kRt)"2/V(S) 


+ j(/2  kRc)'2/3(Rc/Rb)u'(^)]  G(s)  , 


(5.6b) 


I 
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the  longitudinal  component  is 

Ht(Q)  ^ (^-)[v(5)  + (1  - j^)  u(£)  + j(/2  kK^)  ~'^u'(a)1  C(s)  (5. no) 
and  the  function  G(s)  is 


C (s ) 


k~Y . -jks 
0 e J 


2k  i 


ks 


( 3 . bd ) 


1/2  -1 

= (e^/p^)  “ = (120tt)  , v and  u are  defined  in  the  Appendix,  and  v’ 

is  the  derivative  of  v.  The  solution  in  (5.6)  is  largely  based  on  the 
classical  work  of  Fock  [17]. 

Let  us  consider  several  limiting  cases  of  the  solution  given  in 
(5.6).  If  the  radius  of  the  cylinder  becomes  infinite 

kR  -*■  °°  (5.7) 

the  use  of  ( A— 11)  through  (A-15)  in  the  Appendix  in  (5.6)  leads  to 


^(Q)  ^ [l 


ks 


ks 


G (s ) 


kR  - - 


(5.3a) 


H (Q)  ^ (f^-)d  - r^-)G(s)  , kR  - * . (5.3b) 

t ks  ks 

When  (5.8)  is  substituted  into  (5.6a),  we  find  that  H in  (5.6a)  is 
identical  to  the  exAAt  solution  of  the  surface  field  due  to  a magnetic 
dipole  on  a flat  ground  plane  [15], 

The  second  limiting  case  occurs  when 

0 - it/2  . (5.0) 


We  find  from  (5.P)  that  H is  again  given  bv  (5.8b),  but  H,  becomes 

t b 


1/2 


tyw)  mi  - £ - (jL>-  * g (,i 


ij  » t 


(5.  il)a> 


I 


Li 


In  terms  of  tlie  planar  solution  In  (5.8a),  we  mav  rewrite  (5.1()a)  as 

,-Jks 

> ks 


3 /l  2 - j ) :>  •)  i_  e~|kS 

lunar  8 ► 2: r K 0 kR  *— 


8 • - 


H.  IQ)  (H.  (Q)  1 . 

b o p I i 

* l\t» 

(5.10b) 

Uie  result  In  (5.10)  Is  most  Interesting;  and.  In  fact,  somewhat  surprising. 
Hie  surface  rav  traveling  In  direction  8 « it/2  (Fig.  5)  is  a a r r line 


(kK^  * •») . However,  due  to  the  finite  curvature  In  the  binormal  direction 


( * R),  i'n  tlie  cylindrical  surface  differs  from  Its  counterpart  on 
a planar  surface  bv  the  additional  term  In  (5.10b).  At  a large  distance 
away  from  the  source  (ks  *•  •x>)  in  tlie  direction  8 * tr/2,  and  for  a fixed  kR, 


we  find  that  II,  on  a planar  surface  and  that  on  a cvllndric.il  surface  are 
b 


given  by,  respectively. 


Iks 


[H,  (Q)  ] y.  \ — 

b planar  ks 


(5.11) 


- )ks 


-|ks 


H,lQ)  v B 


+ A 


ks 


(5.12) 


where  A and  1)  are  constants  Independent  of  s and  R.  rims , for  large  ks, 

H,  on  a cylinder  is  than  that  on  a plane. 

As  a third  limiting  case,  let 

A -*•  •v  (5.1.)) 

which  occurs  when  observation  point  Q Is  In  the  deep  shadow.  Making  use 
of  ( A- b ) through  (A-10),  we  can  derive  from  (5.n): 


H ,.u  v k ££•  0 l,k|'  1-0.88;.  - + 0.51S  + ks 

‘V  1/3  [/' 

1 l52S(kR) U \ks) 


)] 


H (Q)  ^ IQ) 

t ks  b 


(5.14a) 
( 5 . 1 4 b ) 


Therefore,  In  the  deep  shadow,  the  field  is  a stow  wave  and  decavs 
exponent  la 1 1 v along  the  surface  rav. 


U 





6.  MUTUAL  ADMITTANCE  BETWEEN  SLOTS  ON  A CYLINDER 


Return  to  tile  calculation  of  Y , from  (4.2)  for  two  identical 
circumferential  slots  on  a cylinder  (Figure  6a).  To  calculate  H^  , the 
voltage  excitation  in  slot  1 can  be  replaced  by  an  equivalent 
magnetic  surface  current  density  (Figure  7). 

i£,  = £,  xx  *=  zV.  cos  (^y)  , for  (y,z)  in  slot  1 , (6.1) 

which  radiates  in  a completely  filled  cylinder  [16].  In  (6.1),  y = Rd). 
Making  use  of  the  Green's  function  in  (5.6a),  is  calculated  from  the 
superposition  integral 

■*  f j~T  it 

H,  * i dv,dz,  [V,  /— r cos  (—  y )][b  H,  sin  0 + t H cos  0]  (6.2) 

1 J 1 1 1 / ab  a 7 1 b t 

A1 

where  we  have  written  the  source  point  (v,z)  as  (v^.z^).  Making  use 
of  (6.2)  and  the  electric  field  distribution  of  slot  2 in  (4.2),  we 
obtain  the  final  expression  for  Y between  two  identical  circumferential 
slots  on  a cylinder,  namely. 


_2 

12  = ab 


dy  fdZl 


dy.,dz0  [cos  - y ^ ] [ cos  - (v,,  - R4>Q)]g  (1,2) 


(6.3) 


Here  and  (y0,z,)  are  two  typical  points  in  slots  I and  2, 

respectively.  The  Green's  function  g is 


**(1’2)  = *b 


2 2 
sin*"  0 + H cos-  0 


(6.4) 


where  (H,  ,H  ) are  given  in  (5.h)  with 
b t 


(y 2 ~ yt)“  + (z,  - zt>‘ 


(6.5a) 


tan  l[  (z,  - z^/fy.,  - v^)  ] 
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a 


( 6 . 5b) 


f 


In  a very  similar  manner,  the  mutual  admittance  between  two  identical 
axial  slots  (Figure  ob)  can  be  derived.  The  final  result  reads 

Y12  “ fb  j dvldzl  j dy:dz2  [coa  b 2lI(coS  i(z2  ' V,gz(1*2> 

A.  A0 

(.6. 6) 


where  the  Green's  function  g is 


g (1,2)  * H,  cos”  9 + H sin-  0 
zb  t 


(6.7) 


The  two  surface  integrals  in  (6.3)  or  (6.6)  must  be  evaluated  numerically. 
Extensive  numerical  results  are  given  in  [8],  while  some  representative 
examples  are  quoted  below.  All  values  of  Y are  presented  in  (db,  phase 
in  degrees)  format,  where  db  * 20  log^(|Y^,|  in  mho). 

(i)  Agreement  between  GTD  and  exact  modal  solutions.  Under  the  one- 
modal  approximat ion , an  exact  solution  of  Y , on  a cylinder  can  be  found 
in  terms  of  cylindrical  functions  (the  so-called  "exact  modal  solution") 
[1)  - [3],  [31.  Consider  two  identical  circumferential  slots  with 


parameters 


a = 0.9"  , b - 0.4"  , R = 1.991" 

f » 9 GHz  , \ - 1.3123" 


lb. 8a) 
in .8b  1 


For  various  slot  separations,  values  of  Y^  calculated  by  GTD  solutions 

in  16.  ))  and  by  the  exact  modal  solution  are  presented  in  Table  A.  We 

note  that  they  are  in  excellent  agreement. 

lii)  Effect  of  transverse  curvature  term.  As  explained  in  the 

discussion  of  (5.9)  through  (3.12),  the  rav  travelling  along  the 

generator  of  the  cylinder  is  straight.  However,  the  field  H,  on  it 

0 


is  stronger  than  that  on  a ray  travelling  on  a planar  conducting 
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surface.  Such  a dependence  on  the  surface  curvature  in  the  transverse 

direction  of  the  ray  is  most  interesting,  and  can  be  seen  in  Figure  3, 

where  we  plot  the  ratio 

Yp  on  a cylinder  with  radius  R j 
Yp  on  a plane 

as  a function  of  R for  zQ  = 8"  and  - 0.  We  note  that  the  convergence 

rate  of  the  cylindrical  Yp  to  the  planar  Yp  is  not  as  rapid  as  one 

would  normally  expect.  For  example,  at  kR  = 50,  the  cylindrical  Y,  is 

still  about  10  percent  higher  than  the  planar  one.  The  exact  modal 

solution  in  this  figure  is  truncated  at  kR  = 50,  because  beyond  this 

radius,  it  becomes  extremely  slowly  convergent. 

(iii)  Additional  numerical  results  of  Yp  between  two  identical 

slots  on  a cylinder  are  given  in  Figures  9 to  12.  The  normalized  phase 

is  defined  by  the  phase  of  Yp  expC+iks^),  where  s^  is  the  center-to- 

2 221/2 

center  distance  of  the  slots  and  is  equal  to  (z^  + R 
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7.  GTD  GREEN'S  FUNCTION  FOR  A SURFACE  FIELD  ON  A GENERAL  CONVEX  SURFACE 


To  calculate  the  mutual  admittance  between  slots  on  a general 
convex  surface,  we  have  to  generalize  the  GTD  Green's  function  for  the 
cylinder  in  (5.6).  Referring  to  Figure  13,  let  us  consider  a perfectly 
conducting  convex  surface  E , whose  radii  of  curvature  at  any  point  are 
large  in  terms  of  wavelength.  At  a point  Q^,  described  by  position 
vector  r^  on  E,  there  is  a tangential  magnetic  dipole  source  described 
by  a magnetic  current  density 

K(r)  = M<5  (r  - £ ) (7.1) 

where  M is  the  magnetic  dipole  moment  and  lies  in  the  tangent  plane  of  E. 

The  problem  is  to  determine  a high-frequency  asymptotic  solution  of  H at 

a general  point  Q,,  described  by  position  vector  r.,  on  E.  In  other  words, 

— ^ 

the  GTD  Green's  function  for  the  surface  magnetic  field  for  points  r^  and 
r»  is  to  be  found. 

Before  presenting  the  solution,  let  us  introduce  several  definitions 

and  parameters.  According  to  GTD  |13],  [Id],  the  dominant  high-frequency 

contribution  to  H(r.,)  is  the  field  on  the  surface  ray  from  r ^ to  r,,.  The 

surface  ray  is  a geodesic  of  E.  Some  of  the  geometrical  properties  are 

described  by  (Figure  13)  (i)  the  arc  length  s which  is  chosen  such  that 

s - 0 at  the  source  point  r^  and  s » s at  the  observation  point  r,,;  (ii)  the 

' 

tangent,  normal,  and  binormal,  denoted  bv  (t  , -n  , -b  ) at  r where  n « 1,2; 

n n n n 

and  (iii)  its  two  radii  of  curvature  R (.s),  and  R.  (s)  of  E at  point  s in 

t b 

the  directions  of  tangent  and  binormal,  respectively.  (On  a general 
convex  surface,  both  radii  are  nonnegative.) 

IS 


From  Che  above  parameters,  we  may  calculate  the  following  quantities 
that  are  needed  for  the  solution  of  the  Green's  function: 

(i)  The  large  parameter  in  our  asymptotic  expansion  of  the  Green's 
function  is 


m(s)  - [7  kRt  (s)  ] 1//3 

which  is  a function  of  position  along  the  ray  from  r^  to  r„ 

(ii)  A distance  parameter  from  r^  to  r,  is  defined  by 


(7.2) 


ds 


r ^ 2m“(s) 


(7.3) 


For  the  special  case  when  R^  is  not  a function  of  s (constant  ray  curvature), 

2 

4 is  reduced  to  (ks/2m“),  a well-known  parameter  introduced  first  by 
Fock  [17]. 

(iii)  The  ray  curvatures  at  the  source  and  observation  points  enter 
in  a parameter  defined  by 


[ iiS 1 !/- 

2m(0)  m(s) 


(7.4) 


which  is  positive  real  for  a convex  surface,  and  is  reduced  to  unity  for 
the  special  case  of  a constant  ray  curvature. 

(iv)  Consider  a small  pencil  of  surface  rays  originating  from  r^  and 
propagating  toward  r0  (Figure  13).  The  angle  extended  by  the  pencil  at 
r^  is  di|i^,  and  that  at  r?  is  d^.,.  The  divergence  factor  DF  of  the  pencil 


is  defined  bv 


DF 


sdO/, 


?d^. 


1/2 


(7.5) 


where  p is  the  caustic  distance  of  the  wavefront  at  r,  and  is  always 
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positive.  For  example,  if  E is  a sphere  and  r,  is  the  north  pole,  DF  at 


point  r,  =*  (r,9,J>)  is 


DF 


8 ,1/2 


Lsin  8 

which  varies  from  one  at  the  north  pole  (9  = 0)  to  infinity  at  the  south 
pole  (9  = t)  as  r,,  moves  along  a great  circle. 

(v)  The  "mean"  radii  of  curvature  between  r^  and  r7  are  defined  by 

,1/2 


Rc  = [Rc (0)  Rc (s) ] 


(7.6a) 


,1/2 


R^  = [Rb(0)  1^(3)]"'“  . (7.6b) 

Throughout  this  work,  we  always  assume  that  I is  a smooth  surface  with 
a slowly  varying  curvature.  Then  (R^.R^)  represents  a sort  of  average 
value  of  radii  of  curvature  along  the  ray. 

Return  to  the  electromagnetic  problem  in  Figure  13.  We  assume 
that  m(s)  is  large  and  is  slowly  varying  for  all  s in  the  range 
0 < s < s.  Then  an  approximate  asymptotic  solution  for  the  surface 
magnetic  field  at  r0  due  to  the  dipole  source  in  (7.1)  is  given  by 

H(r,)  = M • (b^,^  + tLt?Hc)(DF)  (7.7a) 

where 


^ = G(s)  { (1  - t^)  rv(£)  - (^-r  T3u(£)  + j(/2kRt)' 


73 


• [tv’(^)  + (Rc/Rb)  t3u’(S)]} 


(7.7b) 


G(3)(liM'V(C)  + (1  ' ki}  t3u(^}  + j(/2kRt)~2/3  I3u'(0]  (7.7 c) 

-jks 


C(s) 


k Yq  e 


ks 


. YQ  = (120n ) 


-1 


(7 . 7d) 
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The  Fock  functions  u and  v and  their  derivatives  u'  and  v'  are  described 


in  the  Appendix.  Several  remarks  about  the  solution  in  (7.7)  are  in  order. 

(i)  It  is  derived  in  an  approximate  manner  from  the  classical 
work  of  Fock  [17)  and  the  recipe  of  GTD,  as  detailed  in  [10].  All 
traditional  GTD  solutions  depend  on  R^,  not  R^.  In  (7.7b),  the  term 
containing  (R^/R^)  was  introduced  through  an  Ansatz  suggested  in  [9]. 

Because  of  the  fact  that  u'  decays  faster  than  v' , this  term  is  important 
only  if  R^  is  very  large  and  R^  is  finite.  An  example  occurs  in  the 
axial  propagation  along  a cylinder,  where  R^  -*■  00  and  R^  is  equal  to  the 
radius  of  the  cylinder.  For  this  particular  example,  it  is  onlv  with  this 
additional  term  that  (7.7b)  agrees  with  the  rigorous  asymptotic  solution 
(derived  recently  by  J.  Boersma  in  an  unpublished  note).  Thus,  the 
Ansatz  in  [9]  is  at  least  partially  verified. 

(ii)  For  the  special  case  that  Z is  a planar  surface  (R  * R^  -*■  °°) , 
(7.7)  recovers  the  known  exact  solution  given  in  (5.8).  When  Z is  a 
cylinder,  (7.7)  is  reduced  to  (5.6). 

(iii)  The  solution  is  valid  for  any  combination  of  r^  and  r,.  In  the 
penumbra  region  (r,  is  close  to  r^  and  Z <K-  1) , (~.7)  gives  approximately  the 
known  planar  solution.  In  the  deep  shadow  (,£;  » 1),  the  residue  series 


representation  of  the  Fock  functions  can  be  used,  and  (7.7)  is  identified  as 


the  creeping-wave  contribution. 

(iv)  Except  for  the  very  simple  surfaces  such  as  a cylinder,  cone  or 
sphere,  no  explicit  parameter  equations  can  be  found  for  the  geodesics  (IS). 


Thus,  for  a general  surface,  one  may  have  to  rely  on  numerical 
for  determining  the  geodesics  and  the  divergent  factor. 


techniques 


3.  GREEN'S  FUNCTION  OF  A CONE 


Let  us  appLy  the  formula  (7.7)  to  the  field  on  an  infinite  cone. 


described  by  the  equations  (Figure  14a) 

x = r sin  0^  cos  4>  , y ■ r sin  0^  sin  ? 


, z = r cos  Q, 


(8.1) 


where  0^  is  the  half-cone  angle  (0  < 0^  < tt/2).  Since  the  cone  is  a 
developable  surface,  the  r3vs  (geodesics)  on  a developed  cone  (Figure  14b) 
are  straight  lines  [18].  Due  to  the  source  at  r^  = (r^,0Q,4>^),  the  main 
contribution  of  the  field  at  r,  = (r,,0~,4O  comes  from  the  shortest  ray 


described  bv 


r^  sin  * r„,  sin 


(8.2) 


As  the  ray  propagates  away  from  the  source  point  r^,  it  reaches  the  highest 
altitude  at  M where  = tt/2.  After  M,  the  ray  travels  downward  away  from 
the  cone  tip.  The  various  parameters  defined  in  Section  7 can  be  simply 
calculated  from  the  cone  geometry,  and  expressed  in  terms  of  coordinates 
(r^,^)  and  (r, ,2>,).  The  arclength  is 


s - { r“  + r"  - 2r^r , cos  [ (^  - 4>0)  sin 


(8.3) 


The  angle  2^  at  r^  is 


_ 1 o , 

H = sin  {—  sin  [ ($  - 4>  ) sin  0 ]} 

1 s 2 10 


(8.4) 
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We  choose  jil  I < tt/2  if  r~  < s~  + r^,  and  |S2  J > tt/2  if  otherwise.  The 


other  parameters  are 


) [*2  * + (1>t  - D^)  sin 


(S.  5) 


»r  r,  tan  0 

l a — = — 

t sin  2^  sin  il,  * 


- _ ‘V;  ta"  -0 

d cos  2^  cos  iT, 


(S.o) 


(8.7) 


1 , , .1/3  | , 2/  3 - 

(-  kr ^ sm  >2,  sin  6Q)  '4>.,  - 5^  cos 


iks/C)  l/:(2k"r1r  ,)"l/6  (sin  fl  sin  .2,  cos  ;2  ) 


1/3 


(S . 8) 


DF  - 1 • (.8.9) 

When  Che  above  parameters  in  (3.3)  through  (.8.9)  are  substituted  into  (7.7). 
we  obtain  an  approximate  solution  for  the  surface  field  on  a cone  due  to 
a direct  surface  ray  contr ibution.  Let  us  consider  a special  observation 
point  r,  such  that 


and  fl,  are  not  close  to  zero 


(8.10) 


After  making  use  of  the  residue  series  representation  for  the  Fock  functions 
(Appendix)  and  keeping  only  the  leading  terms,  then  the  two  components 
of  the  field  in  (7.7)  are  reduced  to 


* 

k“(sin  .2  sin  .2,  cct  9 ) 

% x k T7: T7^—  «P  [-0.S85 


, 2 .1/6  „ . 1/2 
15..3(k  r^r,)  (its) 


H„  -v  0[  (ks) 


-3/2. 


i (jT  + 0.  S 1£1  + ks)  ] (8. 1 la) 


(8.11b) 


which  agree  with  the  rigorous  asymptotic  solutions  given  in  Equations  (50) 

and  (53)  of  [19],  (In  making  the  comparison,  note  the  corresponding 

notations  used  in  [1]  and  here:  -i  •*  i,  A ■*  9,,  L.  -*■  s,  r - r.  , 

c 0 l ' l 

r^  -*  r,,  -*■  t/2  - :2^,  and  -*■  t'  . .)  \s'e  emphasise  that  the  result 

in  (8.  in  or  that  in  [1°]  is  valid  only  under  the  conditions  in  (8.10). 

For  an  arbitrarily  located  observation  point,  (7.')  should  be  used. 

Two  final  remarks  about  the  formula  in  (7.7)  are  in  order.  (i)  For 
a given  source  and  observation  point,  there  are  infinitely  manv  rnvs 
(geodesics)  passing  through  them.  The  contribution  from  each  rav  can  be 
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calculat  from  (7.7),  and  Che  final  field  solution  is  the  superposition 
of  all  r contributions.  In  most  practical  problems  (all  the  numerical 
computat  ns  presented  in  this  paper),  only  the  ray  with  the  shortest 
arclengt  gives  the  significant  contribution  to  the  field  solution,  whereas 
all  othe  rays  may  be  ignored.  (ii)  Depending  on  the  polarization  and 
the  dist  ces  of  the  source  and  observation  points  from  the  cone  tip,  there 
may  be  a ther  significant  contribution  to  the  field  from  the  diffraction 
at  the  t . In  such  a case,  the  total  field  at  any  point  contains  two 
dominant  ontributions : one  from  the  direct  ray  according  to  formula  (7.7), 
and  one  om  the  tip-diffracted  ray.  More  about  the  latter  will  be  given 
in  Secti  9. 
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9.  MUTUAL  ADMITTANCE  BETWEEN  SLOTS  ON  A CONE 

On  the  surface  of  a cone,  let  us  consider  two  arbitrarily  oriented 
slots.  Under  the  assumption  that  the  dimensions  of  the  slots  are  relatively 
small  compared  with  the  radii  of  curvature  of  the  cone  surface,  the  shapes 
of  slots  are  taken  to  be  rectangular  on  a developed  cone. 

Referring  to  Figure  15,  we  describe  the  dimensions  and  the  positions 
of  the  two  slots  by  (a^.b^)  and  [cn>(n  - lH^.ui^l  , n =*  1,2.  Thus,  the 
radial  separation  of  the  two  slots  is  (c0  - c^)  and  the  angular  separation 
is  4>rv  The  angle  u>  measures  the  deviation  of  the  longitudinal  direction 
of  slot  n from  the  radial  direction  of  the  cone.  As  usual,  we  assume  that 
the  slots  are  thin,  and  that  their  lengths  are  roughly  a half-wavelength. 
Then  the  aperture  field  in  each  slot  can  be  adequatelv  approximated  bv 
a simple  cosine  distribution,  which  is  the  "one-mode"  approximation 
described  in  Section  4.  Y^.,  has  two  dominant  high-frequency  contributions: 

one  from  the  direct  rays  going  from  slot  1 to  slot  2,  and  the  other  from 
the  rays  diffracted  at  the  tip  of  the  cone,  viz., 


Y ^ Yd  + YC 
12  1 12  12 


(9.1) 


The  first  term  Y , may  be  explicitly  written  as 
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where 


(cos  y . ) ( cos  y ) 3(1,2) 

al  1 a2 


g(l,2)  * cos  cos  (0^  + H sin  sin 


(9.2a) 


(9.2b) 


Coordinates  (y  ,•  2 ,)  hero  have  their  origin  at  the  center  of  slot  2 (Figure  13), 
not  at  the  center  of  slot  i as  in  the  cylinder  case  (Figure  7). 
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The  Green's  function  components  (H  ,H  ) are  Riven  in  (7.7),  and  angles 

(o),,cjj,)  are  shown  in  Figure  15.  In  evaluating  the  Integrals  in  (9.7a), 
3 -* 

for  two  given  points  (VpZ  ) and  (v,,z.,),  we  must  calculate  some 
geometrical  parameters  appearing  in  H and  Hp  Those  calculations  lead 
to  the  following  results 


2 2 2 fl  2 1/2 

[c“  + v'  + z~  - 2c  /y  + z“  cos  (u)  - a)  . . ) ] 
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(9.3a) 


>(>  a (sin  t),,)  ^ sin  * lt  v“  + z~  r ! sin  (u)  - u>  , , ) ] (9.3b) 
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oJ  , « tan  (2  /v  ) 
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10  - + (ir/2)  - u - sin  + (n  - 1 )<{»,.  sin  0- 

n+2  n n n 0 0 0 


(9. 3c) 
(l).  3d) 


where  n « 1 and  2.  We  evaluate  the  integrals  in  (9.2a)  numerical Iv  with 
the  aid  of  a computer. 

Next,  let  us  consider  ij,,  the  part  of  mutual  admittance  due  to  the 
ravs  diffracted  at  the  cone  tip.  We  approximate  it  by 


Yp  s T sin  iOj  sin  .0 . 


(9.4) 


where  T is  derived  in  [2]  and  is  given  by 

(a.  b^a  ,b0)  tan  “ sin  (kb,/2)  sin  (kb,/2) 


30it  c^c.,  sin  0(. 


I— 


(kb  /2)ikb,/2) 


• exp  J (7-  - kc  - kc.,)  . (9.5) 

4 L Z. 

Here  cj  is  the  zeroth-order  tip  diffraction  coefficient  and  is  a function 
of  the  half-cone  angle  0 . A numerical  table  of  for  several  typical 

values  of  ‘3  1 is  given  in  [2].  We  luive  fitted  those  values  by  a simple 


2n 


expres  on , viz.. 


where 


J j " A exp  ( JB) 


(9.0) 


ll 

\ 

ti 


A - i.30579~l  - 1.755  + 2.7729()  - 1.4599“ 

B » 2.7195  + 1.460890  - 1.12959“  + 0.65669* 

Both  9 and  B are  In  radians.  It  has  been  checked  that  the  numerical 
values  of  o calculated  from  (9.6)  are  in  excellent  agreement  with 
those  tabulated  in  [2]. 

The  final  solutions  for  Yp  (total  mutual  admittance)  and  Y^,  (partial 
mutual  admittance)  are  given  in  (9.1),  (9.2),  and  (9.4).  For  a given 
geometry  of  the  slots  and  cone,  the  two  surface  integrals  in  (9.2a)  are 
evaluated  numerically  bv  choosing  an  integration  grid  roughly  equal  to 
0.05\  v 0.05\.  Unless  specified  otherwise,  all  numerical  computations 
are  based  on  two  identical  slots  with  slot  length  * 0 . 5 \ and  width  = 0.3\. 

a.  "Equivalent"  cylinder.  it  has  been  conjectured  in  [2]  that,  in 
calculating  Y^,  (the  contribution  from  the  direct  ravs)  approximately , 
the  cone  mav  be  replaced  bv  an  "equivalent"  cylinder  with  radius 

R « :j(Cj  + c ,)  sin  0 . (J.7) 

This  conjecture  has  been  quantitatively  checked  out  in  [8].  The 
conclusion  is  that  the  "equivalent"  cylinder  gives  a good  approximation 
for  a small-angled  cone,  e.g.,  0^  - 15°.  However,  the  error  in  y'  , 
calculated  from  the  "equivalent"  cvlinder  can  be  as  large  as  2.5  db  for 
a large-angled  cone  (3  - J0° , for  example). 

b.  Comparison  with  experiments.  A set  of  experimental  data  on 

the  mutual  coupling  between  two  X-band  open-ended  waveguides  (0.°"  ' O.-i") 
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on  a cone  was  reported  in  [2].  As  a function  of  frequency,  measurements 
were  done  on  the  coupling  coefficient  S^.,,  which  is  related  to  Y , 
through  the  formula  in  (3.3).  In  Figures  16  and  17,  three  sets  of  data 
are  presented:  (i)  the  experimental  data;  (ii)  the  theoretical  results 
from  the  present  analysis  in  which  the  calculation  of  Y^.-,  is  based  on 
a cone,  e.g..  Equation  (9.2);  (iii)  the  theoretical  results  from  [2]  in 
which  Y^.,  is  calculated  from  the  exact  modal  solution  of  an  "equivalent" 
cylinder.  Several  observations  can  be  made.  (a)  Both  theoretical  results 
are  in  good  agreement  with  the  experimental  data  (with  the  present  result 
being  slightly  better) . As  explained  in  (9a) , the  "equivalent"  cylinder 
method  works  because  the  cone  angles  (0  t 10°)  are  small.  (b)  The  peaks 
and  valleys  are  caused  by  the  interference  between  Y^?  and  Y^.,,  which  are 
of  comparable  magnitudes  due  to  the  large  angular  separations  (60.3°  and 
30°).  (c)  There  exists  a slight  shift  in  frequency  (Af/f  ~ 3 percent') 

between  the  theoretical  and  experimental  valleys  in  Figure  16.  We 
speculate  that  this  may  be  due  to  a slight  phase  inaccuracy  in  Y^. 

c.  Mutual  admittances  of  circumferential  slots.  In  Figures  13  to  20. 
Y ~ and  for  two  circumf erential  slots  are  displayed  as  functions  of 

angular  separation  and  the  radial  separation  (c^  - c0) . We  note  that 
the  effect  of  Y^0  can  modify  the  curves  of  Y ,,  in  several  different  ways. 
When  the  slots  are  at  the  same  latitude  (Figure  18),  the  direct  coupling 
is  weak.  Thus,  tip  contribution  is  noticeable  even  at  a small  angular 
separation.  As  the  radial  separation  is  increased  (Figure  19),  the  tip 
contribution  is  almost  negligible  for  $ < 65°.  When  the  two  slots  are 

widely  separated  in  the  radial  direction  with  one  slot  near  the  tip 
(Figure  20),  the  tip  contribution  gets  stronger,  the  direct  contribution 
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gets  stronger,  and  the  direct  contribution  becomes  insensitive  to  $ . 
Hence,  the  oscillation  on  the  Y^,  curve  has  a much  larger  period.  In 
fact,  there  is  only  a half-"cvcle"  in  the  range  0 < [>^  < 90°,  and  Y^., 
appears  to  be  shifted  from  Y , by  a fixed  amount. 

d.  Effect  of  slot  orientation  on  mutual  admittance.  Consider  two 
slots  separated  by  1 \ along  the  radial  direction.  The  magnitude  of  Y 
as  functions  of  the  slot  orientation  angles  u>,  and  u),  is  plotted  in 
Figure  21.  As  expected,  the  maximum  value  (-73  db)  occurs  when  both 
slots  are  circumferential  » 90°).  This  value  is  above  1-.  db, 

higher  than  that  when  both  slots  are  radial  (u^  » *0).  The  minimum 

value  (-113  dM  of  Y^0  occurs  when  the  top  slot  is  radial  and  the  bottom 
one  is  c ircumt erent ial . This  resulc  confirms  a common  belief  that  the 
mutual  coupling  between  two  orthogonal  slots  is  generally  negligible. 


10.  SELF-ADMITTANCE  OF  A SLOT  ON  A CYLINDER  OR  A CONE 


The  formula  for  calculating  mutual  admittance  in  (4.2)  can 
be  used  to  calculate  the  self-admittance  Y (or  the  alternative  notation 
Y^),  provided  that  slot  1 coincides  with  slot  2.  However,  in  the  actual 
numerical  evaluation  of  CTD  expressions  such  as  (6.3),  (6.6),  and  (9.2), 
a mathematical  difficulty  arises  as  explained  below. 

Let  us  concentrate  on  (9.2)  with  a^  » b , a,,  = b.,,  c^  » c7,  and 
$ = 0,  which  is  the  direct  ray  contribution  Y3  to  the  self-admittance  of 

a slot  on  a cone.  This  integral  is  divergent , due  to  the  fact  that,  as 
point  approaches  point  Q..,  (Figure  15^, the  Green's  function  in  (9.2b) 
becomes  infinite  as 

g(l,2)  ~ Cs"3  , s - 0 , (10.1a) 

where  s,  defined  in  (8.3),  is  the  distance  between  the  two  points,  and 
the  parameter  C is 

C = (2-3  cos2  fi  ) . (10.1b) 

j2407T“k  1 

It  is  well-known  that  the  singularity  of  cubic  power  is  non-integrable 
with  respect  to  a surface  integral.  This  difficulty  can  be  traced  back 
to  the  derivation  of  the  Green's  function  g.  Strictly  speaking,  g is 
a distribution  and  can  be  written  as 

g ■ Dg  , (10.2) 

where  D is  a second-order  differential  operator  with  respect  to 
coordinates  of  point  Q0,  and  g is  the  Green's  function  of  a vector 
potential  component.  A "legitimate"  expression  corresponding  to  (9.2a) 
should  read 
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(10.3) 


which  is  convergent,  and  Yd  has  a well-defined  finite  value.  However, 
in  writing  (9.2a),  we  have  interchanged  the  differential  operator  D and 
the  second  surface  integration  operator  in  (10.3).  This  interchanging 
is  not  permissible  and,  therefore,  leads  to  the  divergent  integral  in  (9.2a). 

Since  (10.3)  contains  a differential  operator  and  is  not  suitable 
for  numerical  evaluation,  we  prefer  to  work  with  (9.2a),  provided  of 
course  that  we  can  extract  the  correct  finite  part  from  the  divergent 
integral.  To  this  end,  we  rewrite  the  Green's  function  in  (9.2a)  as 


S = S0  + gl  ’ (10.4) 

The  first  term  g^  in  (10.4)  is  the  Green's  function  of  an  infinite  ground 
plane,  and  is  given  by  the  well-known  expression,  c.f.  (5.8), 

g = G(s)  [cos-  f2  + r^— ( 1 - r^~)  (2  - 3 cos-  0, ) I . (10.5) 

0 1 K.S  ks  1 

Note  that,  as  s + 0,  has  exactly  the  same  singular  behavior  in  (10.1) 

as  g.  This  is  expected,  because  in  the  sufficiently  small  neighbornood 

of  a point  source,  the  cone  can  be  approximated  by  its  tangent  plane. 

The  second  term  g^'(g^  = g - g^)  in  (10.4)  is  the  difference  between  the 

Green's  function  of  a cone  and  that  of  a plane.  Near  the  source,  it 

C3n  be  shown  from  (9.2a)  and  (10.5)  that 
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C^s 


-3/2 


s - 0 


(10.6a) 


where 


Cx  = (1920Rt)_I  k"1/2  (1  - j)  (2  - 3 sin2.'^)  . (10. rb) 

When  (10.4)  is  substituted  into  (9.2a),  the  self-admittance  Yd  on  a cone 
is  decomposed  into  two  components,  namely. 
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(10.7) 


vJ  ■ v‘0  * Tj 


The  singularity  of  g^  at  the  source  point  specified  in  (10.6)  is 
integrable.  Thus,  there  is  no  difficulty  in  evaluating  Yd  numerically. 

The  first  term  Yd  is  the  admittance  of  a slot  on  a plane.  It  is 
defined  by  (9.2a)  after  replacing  g by  g^  in  (10.5).  From  (5.6d)  and 
(10.5),  we  recognize  the  following  identity: 


g » d+4  -^r)  ■ < 

k~  Ov” 

Note  that  (10.3)  is  in  the  form  of  (10.2).  Substituting  (10.3)  into 


(9.2a)  and  interchanging  integration  and  differentiation  operators,  we 
obtain 


(10.3) 


f f TTV  _ t t 
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(10.9) 


The  integral  in  (10.9)  is  now  convergent,  and  can  be  considered  as  the 
"finite  part"  of  the  divergent  integral  in  (9.2a).  For  numerical 
evaluation,  (10.9)  in  the  space  domain  is  converted  to  that  in  the 
Fourier  transform  domain.  Following  Rhodes  [20],  it  is  simplified  and 
the  final  result  reads 
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, and 
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where  3 = (k-  - oi“ ) “ , Y = (a"  - k“) 

C(ct)  = CJs2  (aa/1l  . (10.11) 

1 - (aa/-a)“ 

This  is  no  difficulty  in  evaluating  (10.10)  numerically. 

In  summary,  the  direct  ray  contribution  of  for  self-admittance 

of  a slot  on  a cone  as  given  in  (9.2a)  is  divergent,  due  to  an  "illegal" 

interchange  of  integration  and  differentiation  operators  in  the  derivation 

process.  The  (correct)  finite  part  of  the  divergent  integral  is  given 

in  (10.7),  where  Y^  is  given  in  (10.10),  and  Y^  in  (9.2a)  after  replacing 

g by  g^.  The  same  difficulty  arises  in  the  case  of  a cylinder,  and  it 

is  treated  by  the  same  procedure  as  in  the  case  of  a cone. 

We  have  derived  the  self-admittance  Y of  a slot  on  three  types  of 

surfaces:  (i)  For  an  infinite  plane,  the  final  solution  Y = Y^  is  given 

in  (10.10).  (ii)  For  an  infinite  cylinder,  Y = Y^  is  given  in  (10.7). 

For  a circumferential  slot,  Y^  is  given  by  the  integral  in  (6.3)  after 

replacing  g by  (g  - g ),  and  recognizing  that  = (t/2)  - 0 and 
4)  u 1 

A^  = A0.  For  an  axial  slot,  Y^  is  given  by  the  integral  in  (6.6)  after 
replacing  g^  by  (g^  - g^),  and  = 8 and  A^  = A,.  (iii)  For  an  infinite 
cone,  Y has  two  contribut ions  as  described  in  (9.1).  If  one  realized  that 
the  two  slots  in  Figure  15  are  identical,  and  occupy  the  same  position 
on  the  cone,  YC  is  given  in  (9.4)  and  Yd  in  (10.7).  To  calculate  Y^ , we 
use  (9.2a)  after  replacing  g by  g^,  where  g^  can  be  gathered  from  (10.4), 
(9.2b),  and  (10.5).  Numerical  results  of  Y on  the  above  three  surfaces 
are  presented  below. 

Slot  on  a plane.  As  a function  of  slot  length  a,  we  plot  (a/'2b)Y 
in  Figure  22  for  three  different  values  of  slot  width  b.  Those  curves 
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are  practically  linear,  and  can  be  described  for  0.4  < (a/A)  < 0.6  to 


a good  accuracy  by 


[1.029  + j 0 . 596  ] + (3.75  + jB)(f-  - 0.5)}  millimho 

2L  A 


(10.12) 


where  3 = 44,  33,  and  21  for  b = 0.0001X,  0.001X,  and  0.01A,  respectively. 

Relation  to  dipole  impedance.  As  discussed  in  Section  4,  there  is  an 
alternative  definition  for  the  (self-  or  mutual)  admittance  of  a slot. 

This  alternative  self-admittance  Y is  related  to  the  present  Y by  (4.4)  or 

Y = (a/2b)Y  . (10.13) 

From  the  duality  relation  in  Maxwell's  equations,  it  can  be  shown 
(p.  519  of  [11])  that  for  the  special  case  a = X/2, 

Z = -^(120Tr)2(2Y)  , (10.14) 

when  Z is  the  input  impedance  of  a centrally  fed,  half-wavelength  dipole 
radiating  in  the  free  space  (not  in  a half-space  as  in  the  case  of  a 
waveguide-fed  slot).  From  (10.12)  through  (10.14),  we  find  that  for  a 


half-wavelength  dipole, 

Z = 73.12  + j 42.36  ohm 
which  agrees  with  the  results  in  [11],  [20]. 


(10.15) 


a.  Slot  on  a cylinder.  Consider  a circumferential  slot  of  dimension 
0.9"  x 0.4"  on  an  infinitely  long  cylinder  whose  radius  is  3.8". 

Figure  23  shows  Y calculated  by  the  present  GTD  solution  and  that  by  the 
exact  modal  series  solution  in  [2].  These  two  solutions  are  in  agreement 
within  0.52  in  magnitude  and  one  degree  in  phase.  Mote  that,  under  the 
"one-mode  approximation,"  the  modal  series  solution  [2]  is  exact . It  is 
amazing  chat  the  present  GTD  solution  gives  such  an  accurate  result  for 
kR  ~ 18. 
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b.  Slot  on  a cone:  variation  with  radial  distance.  In  all  of  the 
following  cone  calculations , the  slot  has  the  dimensions  of  0.5A  x 0.05X, 
except  when  stated  otherwise.  In  Figures  24  and  25,  the  slot  is 
circumferentially  oriented  on  a cone  with  9^  = 30°,  and  the  variation  of 
Y with  the  radial  distance  c is  presented.  We  observe  two  effects: 

(i)  As  c is  increased,  the  radius  R * c sin  9^  of  the  "equivalent"  cylinder 
becomes  larger  and  larger.  The  magnitude  of  Y decreases  and  approaches 
the  asymptotic  value  of  the  slot  on  a plane.  (ii)  At  c = 2\,  the  tip 
contribution  | Y C | is  less  than  1%  of  the  |y|,  and  this  contribution 
diminishes  as  c increases. 

I 

c.  Slot  on  a cone:  variations  with  cone  angle.  As  9^  is  increased, 
the  cone  surface  becomes  flatter.  Therefore,  Y in  Figure  26  approaches 
its  value  on  a plana. 

d.  Slot  on  a cone:  variations  with  slot  length.  It  is  interesting 
to  observe  from  Figure  27  that  the  minimum  values  of  |y|  for  both  the 

cone  and  plane  cases  occur  roughly  at. a = 0.45A,  not  at  the  resonant  length 
a = 0.5A. 

e.  Slot  on  a cone:  variation  with  slot  orientation  angle.  Figure  28 
shows  that  there  is  about  a 10%  increase  in  Yj'as  varies  from  0 
(radial  slot)  to  t/2  (circumferential  slot). 
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11.  ACTIVE  ELEMENT  PATTERN 


As  explained  in  Sections  2 and  3,  there  are  two  types  of  active 
element  patterns:  P^(6,->)  under  the  short-circuited  condition  in 
Figure  3,  and  Q^(S,4>)  under  the  match-loaded  condition  in  Figure  4. 

Within  the  one-mode  approximation  (Section  4),  we  may  set  l = 0 in 
Figure  3.  Then  P^  becomes  simply  the  radiation  pattern  of  a single  slot 
on  a completely  filled  conducting  surface  (as  if  all  other  slots  were 
absent).  Its  calculation  is  thus  greatly  simplified.  Once  P^,  P^,...^ 
in  an  array  are  known,  the  other  active  element  pattern  can  be  determined 
from  the  matrix  relation  in  (2.10). 

Within  the  one-mode  approximation,  P^(8,J>)  can  be  calculated  by  ray 

-► 

techniques  as  follows:  The  aperture  distribution  E^  of  slot  1 in 
Figure  29a  is  first  replaced  by  an  equivalent  magnetic  surface  current 
s E^  x z in  Figure  29b,  which  radiates  over  a completely  filled 
conducting  surface  Z.  By  using  the  superposition  principle,  the  problem 
of  calculating  pattern  P^(9,$)  due  to  becomes  the  determination  of 
the  Green's  function  due  to  a magnetic  dipole  M on  I (Figure  29c).  Depending 
on  the  location  of  the  observation  point,  the  GTD  calculation  of  the 
Green's  function  is  carried  out  as  follows: 

(i)  Observation  point  A is  in  the  lit  region.  The  field  at  A is 
that  on  the  direct  ray  MA  from  the  point  source.  The  only  effect  of  the 
conducting  surface  Z is  to  image  the  S'Ource,  thus  doubling  its  strength. 

(ii)  Observation  point  B is  in  the  shadow  region.  The  ray  path 
from  the  source  to  3 follows  a geodesic  MD  of  the  surface  until  a point 
of  detachment  D.  From  there  on,  the  ray  continues  to  B along  a straight 


line  DB  tangent  to  Z without  abruptly  changing  direction. 


(iii)  The  lit  and  shadow  regions  are  separated  by  the  tangent  plane 
MT  of  surface  E.  A small  neighborhood  of  MT  is  called  the  transition 
regions,  through  which  the  fields  calculated  from  (i)  3nd  (ii)  should 
bend  smoothly  into  each  other. 

Various  formulas,  which  exist  for  calculating  fields  on  the  rays 
described  above  [14],  [21],  [22],  have  been  applied  to  calculate 
of  a slot  on  a cylinder  [6],  or  a cone  [4],  [22],  Because  of  the  facts 
that  (a)  calculation  steps  are  already  available  in  book  form  [14],  and 
(b)  there  are  still  several  unsettled  (unanswered)  questions  concerning 
the  transition  functions,  torsion  effect,  numerical  accuracy,  etc.,  we 
refer  interested  readers  to  the  literature  for  further  details  concerning 
the  calculation  of  active  element  patterns  by  GTD. 


APPENDIX 


FOCK  FUNCTIONS 


In  this  appendix  we  define  and  list  some  useful  formulas  of  the  functions 
w^it),  w,(t),  v(0,  u(0,  and  v^U).  These  functions  are  commonly  known  as 
Fock  functions. 

(i)  Definition:  For  a complex  t and  a real 
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where  integration  contour  F,  (F,,)  goes  from  ® to  0 along  the  line 

Arg  z = -2i/3  (+2u/3)  and  from  0 to  ® along  the  real  axis.  Because  of 

different  time  conventions,  w (w  ) above  is  equal  to  wn(w.)  defined  in  [17]. 
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(ii)  Residue  scries  representation:  For  real  positive  t, 
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where  C - |t  jexp(-jw/3)  , t'  - |t’  |exp(-jTr/3),  and 
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(iii)  Small  argument  asymptotic  expansion:  For  real  positive  Z,  and  5 •»  0, 
v(£)  'v 
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2.221  x 10_1C5  + ... 


(A-14) 

(A-15) 


(iv)  Numerical  evaluation:  For  E,  _>  £ , the  residue  series  representation 


with  the  first  ten  terms  in  the  summation  can  be  used.  For  £ <_  £q>  the  small 


argument  asymptotic  expansion  with  the  first  five  terms  can  be  used.  It  lias 


been  indicated  in  [7]  that  the  smoothest  crossover  is  obtained  if  Cq  3 0.6. 


In  the  present  study,  we  set  f q = 0.7,  where  the  difference  in  the  two 
representations  is  less  than  0.1%  in  magnitude  and  0.9°  in  phase. 
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Figure  3.  The  radiation  pattern  obtained  by  V = 1,  V = 0 
for  n f 1 is  called  the  short-circuited  active 
element  pattern  P^. 
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between  two  axial  slots  as  a function  of 
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Figure  20.  Mutual  admittance  Y..-,  between  two  circumferential  slots 
on  a cone  calculated  from  GTD  solution. 
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